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ABSTRACT 


This  study  addressed  itself  to  the  development  of  a  Stress/Strength 
Interference  Theory  in  form  of  a  practical  engineering  tool,  to  be  used  for 
designing  and  quantitatively  predicting  the  reliability  of  mechanical  parts 
and  components  subjected  to  mechanical  loading. 

In  our  past  investigation  ("Reliability  Prediction-Mechanical 
Stress/Strength  Interference  Theory"  by  Charles  Lipson,  Narendra  J.  Sheth  and 
Ralph  L.  Disney;  Final  Report  RADC-TR-66-710,  March  1967,  referred  to  in  the 
present  report  as  Reference  1)  the  Interference  Theory  was  developed  for  the 
ferrous  materials.  In  the  present  study,  the  work  was  extended  to  non-ferrous 
materials. 

Early  practices  in  stress-strength  relationship  dealt  almost  entire¬ 
ly  along  the  lines  of  factors  of  safety.  Utilization  of  such  factors  is 
justified  when  they  are  based  on  considerable  experience  with  parts  not  too 
different  from  the  one  under  consideration.  ^However,  when  substantial  changes 
in  the  geometry,  the  processing,  or  the  function  of  the  part  are  contemplated, 
a  major  error  may  result  if  the  old  factors  of  safety  are  projected  to  the  new 
set  of  conditions. 

In  the  present  investigation,  an  approach  was  used  which  attempted 
to  recognize  the  above  limitations.  Instead  of  an  indiscriminate  grouping  of 
all  the  variables  affecting  stress  and  strength  (generally  fatigue  strength) 
into  one  index  (factor  of  safety)  these  variables  were  individually  recognised. 
The  principal  variable  is  the  scatter  in  the  stresses  imposed  on  the  part  and 
*  in  the  strength  of  the  material  resisting  these  stresses. 

The  prevailing  practice  is  to  use  the  mean  values  of  the  calculated 
stress  and  strength,  ignoring  the  natural  scatter  that  each  may  possess.  How¬ 
ever,  the  variability  in  these  two  factors  results  in  a  statistical  distribu¬ 
tion  of  stress  and  strength.  When  these  two  distributions  Interfere,  that  is 
when  stress  becomes  higher  than  strength,  failure  results.  Means  of  expres¬ 
sing  these  distributions,  in  a  practical  engineering  sense,  and  means  of  cal¬ 
culating  the  resulting  interferences,  represent  the  heart  of  the  present  study 

In  our  past  investigation^)  on  ferrous  materials,  the  problem  of 
determining  the  scatter  in  the  fatigue  strength  and, subsequently.  Its  distri¬ 
bution  function  was  resolved  by  a  graphical  method.  Among  the  several  dis¬ 
tribution  functions  that  were  considered,  only  the  Weibull  distribution  was 
used  because  of  its  wide  use  and  the  difficulty  of  handling  other  distribu¬ 
tions  by  a  graphical  method. 

In  the  present  investigation  on  the  non-ferrous  materials,  a  com¬ 
puter  epproach  was  developed  which  mede  the  determination  of  the  distribution 
functions  of  fatigue  strength  considerably  less  difficult  (Section  6).  Con¬ 
sequently,  several  distributions  were  tried:  Weibull,  Largest  Extrema  Value, 
Smallest  Extreme  Value,  Logistic  and  Normal.  The  reasons  for  the  selection 
of  these  particular  distributions  are  given  in  Section  6.3.  It  was  later 
found  that  the  Logistic  distribution  had  similar  characteristics  to  the  Normal 


ill 


distribution  and  therefore  it  was  not  followed  any  further. 

Mechanical  Properties  Data  Center  in  Traverse  City,  Michigan,  was 
found  to  be  a  very  useful  source  of  information  for  the  fatigue  strength  data 
for  various  non-ferrous  materials  under  various  conditions  such  as  the  type  of 
loading, surface  finish,  stress  concentration,  heat  treatment,  temperature, 
environment,  etc.  Considerable  fatigue  data  thus  obtained  were  then  organized 
and  systematized  according  to  materials  and  conditions. 

A  computer  program  was  developed  so  that  when  these  data  (raw  S-N 
data)  were  fed  into  the  program,  the  computer  printed  out  the  degree  of  fit 
and  the  parameters  for  all  the  distribution  functions  tried  for  any  given  set 
of  fatigue  strength  data. 

In  the  first  part  of  this  program,  the  scatter  in  life  was  converted 
into  the  scatter  in  strength.  In  the  second  part,  the  distribution  functions 
(Welbull,  Normal,  Logistic,  LEV  and  SEV)  were  fitted  to  these  strength  data. 
The  best  fitting  distributions  (those  with  the  highest  degree  of  fit)  and 
their  parameters  were  then  chosen,  and  these  are  presented  in  a  tabular  form 
in  Appendix  1  for  various  non-ferrous  materials  and  under  various  conditions. 

In  order  to  present  these  data  in  a  pictorial, and  thus  more  easily 
digestible  form,  the  representative  parameters  were  then  plotted  on  log-log 
scale  against  life  for  various  materials  and  conditions  (Section  6.7).  In 
order  to  put  this  Information  on  the  same  consistent  basis,  all  the  data  were 
expressed  in  terms  of  the  same  distribution  (Normal).  The  plots  were  made  by 
expressing  the  mean  (m)  and  standard  deviation  (0)  as  a  function  of  life 
(Section  6.7).  The  reasons  for  the  selection  of  the  Normal  distribution  for 
these  graphs,  the  description  of  the  computer  approach,  and  the  statistical 
tools  used  in  developing  this  approach  are  given  in  Section  6. 

The  problem  of  stress  distribution  (Section  7)  was  investigated  in 
more  detail  than  in  the  previous  study(^/,but  no  additional  information  was 
located.  When  one  speaks  of  "stress  distribution"  he  usually  refers  to  a 
spectrum  of  loading  or  stresses  to  which  a  part  is  subjected.  Indeed,  most 
of  the  available  data  on  the  subject  Is  expressed  in  this  manner.  In  an  en¬ 
gineering  sense,  this  kind  of  a  distribution  means  the  number  of  times  that  a 
given  part  is  subjected  to  a  given  load  or  stress.  In  the  Interference  Theory 
however,  this  is  not  what  is  wanted.  For  consistency  with  the  strength  dis¬ 
tribution,  the  number  of  parts  subjected  to  a  given  stress  is  required  instead. 

In  the  present  investigation,  the  required  stress  distribution  was 
obtained  In  the  same  manner  as  is  our  past  study.  This  was  done  by  convert¬ 
ing  the  stress  spectrum,  which  generally  has  some  mean  stress,  into  a  spec¬ 
trum  with  taro  mean,  with  the  aid  of  the  Goodman  diagram.  This  facilitated 
the  conversion  of  the  resultant  spectrum  into  an  equivalent  stress  based  on 
taro  mean  stress.  The  conversion  was  accomplished  by  means  of  Miner's  rule. 
The  required  stress  distribution  was  than  expressed  in  terms  of  the  equivalent 
stress.  This  distribution  was  than  compared  with  the  strength  distribution  to 
determine  the  probability  of  interference. 


Most  investigators  in  the  past  have  assumed  both  the  stress  and 
strength  distributions  to  be  normal.  In  those  cases  when  they  were  not  normal, 
a  Monte-Carlo  technique  was  employed.  This  involved  a  sophisticated  means  of 
randomly  selecting  a  sample  from  one  distribution  and  comparing  it  with  a 
sample  from  a  different  distribution. 

In  view  of  the  serious  limitations  of  the  Monte-Carlo  technique, a 
method  of  Integrals  was  developed  (Appendix  3)  and  used  in  the  present  study. 
This  method  involved  developing  an  integral  resulting  from  the  Interference 
of  two  distributions  and  evaluating  this  complex  integral.  Numerical  analysis 
was  carried  out  (Appendix  3)  using  an  IBM  360  computer  to  solve  these  inte¬ 
grals. 


In  our  past  investigation  on  ferrous  materials^),  these  Integrals 
and  hence  the  interference  values  were  evaluated  and  tabulated  for  the  fol¬ 
lowing  combinations: 


Stress  Distribution 
Normal 
Weibull 
Normal 


Strength  Distribution 
Normal 
■Weibull 
Weibull 


In  the  present  investigation  on  non-ferrous  materials,  data  analysis 
(Section  6)  revealed  many  cases  where  fatigue  strength  followed  distributions 
other  than  Weibull  or  Normal.  Hence,  it  was  necessary  to  develop  tables  of 
interference  values  for  the  following  additional  combinations: 

Stress  Distribution  Strength  Distribution 

Normal  Weibull  (Extended-see 

Section  4.2) 

Normal  Largest  Extreme  Value 

Normal  Smallest  Extreme  Value 

These  interference  values  (tabulated  in  Appendix  2)  were  found  to  be 
highly  non-linear.  Therefore,  for  the  purpose  of  interpolating  these  value* 
in  a  given  table  or  between  the  tables,  it  is  recommended  that  a  higher  order 
interpolation  than  linear  be  used.  A  suggested  method  of  interpolation, 
with  some  solved  examples  illustrating  this  method,  is  given  in  Section  A-2.1. 

In  order  to  show  the  application  of  The  Interference  Theory  Tech¬ 
nique  developed  in  the  present  study  two  examples  were  solved,  as  described 
in  Section  9. 
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EVALUATION 


1.  This  study  was  concerned  with  the  development  of  a  practical 
engineering  tool  to  he  used  in  predicting  the  reliability  of  mechanical 
parts  fabricated  of  nonferrous  metals  and  subjected  to  dynamic  loading. 

The  tool  was  to  be  compatible  with  the  Stress-Strength  Interference 
Theory  as  reported  in  RADC-TR-66-710  and  be  usable  by  design  engineers 
with  limited  statistical  background. 

2.  The  study  has  provided  a  prediction  technique  which  is  almost  "cook¬ 
book"  in  nature  and  requires  a  minimum  amount  of  computation  except  as 
required  in  the  interpolation  of  the  tabulated  interference  values.  The 
tabular  values  are  highly  nonlinear;  thus,  linear  interpolation  cannot 
be  used.  Procedures  for  higher  order  interpolation  are  given  and 
illustrated.  Within  the  limitation  that  the  material  must  be  a  nonferrous 
metal,  the  reliability  of  any  part  subjected  to  fatigue  loading  can  be 
predicted.  Examples  of  such  parts  would  include  springs,  fasteners,  shafts, 
beams,  torsion  bars,  vehicle  frames  and  suspensions,  forgings  and  castings. 

The  accuracy  of  the  technique  depends  on  the  knowledge  the  user  has  of  the 

sti’ess  distribution,  l.e.,  measured  or  estimated.  *  ’ 

3>  The  results  of  this  study  and  of  the  earlier  study  of  ferrous  metals 
will  be  combined  and  presented  in  the  RADC  Nonelectronic  Reliability 
Notebook. 


DONALD  W.  FULTON 
Reliability  Engineering  Section 
Reliability  Branch 


SECTION  1  INTRODUCTION 


A  complete  cycle  of  reliability  is  made  up  of  four  stages;  1.  Speci¬ 
fication  of  Reliability;  2.  Prediction  of  Reliability;  3.  Verification  of 
Reliability;  4.  Preservation  of  Reliability.  The  heart  of  the  second  stage, 
namely,  the  Prediction  of  Reliability,  is  the  Interference  Theory. 

The  basic  idea  behind  reliability  is  that  a  given  part  has  certain 
physical  properties  which,  if  exceeded,  will  result  in  failure.  The  factor 
which  may  cause  these  properties  to  be  exceeded  is  the  ' tress  imposed  by  the 
operating  conditions.  Thus,  in  prediction  of  reliability  it  is  not  the  stress 
alone  or  the  strength  alone  that  is  the  determining  factor  but  the  combined 
effect  of  the  two. 

Early  practices  in  stress-strength  relationships  dealt  almost  entire¬ 
ly  along  the  lines  of  factors  of  safety.  Once  a  part  was  designed  and  the 
ratio  of  strength  to  stress  was  in  the  range  of  approximately  5  to  10,  it  was 
considered  to  be  safe  for  service.  In  certain  industries  and  in  certain 
applications,  factors  of  safety  as  high  as  20:1  were  employed. 

The  definition  of  the  factors  of  safety  varied  from  user  to  user, 
depending  on  the  sophistication  and  the  complexity  of  the  problem. 

Some  of  these  definitions,  found  in  literature,  are  listed  below: 
Factor  of  Safety 
Factor  of  Safety 
Factor  of  Safety 
Factor  of  Safety 
Factor  of  Safety 
Factor  of  Safety 
Margin  of  Safety 
Design  Factor 
Factor  of  Utilisation 


»  Ultimate  Strength 
Nominal  Stress 

_  Yield  Strength 
Nominal  Stress 

m  Ultimate  Strength 
Actual  Working  Stress 

m  Yield  Strength 
Actual  Working  Stress 

m  Maximum  Safe  Load 
Normal  Load 

m  Computed  Strength 
Computed  Load 

m  Strength-Stress 
Stress 

.  Strength 
Design  Stress 

•  Stress 

Strength 


1 
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Magnitude  of  Variable 
Producing  Failure 

Functional  Reserve  Factor  -  Magnitude  of  Variable 

at  Operating  Conditions 

where  the  variable  could  be  force, 
power,  torque,  material,  surface 
finish,  fillet  radius,  etc. 


Apparently  the  factor  of  safety  was  meant  to  account  for  all  the 
variables  which  were  known  to  affect  the  stress  and  strength  of  the  member. 

The  utilization  of  a  factor  of  safety  of  this  kind  has  justification  only  when 
its  value  is  based  on  considerable  experience,  with  parts  not  too  different 
from  the  one  under  consideration.  However,  when  substantial  changes  in  the 
geometry,  the  processing,  or  the  function  of  the  part  are  contemplated,  a 
major  error  may  result  if  the  old  factor  of  safety  is  projected  to  the  new 
set  of  conditions. 

This  can  be  illustrated  by  the  problem  of  automotive  axle  shafts 
which  were  failing  in  service  in  large  numbers.  These  shafts  had  been 
fabricated  from  a  steel  with  a  tensile  strength  of  240,000  psi,  and  yet,  the 
operating  stresses  as  measured  in  actual  service  were  fnund  to  be  only  13,000 
psi.  This  produced  an  apparent  factor  of  safety  of  240,000/13,000  ■  18.5. 

This  is  obviously  a  fictitious  value,  since  the  shafts  were  failing  in  ser¬ 
vice,  and  the  true  factor  of  safety  was  less  than  one.  The  explanation  lies 
in  the  fact  that  axle  strength  to  be  compared  with  the  13,000  psi  operating 
stress  should  not  have  been  the  ultimate  strength  of  the  material  (240,000 
psi)  but  the  fatigue  strength  corresponding  to  the  surface  finish  of  the 
shafts,  the  mode  of  loading  to  which  the  shaft  was  subjected,  etc.  When  the 
ultimate  strength  was  reduced  by  these  derating  factors,  the  resultant  value 
was  found  to  be  12,000  psi.  This  strength,  when  compared  with  the  13,000 
psi  stress,  produced  the  realistic  factor  of  safety  of  0.9. 

Examples  such  as  this  lead  to  the  next  phase  in  the  relationship 
between  stress  and  strength,  namely,  to  the  concepts  of  a  significant  stress 
and  e  significant  strength.  By  significant  stress  is  meant  the  actual  stress 
Imposed  on  the  part,  and  it  may  Include  the  effect  of  stress  raisers,  magnifi¬ 
cation  due  to  Impact  loading,  residual  stresses,  etc.  By  significant  strength 
(generally,  fatigue  strength)  is  meant  the  actual  strength  of  the  part  in  its 
fabricated  form  under  actual  operating  conditions.  A  rational  approach  to 
significant  strength  still  employs  ultimate  strength  as  the  basis.  However, 
instead  of  an  indiscriminate  grouping  of  all  the  factors  affecting  the  ulti¬ 
mate  strength  into  one  index,  it  attempts  to  evaluate  quantitatively  the 
effect  of  each  individual  factor  pertaining  to  the  part  and  the  conditions 
under  consideration.  The  result  is  a  value  which  is  strictly  applicable  to 
the  part  under  consideration  and  to  the  set  of  loading  conditions  to  which  the 
part  is  subjected  in  service.  The  principal  factors  affecting  strength  and 
which  must  be  considered  in  determining  the  significant  strength  are:  life 
expectancy,  type  of  loading  (axial,  bending,  torsional,  or  a  combination), 
else  effect,  aurface  finish,  surface  treatment,  notch  effects,  mode  of  loading 
(static,  completely  reversed  dynamic,  or  a  combination). 


These  concepts  of  significant  stress  and  significant  strength  repre¬ 
sent  a  major  step  toward  a  more  realistic  prediction  of  suability  and,  as 
such,  they  have  been  included  in  the  present  investigation.  By  themselves, 
however,  they  are  not  sufficient.  This  is  because  the  prevailing  practice  is 
to  use  the  mean  values  of  the  calculated  strength  and  stress,  ignoring  the 
natural  scatter  that  stresses  and  strengths  may  have. 

The  variability  in  these  two  factors  results  in  the  existence  of  a 
statistical  distribution  function  of  stress  and  strength  (See  Figure  2.1)  and 
is  the  heart  of  the  Interference  Theory.  Thus,  for  proper  prediction  of  relia¬ 
bility,  an  estimate  must  be  made  of  both  the  mean  value  and  the  dispersion 
characteristics  of  both  the  strength  and  stress. 

The  strength  of  the  part,  as  all  properties  of  nors-homogeneous 
materials,  varies  from  specimen  to  specimen,  in  view  of  the  variation  in 
hardness,  surface  finish,  degree  of  stress  concentration,  etc.  In  our  past 
investigation  ("Reliability  Prediction  -  Mechanical  Stress/Strength  Interfer¬ 
ence"  by  Charles  Lipson,  Narendra  J.  Sheth,  and  Ralph  L.  Disney,  Final  Report 
RADC-TR-66-710,  March  1967,  referred  to  in  the  present  report  as  Reference  1  ) 
the  problem  of  determining  the  variation  in  strength  and  subsequently  its 
distribution  function  was  resolved  by  a  graphical  method.  Among  the  several 
possible  distribution  functions  considered,  only  the  Weibull  distribution  was 
used  because  of  its  wide  use  and  the  difficulty  of  handling  other  distribu¬ 
tions  by  a  graphical  method.  In  the  present  investigation  a  computer  approach 
was  developed  by  which  the  determination  of  the  distribution  functions  of 
strength  became  considerably  less  difficult.  (This  is  discussed  in  Section  6.) 
Consequently,  several  distribution  functions  were  tried.  These  were:  Weibull, 
Largest  Extreme  Value,  Smallest  Extreme  Value,  Logistic  and  Normal.  The  rea¬ 
sons  for  selection  of  these  particular  distributions  are  given  in  Section  6.3. 
It  was  later  found  that  the  Logistic  Distribution  had  somewhat  the  same  char¬ 
acteristics  as  the  Normal  Distribution  and,  therefore  it  was  not  followed  any 
further.  ”The  computer  printed  the  degree  of  fit  and  the  parameters  for  all 
of  the  above  distributions  for  a  given  set  of  strength  data.  The  best  fitting 
distributions  and  their  parameters  were  then  chosen,  and  they  are  presented  in 
a  tabular  form  in  Appendix  1  for  various  materials  and  conditions  studied. 

The  representative  ones  were  then  plotted  in  order  to  determine  the  effect  of 
various  conditions  on  the  strength  of  various  materials.  In  order  to  put  the 
information  on  the  same  consistent  basis,  all  the  data  were  expressed  in  terms 
of  the  same  distribution  (Normal)  and  the  plots  were  made  by  expressing  the 
means  (u)  and  standard  deviations  (o)  as  a  function  of  life  (See  Section  6). 
The  reasons  for  the  selection  of  the  Normal  Distribution  for  these  graphs  are 
given  in  Section  6. 

As  is  the  case  of  strength,  the  operating  stresses  vary  too.  These 
stresses  vary  from  time  to  time  in  a  particular  part,  from  .part  to  part  in  a 
particular  deaign,  and  from  environment  to  environment.  Therefore,  both  the 
mean  value  and  the  dispersion  characteristics  of  stress  and  atrength  must  be 
determined. 

Once  the  parameters  of  the  atrength  and  streaa  distributions  are 
found,  percent  Interference  and  thus  probability  of  failure  can  be  determined 
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from  the  interference  area  (shaded  area  in  Figure  2.1).  Means  of  computing 
these  interferences  represent  one  of  the  principal  objectives  of  the  present 
Investigation. 
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SECTION  2  INTERFERENCE  THEORY 


Suppose  there  are  two  barrels  containing  slips  of  paper,  each  having  a 
number  printed  on  it.  The  numbers  in  barrel  Y  are  distributed  according  to 
distribution  Y,  as  in  Figure  2.1,  and  the  numbers  in  barrel  X  are  distributed 
according  to  distribution  X.  If,  at  random,  slips  of  paper  from  each  barrel 
are  selected  and  paired,  they  may  be  classified  into  successes  and  failures. 

A  success  is  constituted  by  a  strength  value  exceeding  a  stress  value,  as  for 
example,  when  xj  >  yj_.  Failure  will  occur  if  x2  <  y2  as  shown.  It  will  be 
noted  that,  although  the  shaded  area  is  a  measure  of  interference,  it  is  not 
interference  itself:  a  pair  of  points  X3  and  y3,  although  in  the  shaded  area, 
will  not  produce  failure.  By  continued  pairing  of  stresses  and  strengths  at 
random,  pairs  will  be  found  where  the  stress  will  exceed  the  strength.  By 
continued  experimentation  a  good  estimate  of  the  probability  of  interference 
can  be  found. 


2.1  TWO  NORMAL  DISTRIBUTIONS 


From  an  exhaustive  survey  of  literature  made  during  the  past  investiga¬ 
tion  (1)*  it  was  found  that  most  studies  have  assumed  both  the  stress  and  the 
strength  distributions  to  be  normal.  This  is  a  natural  assumption  to  make  In 
order  to  solve  a  practical  problem,  as  no  work  was  found  dealing  with  an  ana¬ 
lytical  expression  for  the  Interference  of  two  non-normal  distributions. 

When  the  stress  and  strength  distributions  are  assumed  to  be  normal, 
the  probability  of  interference  can  be  determined  from  the  equation(2); 


z  ■ 

where  u  ■  mean  stress 

y 

yx  ■  mean  strength 
2 

ox  -  strength  variance 
Oy  ■  stress  variance 

z  ■  standardized  normal  variate  determined  from 
standard  tables  (See  Table  2.1) 


(2,1) 


Thus,  if  the  average  stress  la  30  ksl  with  a  standard  deviation  of  3tai 
and  the  average  strength  la  50  kal  with  standard  deviation  of  10  ksl,  z  ■  1.91 
and  from  Table  2.1  a,  which  represents  Interference,  is  found  to  be  .0281. 
Thus,  the  percent  interference  (the  probability  of  failure)  Is  2.81%. 


*  Numbers  In  parenthesis  designate  References  listed  on  Page  149. 


Frequency  of  Occurrence 


Stress  (Strength) 


Figure  2.1  Interference  of  Stress  end  Strength  Distributions 
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NORMAL  DISTRIBUTION 


Tabulation  of  the  values  of  a  versus  KQ  for  the 
Standardized  Normal  Curve .  „ 


a 


P(z  > 


Ka> 


•/ 


1 

/2ir 


dz 


.00 

.01 

.02 

•C3 

.04 

.05 

.06 

.07 

.08 

.00 

0.0 

5000 

.4960 

.4920 

.4880 

.4840 

.4801 

.4761 

.4731 

.4681 

.4641 

0.1 

.4602 

.4562 

.4522 

.4483 

.4443 

.4404 

.4364 

.4325 

.4286 

.4247 

o.a 

4207 

.4168 

.4129 

.4090 

.4052 

.4013 

.3074 

.3036 

.3897 

.3850 

0.3 

.3821 

.3783 

.3745 

.3707 

.3669 

.3632 

.3594 

.3557 

.3520 

.3483 

0.4 

3446 

.3400 

.3372 

.3336 

.3300 

.3264 

.3228 

.3192 

.3156 

.3121 

0.5 

.3085 

.3050 

.3015 

.2981 

.2946 

.2912 

.2877 

.2843 

.2810 

.2776 

0.6 

.2743 

.2709 

.2676 

.2643 

.2611 

.2578 

.2546 

.2514 

.2483 

.2451 

0.7 

.2420 

.2389 

.2358 

.2327 

.2296 

.2266 

.2236 

.2206 

.2177 

.2148 

0.8 

.2110 

.2000 

.2061 

2033 

.2005 

.1077 

.1949 

.1922 

.1894 

.1887 

0.0 

.1841 

.1814 

.1788 

.1762 

.1736 

.1711 

.1685 

.1660 

.1635 

.1011 

1.0 

1587 

.1562 

.1539 

.1515 

.1492 

.1469 

.1446 

.1423 

.1401 

.1379 

l.l 

.1357 

.1335 

.1314 

.1292 

.1271 

.1251 

.1230 

.1210 

.1190 

.1170 

1.2 

.1151 

.1131 

.1112 

.1093 

.1075 

.1056 

.1038 

.1020 

.1003 

.0985 

1.3 

.0968 

.0951 

.0934 

.0918 

.0901 

.0885 

.0869 

.0853 

.0838 

.0823 

1.4 

.0808 

.0793 

.0778 

0784 

.0749 

.0735 

.0721 

.0708 

.3694 

.0681 

1.5 

.0668 

.0655 

.0643 

.0630 

.0618 

.0600 

.0594 

.0582 

.0571 

.0550 

l.fi 

.0548 

.0537 

.0520 

.0516 

.0505 

.0495 

.0485 

.0475 

.0465 

.0455 

1.7 

.0446 

.0436 

.0427 

.0418 

.0409 

.0401 

.0392 

.0384 

.0375 

.0367 

1.8 

.0359 

.0351 

.0344 

.0336 

.0329 

.0322 

.0314 

.0307 

.0301 

.0294 

1.0 

.0287 

.0281 

.0274 

.0268 

.0262 

.0250 

.0250 

.0244 

.0239 

.0233 

2.0 

0228 

.0222 

.0217 

.0212 

.0207 

.0202 

.0107 

.0192 

.0188 

0183 

2.1 

.0179 

.0174 

.0170 

.0166 

.0102 

.0153 

.0154 

.0150 

.0146 

.0143 

22 

.0139 

.0136 

.0132 

.0129 

.0125 

.0122 

.0119 

.0116 

.0113 

.0110 

2.3 

.0107 

.0104 

.0102 

.00990 

.00064 

.00939 

.00914 

.00889 

.00866 

.00842 

2.4 

.00820 

.00798 

.00776 

.00755 

.00734 

.00714 

.00695 

.00676 

.00657 

.00639 

2.5 

.0002 1 

.00604 

.00587 

.00570 

.00554 

.00539 

.00523 

.00508 

.00494 

.00480 

2.6 

.00406 

.00453 

.00440 

.00427 

.00415 

.00402 

.00391 

.00379 

.00368 

.00357 

2.7 

.00347 

.00336 

.00328 

.00317 

.00307 

.00298 

.00289 

.00280 

.00272 

.00264 

2.8 

.00256 

.00248 

.00240 

.00233 

.00226 

.00219 

.00212 

.00205 

.00190 

.00193 

2.0 

COI87 

.00181 

.00175 

.00160 

.00164 

.00159 

.00154 

.00140 

.00144 

.00139 

Table  Id  Normal  Distribution^ 
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In  practical  applications  of  the  Interference  Theory  the  following  pro¬ 
blem  arises:  both  distributions  under  consideration  extend  to  plus  and  minus 
infinity.  It  is  apparent,  therefore,  that  any  two  distributions  will  overlap 
and  cause  interference.  This,  of  course,  is  erroneous  because  some  distribu¬ 
tions,  such  as  strength,  must  have  a  finite  lower  bound  of  zero.  In  many 
situations  the  physical  set-up  and  the  sample  size  adjust  for  this  lower  bound. 
For  example,  suppose  a  part  is  designed  so  that  the  mean  of  the  strength  dis¬ 
tribution  is  placed  6  o  away  from  the  mean  of  the  stress  distribution,  both 
distributions  having  standard  deviation  equal  to  a.  From  Equation  2.1  it  is 
found: 


6  o 

Z  * 


4.24 


and  the  probability  of  interference  comes  out  to  be  .00001.  This  means  that 
only  one  part  will  fail  in  100,000  parts  produced.  If  actually  only  50,000 
parts  are  made,  the  physical  problem  has  effectively  truncated  the  distribu¬ 
tions.  However,  the  probability  of  failure  of  one  part  remains  .00001.  This 
means  that, due  to  sample  size,  the  extreme  portions  of  these  distributions  are 
no  longer  important  since  the  sample  size  is  such  that  not  even  one  failure 
can  be  expected. 


For  the  ease  of  application  of  the  Interference  Theory  for  the  case  of 
two  Normal  distributions,  a  graph  was  constructed  as  shown  in  Figure  2.2.  The 
example  solved  previously  through  Equation  (2.1)  now  yields: 


_  50ksi-30k8i  -  6.67  and, 

a  .  3ksi 


a 

max 

o  . 

min 


lOksi 
3  ksi 


3.3 


and  the  percent  interference  (the  probability  of  failure)  comes  out  approx¬ 
imately  yf>  as  before. 

2.2  CONSISTENCY  OF  TWO  DISTRIBUTIONS 


-  In  the  application  of  the  Interference  Theory,  the  following  important 
point  must  be  conaidered:  the  distribution  of  stress  and  the  distribution  of 
strength  must  be  consistent  with  each  other.  In  a  fatigue  teat  or  in  actual 
application  in  service,  a  single  part  has  a  single  fatigue  life  for  a  given 
loading  condition.  Subsequent  tasting  of  additional  parts  under  the  same  load 
will  show  a  scatter  in  life  leading  to  a  life  distribution.  Through  more  ex¬ 
tensive  tasting  a  strength  distribution  for  a  given  life  can  be  obtained,  such 
as  the  distribution  in  Figure  2.1  (the  method  of  obtaining  a  strength  distri¬ 
bution  from  a  life  distribution  is  described  in  Section  6). 


It  follows  then  that,  for  a  consistent  development,  the  stress  distribu¬ 
tion  must  be  of  the  same  nature  as  the  strength  distribution.  That  is,  it 
should  represent  the  plot  of  the  frequency  of  occurrence  of  an  applied  stress 
versus  the  applied  stress.  This  is  not  the  same  as  the  stress  distribution  con¬ 
ventionally  derived  from  the  spectrum  of  loading  acting  on  the  part.  The  con- 


Figur*  2.2  Probability  of  Interference  of  Two  Normal  Distributions. 


version  of  one  to  the  other  must  be  accomplished  through  the  use  of  the  equi¬ 
valent  utress  (Sequ)  and  Miner's  or  Corten-Dolan's  Rules.  This  is  the  method 
used  in  the  present  investigation,  as  described  in  Section  7. 


2.3  NON-NORMAL  DISTRIBUTIONS 

So  far  the  discussion  has  been  limited  to  the  cases  when  both  the  stiess 
and  strength  distributions  can  be  assumed  to  be  normal.  In  cases  when  either 
one  or  both  are  not  normal  the  problem  is  much  more  involved.  For  example, 
the  intersection  of  a  Normal  and  a  Weibull  distribution  produces  a  distribu¬ 
tion  of  an  unknown  origin. 

In  the  past,  problems  such  as  this  were  solved  largely  through  "brute 
force,"  by  a  method  commonly  referred  to  as  the  Monte-Carlo  Technique.  Essen¬ 
tially,  the  Monte-Carlo  Technique  consists  of  a  sophisticated  means  of  random¬ 
ly  selecting  a  sample  from  one  distribution  and  comparing  it  with  a  random 
sample  taken  from  a  different  distribution.  This  is  accomplished  with  the  aid 
of  Tables  of  Random  Numbers.  The  resultant  paired  data  are  plotted  as  a  Cumu¬ 
lative  Distribution  Function  on  Normal,  Weibull,  etc.  probability  papers  and 
the  percent  interference  is  read  from  the  graph. 

2.4  THE  INTEGRAL  METHOD 

In  the  present  investigation  a  Method  of  Integrals  was  used  in  prefer¬ 
ence  to  the  Monte-Carlo  Technique.  This  method  involves  determining,  as  an 
integral,  the  expression  for  the  interference  of  the  two  distributions  under 
consideration  and  establishing  percent  interference  from  this  integral. 

The  advantages  of  the  Integral  Method  are: 

1.  For  some  distributions  the  integrals  have  been  already  tabulated 
and  percent  interference  can  be  read  directly  from  the  table. 

2.  In  those  cases  where  the  integrals  have  not  been  already  tabulated, 
they  can  be  evaluated  by  Numerical  Analysis  as  done  in  the  present 
investigation. 

3.  The  major  shortcoming  of  the  Monte-Carlo  Technique  is  that  it 
requires  a  very  large  sample  size  for  any  accuracy.  This  short¬ 
coming  is  avoided  when  the  Method  of  Integrals  is  used. 

4.  One  of  the  objectives  of  the  present  study  is  to  develop  and 
evaluate  an  analytical  expression  for  Interference  of  any  two 
distributions.  Such  an  expression  is  possible  when  the  Method  of 
Integrals  is  used  but  not  when  the  Monte-Carlo  Technique  is  employed 
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SECTION  3  OBJECTIVES  OF  THE  PRESENT  STUDY 

The  objectives  of  the  study  described  in  this  report  were: 

1.  To  refine  and  to  reduce  to  practice  the  Stress/Strength  Inter¬ 
ference  Theory  technique  for  designing  and  predicting  the  quanti¬ 
tative  reliability  of  mechanical  parts  and  components  fabricated 
from  various  non-ferrous  materials  under  mechanical  loading. 

Maximum  use  was  to  be  made  of  empirical,  practical  engineering 
values  as  well  as  a  sound  theory. 

2.  To  study  the  effect  of  such  factors  as  type  of  loading,  surface 
finish,  temperature,  heat  treatment,  stress  concentration,  surface 
treatment,  manufacturing  processes,  environment,  etc.  on  the 
statistical  distribution  of  fatigue  strength. 

3.  To  determine,  from  the  existing  available  empirical  data,  the 
distribution  of  the  fatigue  strength  under  the  effect  of  each 
of  the  above  factors. 

4.  To  develop  a  computer  approach  [as  contrasted  with  the  graphical 
approach  of  our  past  investigation(l) ]  for  determining  the  sta¬ 
tistical  distribution  function  of  fatigue  strength  at  a  given  life. 
This  method  has  the  advantage  of  high  accuracy  and  time  saving. 

5.  To  develop  the  means  of  synthesizing  the  strength  distribution 
function  when  such  function  is  non-time  variant,  i.e.  infinite  life 
design,  and  when  such  function  is  time  variant,  i.e.  finite  life 

>  design. 

6.  To  develop  an  analytical  expression  of  the  distribution  of  inter¬ 
ference  for  the  general  case  where  the  two  interfering  distribu¬ 
tions  are  different.  In  the  present  investigation  the  distributions 
studied  were: 


Stress 

Strength 

Normal 

Welbull 

Normal 

Largest  Extreme  Value 

Normal 

Smallest  Extreme  Value 
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SECTION  4  STUDY  APPROACH 


4.1  LITERATURE  SEARCH 

In  the  past  investigation  on  ferrous  materials (1) ,  an  exhaustive  liter¬ 
ature  search  was  made  to  determine  the  State  of  the  Art  in  the  field  of  Relia¬ 
bility  Prediction  -  Mechanical  Stress/Strength  Interference.  Specific  topics 
covered  were:  Interference  Theory,  Mathematical  Tools  as  related  to  the  Inter¬ 
ference  Theory,  Monte-Carlo  Technique,  Reliability  Prediction,  Fatigue  of 
Metals,  Cumulative  Damage,  etc.  As  these  topics  were  common  to  the  past  and 
the  present  investigations,  it  was  not  necessary  to  repeat  them  in  the  current 
study.  Therefore,  in  the  present  project,  no  literature  search  was  made. 

4.2  THEORETICAL  ANALYSIS 


One  of  the  objectives  of  this  investigation  was  to  develop  and  evaluate 
an  analytical  expression  for  the  interference  of  two  distributions.  When  the 
two  distributions  are  Normal  the  interference  can  be  simply  expressed  by  a  z- 
distribution,  as  described  in  Section  2.  From  an  extensive  survey  of  litera¬ 
ture  made  in  our  past  investigation,  no  work  was  found  dealing  with  the  analy¬ 
tical  expression  for  interference  when  the  two  interfering  distributions  are 
not  Normal.  The  purpose  of  this  phase  of  the  investigation,  then,  was  to 
develop  such  expressions. 

In  the  past  investigation. the  analytical  expressions  and,  hence,  the 
tables  of  interference  values,  were  developed  by  means  of  the  Method  of  Inte¬ 
grals  as  discussed  in  Section  2.4  of  Reference  1.  These  included: 


Stress  Distribution  Strength  Distribution 

Normal  Normal 

Weibull  Weibull 

Normal  Weibull 

In  the  present  investigation,  data  analysis  (Section  6)  revealed  cases 
where  fatigue  strength  data  followed  distributions  other  than  Weibull  or  Nor¬ 
mal.  Hence  it  was  necessary  to  develop  tables  of  interference  values  for 
cases  where  stress  distribution  was  Normal  and  strength  distribution  was  Lar¬ 
gest  Extreme  Value  and,  also,  stress  Normal  and  strength  Smallest  Extreme 
Value.  This  was  done  by  evaluating  the  complex  integrals  expressing  the  dis¬ 
tribution  of  Interference.  Numerical' analysis  was  carried  out  (Appendix  3) 
using  an  IBM  360  computer  to  solve  these  Integrals.  Tables  were  then  prepared 
for  the  interference  values  as  a  function  of  the  distribution  parameters 
(Appendix  2). 


In  addition,  in  the  present  investigation  it  was  necessary  to  extend 
the  scope  of  study  for  the  case  when  the  stress  distribution  is  Normal  and  the 
strength  distribution  is  Weibull.  In  the  past  investigation  on  ferrous  mate¬ 
rials  vl),  the  study  was  limited  to  the  cases  where  the  stress-strength  inter- 
feme.  P«— Mr,  C  „  _  m.  _  .t„„.ch  ind... 

'  o  stress  index  * 
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was  between  10  and  100.  This  was  because  the  strength  data  on  ferrous  mate¬ 
rials  fell  into  this  range  of  C  values.  In  the  present  investigation  on  the 
strength  data  of  non-ferrous  materials,  C  values  fell  between  1.0  and  10.0  in 
many  cases.  Hence,  it  was  necessary  to  extend  the  interference  tables  t|o  in¬ 
clude  values  of  C  between  1.0  and  10.0.  Thus,  tables  of  interference  vajlues 
were  prepared  to  include  the  combinations: 


Stress  Distribution 
Normal 
Normal 
Normal 

These  tables  are  given  in  Appendix  2. 


Strength  Distribution 
Weibull  (Extended) 
Largest  Extreme  Value 
Smallest  Extreme  Value 


4.3  EMPIRICAL  DATA 


I 


Interference 

and 


The  tables  of  interference  values  are  the  heart  of  the 
Theory  as  applied  to  engineering  practice.  Once  the  stress  distribution 
strength  distribution  parameters  are  known,  the  percent  interference,  an<jl 
the  probability  of  failure,  can  be  read  directly  from  these  tables  (for 
cedure  see  Section  9  and  for  the  tables  see  Appendix  2) . 


ti 


thus 
e  pro- 


Over  250  articles  from  literature  and  other  sources  were  examined I in 
the  past  investigation,  and  practically  no  data  were  found  concerning  the  sta¬ 
tistical  distribution  of  stress.  In  the  present  investigation  some  further 
work  was  done  in  this  area,  but  no  additional  information  was  found.  The  only 
data  located  referred  to  spectrum  of  loads  or  stresses,  which,  as  pointed  out 
in  Section  7,  does  not  represent  the  stress  distribution  required  for  the 


Interference  Theory, 
for  given  dispersion 
found.  The  range  of 
engineering  practice 


Interference  Tables,  therefore,  were  constructed /so  that, 
characteristics  of  stresses,  a  percent  interference  can  be 
these  characteristics  chosen  here  and  corresponding  to 
are  (if  the  stress  distribution  is  Normal) :  t 


.01 


<  £  < 
U 


.10 


To  determine  the  distribution  of  strength  it  was  necessary  to  collect  a 
great  deal  of  data  in  order  to  arrive  at  a  meaningful  distribution.  To  sys¬ 
tematize  the  effort  of  collecting  data,  a  format  was  prepared  which  Included 
the  factors  which  are  known  to  affect  the  final  distribution  bf  strength.  An 
attempt  was  made  to  collect  data  in  different  areas  in  order  to  determine  the 
effect  of  such  factors  as  type  of  loading,  size,  processes;  surface  conditions, 
heat  treatment,  surface  environment,  temperature,  surface  treatment,  stress 
concentration,  etc. 


The  Mechanical  Properties  Data  Center  in  Traverse  City,  Michigan,  was 
found  to  be  a  very  useful  source  of  information  for  the  fatigue  strength  data. 
They  have  been  very  cooperative  in  providing  the  necessary  information.  An 
exhaustive  collection  of  data  on  the  fatigue  strength  of  non-ferrous  materials 
was  made  from  the  Mechanical  Properties  Data  Center.  These  data  were  than 
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systematized,  evaluated  in  terms  of  th(f  distribution  parameters  (Section  6), 
tabulated  (Appendix  1) ,  and  the  representative  ones  plotted  (Section  6) . 

| 

4.4  FACTORS  AFFECTING  THE  STATISTICAL  DISTRIBUTION  OF  FATIGUE  STRENGTH 

I 

Since  fatigue  strength  represents  tthe  major  interest  in  the  engineering 
applications  of  the  Interference  Theory,  this  problem  was  studied  in  some  de¬ 
tail.  The  statistical  distribution  of  the  fatigue  strength  of  a  mechanical 
component  is  a  function  of  a  number  of  factors,  such  as  type  of  loading,  sur¬ 
face  finish.,  stress  concentration,  heat  treatment,  temperature,  environment, 
size  of  specimen,  frequency  of  loading,  surface  treatment,  manufacturing  pro¬ 
cesses,  grain  size,  impurities,  and  time.  Each  shows  variability  which  is 
characterized  by  some  form  of  a  distribution.  The  effects  of  these  factors  on 
the  statistical  distribution  of  strength  were  studied  in  the  present  investi¬ 
gation.  v 

Fatigue  strength  can  be  defined  as  the  maximum  stress  that  can  be  sus¬ 
tained  for  a  specified  number  of  cycles  without  failure,  the  stress  being  com¬ 
pletely  reversed  within  each  cycle.  In  the  case  of  steels,  a  component  is 
said  to  have  finite  fatigue  strength  if  it  fails  between  10^  and  10&  or  10? 
cycles  due  to  a  given  magnitude  of  cyclic  load.  For  non-ferrous  materials  this 
is  generally  extended  to  5  x  10$  cycles.  In  the  present  investigation  the 
data  were  extrapolated  to  10&  cycles. 

Type  of  Loading:  The  three  major  types  of  fluctuating  load  encountered 
in  designing  parts  are  axial,  bending,  and  torsion.  Experimentally  determined 
values -of  the  ratio  of  average  fatigue  strength  for  axial  loading  as  compared 
to  bending  load  were  reported  in  literature  as  ranging  generally  from  0.75  to 
1.0. (4»5)  Although  a  great  deal  of  work  has  been  done  to  obtain  precise  values 
for  this  ratio,  no  detailed  study  has  ever  been  made  as  to  the  statistical  as¬ 
pects  of  these  strengths.  Investigations  have  been  conducted  to  find  statis¬ 
tical  distributions  (Normal,  Exponential,  Weibull,  etc.)  of  fatigue  strength 
tested  under  a  given  type  of  loading,  such  as  bending.  No  work  was  done  to 
determine  the  effect  on  the  distribution  if  the  loads  were  other  than  bending. 
In  the  present  investigation,  an  attempt  was  made  to  study  the  effect  of  dif¬ 
ferent  loads  on  the  statistical  distribution  of  the  fatigue  strength.  The 
statistical  parameters  of  the  distribution  for  various  materials  under  differ¬ 
ent  loads  were  determined,  tabulated  according  to  materials  (Appendix  1),  and 
plotted  in  Section  6. 

Surface  Finish:  The  surface  finish  of  a  part  does  affect  its  endur¬ 
ance  strength.  Hence,  the  condition  of  finish  should  be  taken  into  account 
when  the  design  is  based  on  fatigue.  Surfaces  which  have  an  effect  on  the 
significant  strength  can  be  classified  into  five  broad  categories:  polished, 
ground,  machined,  hot -rolled,  and  as-forged.  The  worse  the  surface  condi¬ 
tion  the  lower  will  be  the  mean  fatigue  strength  but  the  higher  will  be  the 
scitttr.  A a  a  result,  the  degree  .of  interference  is  likely  to  be  pronoun¬ 
cedly  affected  by  the  type  of  surface  finish  Imparted  to  the  member.  Dif¬ 
ferent  surface  affects  were  studied  in  the  present  investigation  and  the 
fatigue  strength  parameters  were  tabulated  (Appendix  1)  and  plotted  in 
Section  6. 


Stress  Concentration:  A  notch  or  a  stress  raiser  in  a  part  subjected 
to  fatigue  loading  can  be  regarded  as  a  factor  causing  a  local  increase  in 
stress  or  a  reduction  in  strength.  For  example,  a  notch  with  a  stress  concen¬ 
tration  factor  of  2  can  be  thought  of  as  doubling  the  jtress  or  as  halving  the 
strength.  In  the  present  investigation  this  factor  was  taken  as  a  strength 
reduction  factor. 

If  all  parts  were  made  of  materials  which  are  completely  homogeneous 
and  have  perfectly  polished  surface  finishes,  the  effect  of  a  notch  would  be 
to  increase  the  stress  by  the  factor  Kt.  Since  actual  materials  are  not  per¬ 
fectly  homogeneous  and  actual  surfaces  are  seldom  perfectly  polished,  there 
exist  internal  and  surface  stress  raisers.  For  this  reason,  the  addition  of  a 

notch  to  a  part,  already  having  stress  concentration  due  to  geometry,  generally 
produces  a  smaller  effect  than  would  be  predicted  from  the  theoretical  stress 
concentration  factor,  Kt.  The  extent  to  which  a  notch  reduces  the  endurance 
limit  of  a  part  is  referred  to  as  the  fatigue  stress  concentration  factor,  or 
the  fatigue  strength  reduction  factor,  and  is  designated  by  the  symbol  Kf. 

This  is  defined  as : ( 6) 

K  m  endurance  limit  of  specimen  without  the  notch 
f  endurance  limit  of  specimen  with  the  notch 

In  this  study,  an  attempt  vi<s  made  to  determine  the  effect  of  stress 
concentration  on  the  statistical  scatter  of  the  fatigue  strength.  More  speci¬ 
fically,  the  objective  was  to  find  out  whether  this  factor  changes  the  mean 
strength  only,  whether  it  has  an  effect  or.  the  standard  deviation,  or  whether 
it  completely  changes  the  nature  of  the  distribution  Itself.  The  data  were 
collected  for  various  materials,  various  stress  concentrations,  and  at  differ¬ 
ent  temperatures.  Changes  in  the  parameters  of  the  statistical  distribution 
of  the  strength  due  to  the  effect  of  stress  concentration  for  different  mate¬ 
rials  at  various  testing  temperatures  are  tabulated  in  Appendix  1  and  plotted 
in  Section  6. 

Heat  Treatment:  Different  heat  treatments  such  as  annealing,  solution 
heat  treating,  normalizing,  quenching,  tempering,  aging,  etc.,  can  be  imparted 
to  materials  to  change  their  mechanical  properties.  Heat  treatment  may  change 
the  average  fatigue  strength  but  also  the  statistical  scatter.  Pertinent 
parameters  are  tabulated  in  Appendix  1  and  plotted  in  Section  6. 

Temperature:  In  a  similar  manner  the  effect  of  temperature  on  the  sta¬ 
tistical  scatter  of  fatigue  strength  is  shown  in  Appendix  1  and  Section  6. 

Environment:  The  effect  of  environment  on  the  fatigue  strength  was 
confined  in  this  investigation  to  the  corrosive  environment  since  most  of  the 
empirical  strength  data  was  In  this  area.  It  is  known  that  a  corrosive  envi¬ 
ronment  has  a  significant  effect  on  the  fatigue  strength,  particularly  at  a 
higher  life.  The  statistical  distribution  parameters  for  various  non-ferrous 
materials  under  different  corrosive  environments  were  determined,  tabulated 
according  to  materials  (Appendix  1) ,  and  plotted  (Section  6) . 
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Other  Factors:  In  the  present  study  it  was  observed  that  size  and  fre¬ 
quency  of  loading  do  not  have  any  consistent  effect  on  the  distribution  of  the 
fatigue  strength,  whereas  surface  treatment,  manufacturing  processes,  and  im¬ 
purities  do  exhibit  consistent  effects.  The  reason  is  that  the  latter  may 
change  the  microscopic  structure  which,  in  turn,  is  related  to  strength.  The 
statistical  distribution  parameters  for  various  non-ferrous  materials  under 
the  above  conditions  were  determined,  tabulated  (Appendix  1)  and  plotted 
(Section  6) . 

4.5  COMPUTER  ANALYSIS  OF  STRENGTH  DATA 

Data  collected  during  this  phase  of  the  investigation  were  organized 
and  systematized  according  to  materials  and  conditions.  A  computer  approach 
was  developed  to  analyze  these  data  (Section  6) .  In  the  past  investigation 
using  the  graphical  approach,  only  the  Weibull  distribution  was  used  in  the 

data  analysis  for  the  reasons  discussed  in  Reference  1.  In  the  present  inves¬ 
tigation,  with  the  aid  of  the  computer,  other  distributions  besides  Weibull, 
such  as  Largest  Extreme  Value,  Smallest  Extreme  Value,  and  Normal,  were  tried. 
The  reasons  for  the  selection  of  these  particular  distributions  are  given  in 
Section  6.3.  The  best  fitting  distribution  and  its  parameters  were  determined 
for  various  materials  and  conditions  and  these  are  tabulated  according  to 
materials  in  Appendix  1. 


SECTION  5  ANALYTICAL  EXPRESSIONS  FOR  INTERFERENCE 

5.1  INTRODUCTION 

5.1.1  Interference  Probabilities 

In  interference  theory  one  supposes  that  the  strength  of  a  manufactured 
part  is  not  known  with  certainty  prior  to  performing  some  test  on  it  and  that 
the  stress  induced  by  a  load  is  not  known  with  certainty  prior  to  actually 
loading  the  part.  Thus,  for  example,  one  does  not  know  with  certainty  that 
the  strength  of  a  part  is  exactly  50  ksi.  He  may  know  that  the  part  cannot 
have  a  strength  greater  than  58  ksi  or  less  than  40  ksi.  Or  he  may  know  that 
the  average  strength  that  has  been  obtained  in  previous  tests  on  these  parts 
is  49  ksi.  He  may  have  some  measure  of  how  dispersed  the  strength  measures 
are  around  this  average  strength.  The  point,  of  course,  is  that  this  type  of 
knowledge  is  quite  different  from  knowing  precisely  what  the  strength  is  prior 
to  testing.  For  a  multiplicity  of  reasons,  strengths  of  seemingly  identical 
parts  are  not  exactly  the  same, and  precisely  what  strength  a  part  will  have 
cannot  be  known  until  some  type  of  strength  test  is  performed.  In  the  theory 
of  probability  one  says  that  the  strength  of  a  part  is  a  random  variable.  Cer¬ 
tainly  the  same  type  of  reasoning  applies  to  the  stress.  Thus  for  a  mathema¬ 
tical  theory  of  interference  one  starts  with  the  idea  that  strength  is  a  ran¬ 
dom  variable,  say  X,  and  stress  Is  a  random  variable,  Y. 

In  describing  the  properties  of  random  variables,  since  their  values 
are  not  known  exactly,  one  supposes  that,  associated  with  every  set  of  values 
that  the  random  variable  can  take,  there  is  a  real  number  called  the  probabil¬ 
ity  that  the  random  variable  takes  values  in  the  set.  These  probabilities 
are  non-negative  real  numbers,,  they  are  all  less  than  1  and  in  the  sense  given 
below  they  "sum"  to  1. 

If  x  is  any  real  number  then  there  is  a  probability  that  the  random 
variable  takes  some  value  less  than  or  equal  to  x.  Symbolically, 

Pr(X  *  x) 

is  a  number  such  that  0  *  Pr(X  *  x)  5  1.  Surely  (i.e.  with  probability  1) 

X  -  «  so  that  * 

Pr  (X  -  »)  -  1  , 
and 

Pr (-»  >  X)  ■  0* 

Clearly  Pr (X  *  x)  depends  on  the  real  number  x.  Consequently  one 
defines  a  probability  distribution  function.  P(x),  by  the  relation 


One  sees  Immediately  that 

0  *  F(x)  *  1, 

F(— )  “  0 
F(«)  -  1 

and  that  F(x)  Is  a  non-decreasing  function. 


In  most  engineering  applications  F(x)  has  a  derivative  for  every  value 
of  x,  and  one  defines  the  probability  density  function,  f (x) ,  by 


f(x) 


ills! 

dx 


One  takes  f(x)dx  ■  Pr(x  <  X  *  x  +  dx). 

(i.e.  the  probability  density  function  multiplied  by  dx  is  the  probability 
that  X  takes  values  ir.  the  neighborhood  of  x.)  Since 

f (x)dx  -  dF(x)  , 

one  has 

00 

£ 

the  probabilities  given  by  f(x)dx  “sum"  to  1. 

In  the  mathematical  theory  of  interference  one  assumes  that  the  proba¬ 
bility  density  function  for  the  random  variables  X (strength)  and  Y  (stress) 
are  known.  Sections  6  and  7  following  show  how  these  functions  can  be  found 
from  engineering  data.  Thus  the  "givens"  of  the  mathematical  theory  of  inter¬ 
ference  use  the  random  variables  X  and  Y,  the  set  of  values  that  they  each  can 
take  (usually  the  non-negative  real  line),  and  the  probability  density  or  dis¬ 
tribution  function  F(x)  (or  f (x) ]  and  G(y)  [or  g(y)]. 

The  problem  to  which  Interference  theory  addresses  itself  is  that  of 
finding  the  probability  of  failure.  Failure  is  said  to  occur  whenever  the 
stress  exceeds  strength.  Thus  from  the  known  probabilities  for  the  X  and  Y 
random  variables  one  wishes  to  find 

Pr(Y  *  X)  , 

which  is  the  probability  that  strass  excaeds  strength  or  the  probability  of 
failure. 

5.1.2  Calculation  of  Probabilities  of  Failure 

There  are  two  useful  ways  of  determining  the  probability  of  failure 
from  the  known  properties  of  X  and  Y. 


f(x)dx 


I 


dF(x)  -  F(x) 
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(a)  Since  one  wishes  to  find  Pr(Y  *  X)  it  is  convenient  in  some  cases 
(e.g.  when  stress  and  strength  are  normally  distributed)  to  define  a  new  ran¬ 
dom  valuable,  Z,  by  the  relation 

Z  -  X  -  Y  . 

Then  if  one  can  find  the  probability  density  function  of  Z,  h(z),  the  proba¬ 
bility  of  failure  will  be  simply  the  probability  that  z  s  0,  In  terms  of  h(z) 
this  is  found  by 

Pr (failure)  - 


The  problem  in  general  is  then  to  find  h(z)  from  the  known  probability  density 
functions  f(x),  g(y).  Ai  complete  discussion  of  this  method  and  its  applica¬ 
tions  is  found  in  Sectioii  A-3.2  of  Reference  [l]. 

In  the  important  special  case  in  which  both  stress  and  strength  are 
normally  distributed  random  variables  it  is  well  known  that  Z  is  also  normally 
distributed  with  parameters. 

yZ  "  PX  -  PY 

2  2  2 
and  °Z  "  °X  +  aY 

Consequently  the  probability  of  failure  can  be  found  directly  from  tables  of 
the  normal  curve  areas.  One  wants  the  area  from  -*  to  0  from  these  tables. 

A  complete  discussion  of  how  to  do  this  is  given  in  Section  2. 

(b)  For  most  applications  method  (a)  above  is  unnecessarily  complex 
because  one  must  first  find  the  entire  density  function  of  the  random  variable 
Z  before  finding  the  probability  of  failure.  Since  the  random  variable  1  is 
of  no  practical  value  for  Z  <  0,  the  approach  in  part  (a)  is  unduly  long.  The 
methods  described  in  this  section  are  more  direct  and,  from  our  experience, 
more  useful  in  general. 

One  can  derive  the  probability  of  failure  as  follows.  Suppose  we  sup4*» 
impose  the  stress  and  strength  density  function  on  the  same  graph  as  shown  in 
Figure  5.1. 

Although  Y  is  a  random  variable,  let  us  fix  attention  on  a  particular, 
small  interval  that  Y  can  take  values  in.  Let  us  fix  y  <  Y  *  y  +  dy.  Then 
let  us  find  the  probability  that  the  random  variable  X  takes  values  lass  than 
this  fixed  Y.  One  can  show  that  this  probability  is 

Pr(X  -  y  |  y  <  Y  s  y  +  dy)  f(x)dx  . 

The  left  hand  side  of  this  expression  is  called  a  conditional  probability.  It 
ia  the  probability  that  tha  random  variable  X  takas  values  less  than  the  real 
number  y  whan  it  is  known  ("given  that")  the  random  variable  Y  is  "nearly"  y. 


£ 


h(z)dz. 
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By  the  definition  of  f(x),  this  probability  is  obviously  the  same  as  the  right 
hand  side  of  the  expression.  If  now  we  multiply  this  conditional  probability 
by  -  \ 

*  \  Pr(y  <  Y  -  y  +  dy)  ■  g(y)dy  , 

we  obtain  the  joint  probability  that  X  *  y  and  y.  <  Y  -  y  +  dy  which  symboli¬ 
cally  is 

Pr(X  *  y  |  y  <  Y  5  y  +  dy)  Pr(y  <  Y  5  y  +  dy) 

»  Pr(X  -  y;  y  <  Y  *  y  +  dy)  .  , 

This  probability  is  given  by  the  integral 

f (x)dx  g(y)dy  . 


From  the  joint  probability  one  obtains  the  probability  of  failure  by: 


Pr (failure) 


■H 


f(x)  g(y)  dx  dy  . 


o  o 


Thus  this  double  integral  gives  the  probability  of  failure  directly. 

If  one  recalls  the  relation  between  f(x)  and  F(x),  it  is  clear  that  the 
double  integral  is  easily  reduced  to  the  single  Integral 

Pr(failure)  *  F(y)  g(y)  dv  . 
vo 

If  F(x)  is  easily  obtainable  then  this  expression  is  easier  to  work  with  than 
the  double  integral.  For  example,  one  knows  that  for  a  strength  with  a  Weibull 
probability  density  function  the  probability  distribution  function  F(x)  is 

*lven  by  -}x-xq  m. 


F(x)  -  1  -  e 


In  these  cases  the  probability  of  failure  lo  then  given  by 

g-»0  I b* 

Pr(failure)  ■  /  (1-e  )  g(y)  dy 


-{ 


“*ol  b. 


g(y)  dy  . 


The  latter  expression  follows  because 


g(y)dy  •  1 


if  the  random  variable  Y  takes  only  positive  values  which  is  the  usual  case  in 
interference  theory. 

5.1.3  Interference  Tables.  Appendix  2 

In  Section  5.1.2  it  was  shown  that  the  probability  of  failure  could  be 
expressed  as  an  integral  involving  the  known  probability  density  or  distribu¬ 
tion  functions.  In  certain  cases  this  integral  can  be  evaluated  in  closed 
form  (e.g.  when  f(x)  and  g(y)  are  both  exponential  functions).  In  some  cases 
this  Integral  can  be  evaluated  in  terms  of  other  well  known  and  tabulated 
functions  (e.g.  when  f(x)  and  g(y)  are  both  gamma  functions  or  when  f(x)  and 
g(y)  are  both  normal  functions).  In  general  it  is  not  to  be  expected  that  the 
integral  for  the  probability  of  failure  can  be  evaluated  in  closed  form  or  in 
a  form  involving  other  well  known  functions,  (e.g.  when  f(x)  and  g(y)  are 
both  Weibull  functions  or  when  f(x)  is  a  Weibull  function  and  g(y)  is  a  normal 
function.)  In  those  cases  one  must  resort  to  numerical  evaluation  of  the 
integral. 

Since  the  integrals  giving  the  probability  of  failure  cannot  be  ex¬ 
pressed  in  terms  of  well  known  functions,  in  general,  we  have  evaluated  the 
integral  numerically.  Tables  of  the  probability  of  failure  are  given  in  Sec¬ 
tion  A-2  of  [l] .  A  full  discussion  of  the  numerical  methods  used  and  the 
errors  of  approximation  appropriate  to  the  tables  are  given  in  Section  A-4  of 
[1]. 


Additionally,  one  is  interested  in  the  case  in  which  the  strength  dis¬ 
tribution  is  of  the  form 

S(x-M) 

1  -  e  when  &  >  0, 

—to  —  X  -  +oo  , 


called  the  Smallest  Extreme  Value  probability  distribution  function.  One  is 
also  Interested  in  the  case  in  which  strength  is  distributed  as 

-B(x~M) 

e  e  when  B  >  0, 

-»  <  x  <  +*  » 


which  is  called  the  Largest  Extreme  Value  probability  distribution  function. 

Section  A-2. 2  gives  the  probability  of  failure  for  the  case  in  which 
strength  is  Weibull  distributed  end  stress  is  normally  distributed. 

Section  A-2. 5  gives  the  probability  of  failure  for  the  case  in  which 
strength  is  distributed  with  the  Smallest  Extreme  Value  distribution  function 
and  stress  is  normally  distributed. 
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Section  A-2.4  gives  the  probability  of  failure  for  the  case  in  which 
strength  is  distributed  with  the  Largest  Extreme  Value  distribution  function 
and  stress  is  normally  distributed. 

5.2  USE  OF  INTERFERENCE  TABLES  IN  APPENDIX  2 

5.2.1  Parameters  for  the  Weibull  Distributed  Strength.  Normal  Distributed 
Stress 

The  form  for  the  integral  involved  in  finding  the  probability  of  failure 
when  the  strength  is  Weibull  distributed  and  the  stress  is  normally  distribu¬ 
ted  is  given  in  Section  A-3.3.5  of  [lj.  Tables  of  these  probabilities  are 
given  in  Section  A-2.1  of  [ l] .  A  discussion  of  the  numerical  analysis,  error 
and  accuracy  of  the  tables  is  given  in  Section  A-4.1.7  of  [lj. 

The  form  of  the  distribution  of  the  strength  has  been  given  in  Section 
5.2.1  of  [1]. 

In  the  past  investigation,  extensive  tables  of  Interference  values 
(Section  A-2.1. 2  of  [ 1 ] )  were  prepared  for  the  parameter  C  (as  defined  below) 
equal  to  10  or  above.  In  the  present  investigation  on  non-ferrous  materials, 
the  values  of  parameter  C  much  lower  than  10  were  needed.  Hence,  the  Tables 
A-2.1. 2  of  [l]  were  extended  to  cover  the  range  of  1  s  C  <  10.  Therefore,  the 

t  list  probabilities  of  failure  for 

-  1,  1.2,  1.3,  ...  3.2 

■  1.,  2 . 10. 

■  0,  .2,  .4,  ...  2.8  and  -.2  to  -10.0 

•  the  slope  of  the  strength  distribution.  In  the 
tables  this  is  called  B(x). 

-  the  ratio  of  the  difference  of  the  characteristic 
strength  and  the  truncation  parameter  (or  lower 
bound  of  strength)  to  the  standard  deviation  of  the 
stress.  In  the  tables  this  is  called  C  for  typo¬ 
graphical  simplicity. 

-  the  difference  between  the  strength  truncation  pan¬ 
meter  and  the  mean  stress  divided  by  the  standard 
deviation  of  the  stress.  In  the  tables  this  is 
called  A  for  simplicity. 

body  of  the  table  are  the  probabilities  of  failure 
at  the  heading  of  the  tables.  From  the  discussion 
of  [l]  these  probabilities  are  correct  to  *5  x  10“^. 


Tables 


where 


of  this  repor 

B(x) 

C 

A 


(0  -x  >/o 

X  o 


(xo-u)/o 


The  values  in  the 
for  the  parameters  given 
given  in  Section  A-4.1.7 


5.2.2  Parameters  for  the  Smallest  Extreme  Value  Distributed  Strength. 

*Nonnal  Distributed  Stress 

In  this  case  it  has  been  found  that  the  integral  giving  the  probability 
of  failure  can  be  expressed  in  terms  of  two  parameters  which  are  denoted  by 
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a  and  y.  In  the  tables  these  parameters  are  given  by 

y  ■  $(y-M)  ■  the  difference  in  the  location  parameter  (y)  for 
the  normally  distributed  stress  and  the  mode  (M)  for  the 
Smallest  Extreme  Value  distributed  strength.  This  differ¬ 
ence  is  multiplied  by  6,  the  slope  parameter  of  the  extreme 
value  distribution. 

o  ■  So  ■  the  product  of  the  standard  deviation  of  the  Normal 
distribution  for  the  stress  variable  and  the  slope  para¬ 
meter  for  the  Smallest  Extreme  Value  distributed  strength. 

The  tables  give  the  probability  of  failure  for 

a:  .001,  .005,  .025,  .05,  .075,  and  .1  to  3.0  in 
steps  of  0.1  and  3.5  to  10  in  steps  of  0.5. 

y:  0  to  -A0  in  steps  of  1.0 

The  values  in  the  tables  give  the  probability  of  failure  corresponding 
to  the  values  of  a,y  given  at  the  column  and  row  heads.  These  probabilities 
are  correct  to  within  *2  x  10“4  according  to  the  discussion  given  in  Section 
A-3. 

5.2.3  Parameters  for  the  Largest  Extreme  Value  Distributed  Strength. 

Normal  Distributed  Stress 

In  this  case,  the  integral  expressing  the  probability  of  failure  is 
almost  identical  to  that  in  the  previous  case  (Section  5.2.2).  The  defini¬ 
tions  of  the  parameters  a  and  y  are  identical.  The  tables  give  probabilities 
of  failure  for  the  ranges  of  a  and  y  listed  in  Section  5.2.2.  These  proba¬ 
bilities  are  also  correct  to  *2  x  10-4  as  discussed  In  Section  A-3. 

5.2.4  Use  of  the  Tables.  Explanation  of  Missing  Values  and  Interpolation 

Numerical  examples  of  the  use  of  the  tables  are  given  in  Section  9.  In 
general  the  user  will  enter  the  table  with  known  parameters,  for  example,  b  , 
©x,  Xq,  and  the  appropriate  parameters  for  the  stress  distribution,  and  wlshx 
to  find  the  probability  of  failure.  This  is  a  direct  table  look-up.  In  some 
design  problems  the  user  will  have  a  given  probability  of  failure  to  achieve 
and  will  know  the  general  shape  of  the  distribution  of  stress  and  strength 
appropriate  to  the  material  that  he  is  using.  The  table  will  then  give  him 
the  relative  parameters  (there  may  be  many  of  these)  to  design  for.  It  would 
be  expected  that  a  coat  analysis  would  give  the  acceptable  parameter  values 
for  each  distribution.  As  long  as  the  relative  values  are  as  given  in  the 
table,  the  probability  of  failure  will  be  the  same  no  matter  what  the  values 
for  each  distribution  are. 

In  the  Welbull  Strength-Normal  Stress  tables,  the  missing  values  within 
the  ranges  of  the  parameters  used  are  nearly  zero  and  hence  have  not  been  tabu¬ 
lated.  When  A  ■  (x0-y)/o  >3.5,  the  probability  of  failure  is  less  than  1  x 
IQ-*  since  the  area  under  the  normal  curve  from  3.5  to  *  is  less  than  3  x  10~4. 
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For  the  Normal-Extreme  Value  cases,  the  probabilities  of  failure  are 
zero  (to  4-places)  whenever  Y  <-39  and  a  <10. 

When  o  «  go  =  o  (or  o  =  0) ,  then  the  stress  distribution  for  either 
the  Normal-Smallest  Extreme  Value  or  Normal -Largest  Extreme  Value  case  is  a 
straight-line  distribution,  and  the  probabilities  of  failure  can  be  computed 
as  in  Appendix  2.  These  formulas  are  repeated  here  since  they  give  very  good 
approximations  for  the  probabilities  even  when  a  ^  0  but  is  small  (say  a<.00l). 


For  a  =  0: 

Pr  (failure^}  =  1  -  e~e 

(Normal-Smallest) 

„Y  where Y  =  8  (M-M) 

Pr[failur0  »  e"e 

(Normal -Larges t) 

It  can  be  seen  that  the  tables  are  non-linear  for  almost  all  values 
of  the  parameters.  This  can  cause  inaccuracies  when  the  tables  are  inter¬ 
polated.  The  absolute  value  of  the  interpolation  error  depends  on  which 
tables  are  interpolated.  For  precise  values  the  user  should  use  a  higher 
order  interpolation  formula  (as  given  in  Section  A-2.1)  rather  than  linear 
interpolation.  We  have  not  explored  the  relative  errors  of  interpolation 
closely.  In  those  cases  checked,  the  relative  errors  are  small. 
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SECTION  6  STATISTICAL  DISTRIBUTION  OF  FATIGUE  STRENGTH 
BY  COMPUTER  APPROACH 

Most  fatigue  testing  involves  subjecting  specimens  or  parts  to  a  fluctuating 
stress  to  failure  and  repeating  this  process  at  various  stress  levels.  The 
data  thus  obtained,  known  as  life  data,  are  used  to  construct  the  conventional 
S-N  diagram.  In  this  case,  the  scatter  obtained  is  the  scatter  in  life  at  a 
given  stress.  In  the  present  investigation  the  attention  was  focused  on  the 
nature  of  the  scatter  in  fatigue  strength  at  a  given  life.  In  the  past  inves- 
tigatlon(l),  where  the  study  was  made  mainly  on  the  ferrous  materials,  the 
problem  of  determining  the  scatter  in  the  fatigue  strength  and,  subsequently, 
its  distribution  function  was  resolved  by  a  graphical  method  discussed  in 
Reference  1.  Among  the  several  possible  distribution  functions  that  were  con¬ 
sidered,  only  the  Weibull  distribution  was  used  because  of  its  wide  use  and 
the  difficulty  of  handling  other  distributions  by  a  graphical  method.  In  the 
present  investigation  on  the  non-ferrous  materials  a  computer  approach  was 
developed  by  which  the  determination  of  the  distribution  functions  of  strength 
became  considerably  less  difficult.  Consequently,  it  was  possible  to  try  out 
such  other  distributions  as  Largest  Extreme  Value,  Smallest  Extreme  Value, 
Logistic,  and  Normal  besides  the  Weibull  which  was  tried  in  the  past  investi¬ 
gation.  The  reasons  for  selection  of  these  particular  distributions  are  given 
in  Section  6.3.  The  computer  program  was  so  developed  that  when  the  raw  data 
(conventional  S-N  type  data)  were  fed  into  this  program,  the  computer  printed 
the  degree  of  fxt  and  the  parameters  for  all  the  above  distributions  for  a 
given  set  of  strength  data.  In  the  first  part  of  this  program  the  scatter  in 
the  life  was  converted  into  the  scatter  in  the  fatigue  strength.  In  the  se¬ 
cond  part,  the  distribution  functions  (Weibull,  Normal,  Logistic,  Largest  Ex¬ 
treme  Value,  and  Smallest  Extreme  Value)  were  fitted  to  these  strengths  data 
simultaneously,  and  the  one  with  the  highest  degree  of  fit  was  the  best  fitting 
distribution.  The  detailed  discussion  on  this  computer  approach  and  the  ne¬ 
cessary  statistical  tools  that  were  used  in  developing  the  approach  are  given 
in  the  following  sections.  The  flow  chart  of  the  computer  program  u3ed  to 
determine  the  fatigue  strength  distribution  is  given  in  Aopendlx  4. 


6.1  LEAST  SQUARES  METHOD  OF  FITTING  A  LINE  TO  S-N  TYPE  DATA 

6.1.1  Method  of  Least  Squares 

The  principle  underlying  the  "best  fit"  of  a  line  by  the  least  squares 
method  is  that  the  sum  of  the  squares  of  the  deviations  in  the  y-direction  of 
the  data  points  (Zi,  Yj,)  from  the  most  probable  line  is  minimum  (Sea  Figure 
6.1).  The  equation  of  this  line  is: 

Y  ■  o  +  $X  (6.1) 

where  a  is  intercept  on  Y-axis 

and  6  is  the  slope  of  the  line. 

The  parameters  o  end  0  are  found  from  the  following  relationships (2) : 


76 
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Figure  6.1  Least  Squares  Line 
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(6.2) 


(6.3) 


where  n  -  the  total  number  of  points 
■  abscissas  of  points 
Yi  -  ordinates  of  points 


6.1.2  Application  of  Least  Squares  Method  to  S-N  Data 

When  stress-life  (S-N)  data  are  plotted  on  a  log-log  paper  or  log  (S) 
and  log  (N)  are  plotted  on  a  Cartesian  paper,  an  approximately  straight  line 
results.  Hence,  the  S-N  data  are  related  by  the  equation: 

In  (S)  -  a  +  0  In  (N) 

When  In  (S)  and  In  (N)  are  substituted  for  Y  and  X 
and  (6.3),  they  reduce  to: 

n|  I  (In  Nt)  (In  S^j  -|J  In  Njj  J  J  In  S 

■(£  <id  Ni>1  -  LI 10  Hi  1 

n  n 

l  (In  Sj)  -  el  In  N. 

..at - ±=± — : 


where  a  end  8  ere  the  Intercept  end  the  slope  of 
Figure  6.2) 

n  -  the  total  number  of  test  points 
%  -  the  ith  life,  corresponding  to  the  ith  stress, 


in  the  Equations (6.2) 

J 

(6.4) 

(6.5) 

the  S-N  line  (See 
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6.1.3  Correlation  Coefficient 

Correlation  Coefficient  (R)  is  an  index  of  the  goodness  of  fit  of  the 
least  squares  line  to  a  set  of  data.  In  the  previous  section  it  was  pointed 
out  that  a  straight  line  can  be  fitted  to  a  set  of  data  with  known  coordinate^ 
for  example,  In  (S)  and  In  (N) .  This  section  deals  with  a  computation  of  cor¬ 
relation  coefficient  (R)  which  is  a  measure  of  the  goodness  of  fit.  The  cor¬ 
relation  coefficient  is  defined  as: 

R  -  (S)  -f  (6.6) 

S 

y 


where  0  ■  the  slope  of  the  fitted  line. 
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(6.7) 


(6.8) 


If  the  correlation  coefficient  (R)  la  1.0,  every  data  point  falls  on 
the  least  square  line  and  the  data  has  the  beat  fit.  On  the  other  hand,  if 
the  correlation  coefficient  (R)  is  equal  to  zero,  all  the  data  are  around  a 
circle  and  the  data  has  very  poor  fit.  (Sec  Plgure  6.3) 


\ 
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Five  distributions  (Weibull,  Normal,  Logistics,  Largest  Extreme  Value, 
and  Smallest  Extreme  Value)  were  fitted  to  each  set  of  the  fatigue  strength 
data,  and  the  values  of  the  correlation  coefficient  (R)  for  each  fitted  distri¬ 
bution  were  computed  from  Equation  (6.6).  The  distribution  having  the  largest 
value  of  R  was  taken  as  the  best  fitting  distribution. 

6.2  CONVERSION  OF  LIFE  DATA  TO  STRENGTH  DATA 

In  the  past  investigation^  on  ferrous  materials  the  conversion  of  life 
data  to  strength  data  was  accomplished  graphically  where  the  data  were  plotted 
on  the  conventional  S-N  diagram  and  the  least  squares  S-N  curve  was  then  fit¬ 
ted  to  the  test  points.  Passing  through  each  point  an  S-N  curve  parallel  to 
the  least  squares  curve  was  drawn.  These  made  a  set  of  parallel  S-N  lines. 

(See  Figure  6.4).  A  vertical  line  was  drawn  at  a  given  life,  say  N^,  inter¬ 
secting  the  family  of  S-ji  curves.  The  points  of  intersection  Sj_,  S2,  S3  ... 
represent  the  scatter  of  the  fatigue  strength  at  the  life  Ni  (See  Figure  6.4). 
Here,  it  was  assumed  that  to  each  specimen  of  the  population  can  be  attributed 
an  individual  S-N  curve,  and  that  there  exists  for  any  population  of  specimens 
(at  fixed  test  conditions)  a  family  of  non-intersecting  S-N  curves  which  can 
be  determined  with  any  desired  accuracy,  each  curve  corresponding  to  a  given 
probability. 

In  the  present  investigation  on  non-ferrous  materials,  this  basic  as¬ 
sumption  was  held  still  true,  and  for  the  purpose  of  developing  the  computer 
approach,  it  was  necessary  to  develop  an  analytical  approach  for  the  above 
graphical  method.  This  is  given  in  the  following  sections. 

6.2.1  Computation  of  Parallel  S-N  Lines 

The  general  form  of  the  equation  of  the  parallel  S-N  curves  is: 

In  S  -  0^  +  8  In  N  (6.9) 

where  6  is  the  slope  of  the  lines  as  computed  from  Equation  (6.4) 

ou  is  the  Intercept  of  the  line  passing  through  a  given  point 
(Ni,  St) 

and  1  is  the  subscript  which  takes  values  from  1  to  n  where  n  is 
the  number  of  parallel  lines 

In  order  to  determine  the  equation  of  each  line,  it  was  necessary  to 
compute  the  corresponding  values  of  from  the  following  equation: 


a±  -  (In  St)  -  6  (In  Hj)  .  (6.10) 

By  substituting  the  values  of  S-N  data,  (S^,  N^,  (S2*  N2),  (S3,  N3)...,  in 
Equation  (6.10),  a^,  <g,  03...  can  bs  computed.  Hence,  using  Equation  (6.9), 
a  set  of  equations  for  each  parallel  line  was  obtained. 


N 


1 


Life,  cycles  (log  scale) 


Figure  6. A  S-N  Diagram  for  converting  Life  Data 
to  Strength  Data 
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(6.11) 


InS  *  +  3  InN 
InS  ■  02  +  6  InN 
InS  *03+6  InN 

6.2.2  Computation  of  the  Scatter  of  Strength  at  a  Given  Life 

A  life,  say  N^,  was  chosen  at  which  the  scatter  of  strength  was  desired 
By  substituting  In  N^  in  Equation  (6.11),  a  set  of  values  of  In  Si,  In  S2, 

In  S3  ...,  were  obtained.  By  taking  antilog  of  these  values,  the  scatter  of 
fatigue  strengths  Si,  S2,  S3  ...  for  the  life  Ni  was  obtained.  This  was  re¬ 
peated  for  several  values  of  life  Nj ,  and  the  resulting  scatters  of  fatigue 
strengths,  S^,  were  stored  by  computer  for  further  analysis. 

Median  ranks (2)  were  then  assigned  to  the  scatter  of  the  fatigue 
strength  data,  Sij ,  arranged  in  increasing  order.  The  analytical  expression 
for  the  median  ranks  that  was  used  for  the  computer  programming  was (7): 

Median  Rank  *•  ^  r~v 
n  +  .4 


where  j  »  order  number  and 

n  -  total  number  of  dat?  p-'ints. 

6.3  DISTRIBUTION  FUNCTIONS 

Hie  most  common  probability  distributions  that  were  found  from  litera¬ 
ture.  for  expressing  the  fatigue  strengths  are:  Weibull,  Largest  Extreme  Value, 
Smallest  Extreme  Value,  Normal,  and  Logistic.  These  cover  a  considerable 
variety  of  distribution  patterns.  For  example,  Weibull  distribution  generates 
a  wide  variety  of  distributions  by  choosing  different  values  of  its  parameters, 
Xo,  9,  and  b,  (this  is  discussed  in  the  following  section).  Extreme  Value 
distribution  (Largest  and  Smallest)  takes  care  of  some  of  the  unexplainable 
phenomenon  observed  in  connection  with  the  fracturing  of  materials  under  the 
applied  dynamic  load  where  the  puzzling  question  1st  Why  is  the  strength  of  a 
specimen  considerably  smaller  or  largsv  than  its  expected  value?  The  reason 

for  this  difference  in  the  values  lies  in  the  fact  that  there  may  exist  flaws 
in  the  specimen  that  will  weaken  it.  Normal  distribution  vss  used  for  the 
obvious  reason  that  it  esn  be  justly  called  the  most  common  and  well  known 
distribution.  Logistic  distribution,  which  has  somewhat  the  seme  character¬ 
istics  as  the  Normal  distribution,  was  tried,  but  later  it  was  found  from 
the  computer  results  that  Normal  fits  the  data  just  as  well  as  Logistic. 

Hence,  in  those  cases  where  Logistic  fitted  the  data  best,  the  Normal  was 
used,  Instead,  to  express  the  fatigue  strength.  Besides,  Normal  distribution 
is  widely  used  and  is  easier  to  understand  for  engineers  with  little  statis¬ 
tical  background.  It  la  for  these  reasons  that  the  discussion  on  the  Log¬ 
istic  distribution  does  not  appear  in  the  following  section.  Other  distrib¬ 
utions  are  discussed  below. 
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6.3.1  Weibull  Distribution 

The  Weibull  equation  is  a  three  parameter  function.  The  general  expres¬ 
sion  for  the  Weibull  density  function  is: 


x-Xob-i 


(6.12) 


X  <  x  ^  “ 
o 


and  the  general  expression  for  the  cumulative  distribution  function  is: 


l*-X ol  b 


P(x)  “  1  -  e 


where,  as  used  in  this  study, 


,  X 
’  o 


(6.13) 


x  is  the  independent  variable  (fatigue  strength) 

Xq  is  the  lower  bound  of  fatigue  strength 

0  is  the  characteristic  fatigue  strength,  where  63.2%  of  the 
population  have  strengths  less  than  or  equal  to  this  value (2) 

b  is  the  Weibull  slope 
6.3.2  Normal  Distribution 

Normal  distribution  is  a  two  parameter  function.  The  general  expression 
for  Normal  density  function  is: 


1  2o 

e 


,  -«  <  x  <  +“ 


(6.14) 


where  x  is  the  independent  variable  (fatigue  strength) 
U  is  the  population  mean 

a  is  the  population  standard  deviation 
By  e  simple  transformation,  this  can  be  reduced  to: 


p(y)  -7^2  6 


(6.15) 
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where  z  *  standardized  normal  variate 
and  P(y)  -  probability  of  z  *  y. 

6.3.3  Extreme  Value  Distributions 


2=11 

0 


Gumble(8)  designed  what  is  known  as  the  Extreme  Value  Distribution 
which  is  generally  used  to  estimate  the  tail-ends  of  the  frequency  distribu¬ 
tions.  This  distribution  determines  the  probability  of  occurrence  of  an  event. 
Largest  values  or  smallest  values  of  an  event  can  be  estimated  by  Largest  Ex¬ 
treme  Value  distribution  or  Smallest  Extreme  Value  distribution  respectively. 
The  general  expression  for  the  density  function  is: 


p(x)  -Be 


Largest  Extreme  Value  (L.E.V.): 
-  B  (x-M) 


3  e 


-  3 (x-M)  -  e 

-  3  (x-M) 


-e 


-3 (x-M) 


,  -—  —  x  *  +* 


Smallest  Extreme  Value  (S.E.V.): 


p(x)  -  3  e 


+3  (x-M)  -  e 


+B  (x-M) 


B  e+e  (x"M)  •  e'e 


+3  (x-M) 


*  x  *  +» 


(6.16) 

(6.17) 


(6.18) 


The  general  expression  for  the  cumulative  distribution  function  is: 

Largest  Extreme  Value  (L.E.V.): 

-6  (x-M) 

P  -  e”e  (6.19) 


Smallest  Extreme  Value  (S.E.V.): 


+fl  (x-M) 

P  -  1  -  e”6 


(6.20) 


where  B  -  intensity  function  (slope) 


M  -  Extreme  Value  (mode) 


3  and  M  are  also  called  the  extremal  parameters.  Comparison  of  characteris¬ 
tics  of  Extreme  Value  distribution  with  Normal  distribution  is  given  in  Table 
6.1  which  is  self-explanatory. 
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Characteristic 


Largest  Value  Smallest  Value  Normal  Varial 


A.  Distribution 


1. 

Cumulative  probability  func- 

-y 

V 

z  J. t2 

tion  for  reduced  variate,  4> 

-e 

e 

1-e 

i  737  *  d‘ 

2. 

Reduced  variate  - 

y*  a(x-u) 

y=  a(x-u) 

Z-i(x-u) 

3. 

Original  variate,  x,  in 
terms  of  reduced  variate  - 

x=u+y/ a 

x-u+y/ a 

x=u4<JZ 

B.  Reduced  Variate  y 


4. 

Moment  generating  function 
G(t) - 

r(i-t) 

r(i+t) 

e^2 

5. 

Mode  J  - 

0 

0 

z=0 

6. 

Mean  y  - 

Y-0.57722 

-0.57722 

z”0 

7. 

Median  y  - 

-lglg2-036651 

-0.36651 

z»0 

8. 

Standard  deviation  o  — - - 

W//ST-1.28255 

1.28255 

1 

9. 

Skewness  0^  - 

1.29857 

-1.29857 

0 

10. 

Kurtosis  0£ - - 

27/5 

27/5 

3 

C.  Deviates  y$  Corresponding  to  Given  Probability  Points  $ 


Table  6.1  Comparison  of  Characteristics  of  Extreme  Value  Distribution 
and  Normal  Distribution^) 


Characteristic  Largest  Value  Smallest  Value  Normal  Variate 


D.  Area  Pkg  Included  Within  Mode  *  Ka 


II  _ 1///////, 

-K  r 

'A//////// 1 _ 

0  +K*- 

y 

13.  K-l;  PKo-*(o)-*(-o) - 

0.73064 

0.73064 

0.68269 

14.  K-2;  PKa=<K2a)-$(-2a) - 

0.92597 

0.92597 

0.95450 

15.  K-3;  PKa=<t(3a)-4>(-3a) - 

0.97890 

0.97890 

0.99730 

16.  Variate  y  for  which  Py»0. 68269 

1.14071 

1.14071 

1.00000 

17.  Variate  for  which  Py-0^95450— 

3.06685 

3.06685 

2.00000 

E.  Original  Variate, 

X 

18.  Mode  X - 

u 

-u 

U 

19.  Mean  X - 

u-Hy  /  o 

-u-y/a 

y+aZ*y 

20.  Median  x - 

u+0. 36651/ a 

-u-0. 36651/ a 

y+0Z*y 

21.  Standard  deviation  ax - 

1  „  1  ir 

1  1  IT 

a  °  a# 

a-1 

Table  6.1  (Continued)  Comparison  of  Characteristics  of  Extreme 
Value  Distribution  and  Normal  Distribution^®' 
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6.4  THE  LINEARIZING  TRANSFORMATION  EQUATIONS  OF  THE  DISTRIBUTION  FUNCTIONS 


When  x,  the  independent  variable  (in  this  case  the  fatigue  strength), 
and  F,  the • dependent  variable  (in  this  case  the  cumulative  probability  of  the 
fatigue  strength),  are  plotted  on  a  Cartesian  graph  paper,  the  resultant  plot 
is  non-linear  (sigmoidal  or  S-shaped  plot)  as  shown  in  Figure  6.5. 

Such  a  curve  is  difficult  to  work  with.  To  fit  this  curve  to  a  set  of 
data  it  is  necessary  to  use  the  curvilinear  regression  analysis.  This  analy¬ 
sis  is  quite  involved  as  compared  to  the  linear  regression  analysis  which  can 
be  used  to  fit  only  the  straight  lines.  Hence,  if  the  non-linear  function  can 
somehow  be  transformed  into  a  linear  function,  the  linear  regression  analysis 
can  be  employed.  This  can  be  done  by  transforming  the  dependent  (P)  and  inde¬ 
pendent  (x)  variables  of  the  function  to  a  new  set  of  dependent  [T(P) ]  and 
independent  [G(x)]  variables.  This  process  is  called  the  Linearizing  Trans¬ 
formation.  The  functional  relationships  between  the  variables  (independent 
and  dependent)  of  the  sigmoidal  form  of  the  distribution  (Figure  6.5)  and 
variables  of  the  linear  (Figure  6.6)  form  of  the  distribution  are  called  the 
linearizing  transformation  equations.  These  are  given  as: 

X  ■  G(x)  -  independent  variable  plotted  on  abscissa 

Y  ■  T(P)  -  dependent  variable  plotted  on  ordinate 

The  linearizing  transformation  is  used  in  many  engineering  applications.  For 
example,  in  strength-life  (S-N)  diagrams  with  the  ordinate  and  the  abscissa 
having  logarithm  scale; 

Y  -  In  S 

and  X  ■  In  N  are  the  linearizing  transformation  equations  where 
In  S  and  In  N  are  called  the  linear  mode  variables  and  S  and.N  are  the  origi¬ 
nal  variables.  Or,  in  case  of  Weibull  distribution,  )>nln  j_p(x)  and  In  (x-X0) 
are  the  linear  mode  variables  where  as  P(x)  and  x  are  the  original  (or  non¬ 
linear  mode)  variables.  In  the  analysis  of  fatigue  strength  data,  lineariza¬ 
tion  transformation  can  be  stated  as: 

Y  -  T(P)  (6.21) 

X  -  G(x)  (6.22) 

where  T  and  G  are  the  linearizing  transformation  equations  relating  original 
variables  P  and  x  to  linear  mode  variable  Y  and  X  respectively.  The  functions 
G  and  T  should  be  such  'that  the  inverse  transform  operation  on  them  gives  the 
original  variables .  These  inverse  transformation  equations  are: 

P  -  T”1  (Y) 

x  -  G"1  (X) . 

The  inverse  transformation  equations  are  important  for  finding  the  parameters 
of  original  distribution  functions.  This  will  be  discussed  in  the  section 
that  follows. 


1.00  r  I .  "I .  I  1  0.975 


cm  + 

CM  <H 


m 

o 

o 

o 

o 

o 

o 

o 

o 

o 

in 

m 

<N 

ON 

ON 

oo 

vO 

m 

CO 

CN 

•H 

o 

O 

• 

• 

* 

• 

• 

• 

• 

• 

• 

«. 

. 

• 

o 

o 

o 

o 

o 

o 

o 

© 

o 

o 

o 

o 

ax®aS  .iBauTTuoN  -  paxT®i  ^uaoaaj 


aT®°S  «a«H  -  (d)  l  «  i 


cm  f 

CM  •* 


F? 

e 

iH  4J 
CM  CO 


o  o  o  o  o 
9X«3S  Jvaupi  -  paXT«i  ^uaoaaj 


39 


Figure  6.5  Sigmoidal  Mode  of  Cumula-  Figure  6.6  Straight  Line  Mode  of  Cumu- 

tive  Distribution  Function  lative  Distribution  Function 


Either  or  both  the  abscissa  scale  and  the  ordinate  scale  may  be  trans¬ 
formed  for  linearizing  the  frequency  distribution  functions.  In  the  following 
section  the  linearizing  transformation  equations  are  derived. 

6.4.1  Weibull  Distribution 


The  general  expression  for  the  Weibull  cumulative  distribution  function. 


■X-X0. b 


Equation  6.13,  is 

p  -  1  -  e 
with  original  variables 

P  ■  probability  of  failure 
x  «  strength 
or, 

1 

- - -  e 

1  -  p 


XQ  *  x  *  00 


lnln  (~^)  -  b  In  (x-XQ)  -  b  In  (0-XQ)  . 

(6.23) 

i  a  form  of  Y  -  Bx  +  a  where: 

Y  -  lnln 

1-P 

(6.24) 

X  -  In  (x-XQ) 

(6.25) 

6  *  b 

(6.26) 

a  -  -b  In  (0-XQ) 

(6.27) 

The  Equations  (6.24)  and  (6.25)  respectively  for  ordinate  and  abscissa  scales, 
are  the  linearizing  transformation  equations  for  Weibull  distribution.  Hence, 
if  the  new  scales  (ordinate  and  abscissa)  are  constructed  corresponding  to 
these  equations,  the  points  following  Weibull  distribution  should  plot  as  a 
straight  line  on  these  new  scales. 

6.4.2  Normal  Distribution 

In  the  case  of  Weibull  distribution,  the  density  function  (Equation 
6.12)  can  be  integrated  in  the  closed  form  to  give  the  cumulative  distribution 
function  (Equation  6.13), 


P(x)  -  1  -  e 


x-X  b 

<<wT> 


O 


i 
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whereas  in  the  case  of  Normal  distribution,  it  is  not  possible  to  integrate 
the  density  function  in  the  closed  form  to  give  Normal  cumulative  distribution 
function.  Hence,  this  is  expressed  as: 


P(x) 


dx 


Whereas,  in  the  case  of  Weibull  distribution,  the  Equation  6.13  can  be  trans¬ 
formed  into  a  straight  line  function  (See  Equation  6.23),  there  is  no  simple 
transformation  for  the  case  of  Normal  distribution  that  can  be  used  to^relate 
the  original  variable  P(x)  to  a  linear  mode  variable, which  is  lnln  l-pY^y  in 
the  case  of  Weibull.  Hence,  the  following  method  was  used  to  linearize  the 
Normal  distribution. 

The  general  equation  of  jhe  linear  form  (straight  line)  of  a  distribu¬ 
tion  is : 

Y  «  a  +  8X  (6.28) 

where  X  is  some  function  of  strength  x  [for  example,  X  ■  In  (x-XtJ)1for  Weibull] 
and  Y  is  some  function  of  probability  P(for  example,  Y  *  lnln  T-p)*  Since 
Y  is  not  a  simple  function  of  probability  (P),  one  needs  to  find  a  new  variable 
which  should  satisfy  the  Equation  (6.28),  and  it  should  be  related  to  P  In 
such  a  manner  that  it  can  be  used  to  find  necessary  values  of  this  new  variable 
for  any  given  corresponding  value  of  P. 

Such  a  variable  can  be  the  standardized  normal  variate,  z.  The  stan¬ 
dardized  normal  variate,  z; 

1.  is  related  to  P  values  and  could  be  found  directly  from  the  Tables 
of  Normal  distribution,  such  as  given  in  Section  2,  Table'  2.1. 

2.  satisfies  the  Equation  (6.28).  This  is  shown  below: 
z  is  expressed  as: 

z  .S* 
a 

or  z  -  (k  X  -  (J) 

0  0 

This  has  the  same  form  as  Equation  (6.28),  where 

Y  •  z,  —  ■  $,  ■  a,  and  x  -  strength. 

oo 

Hence,  the  linearizing  transformation  equations  for  Normal  distribution  are: 

Y  -  z,  the  Standardized  Normal  Variate  (6.29) 

X  -  x  (6.30) 
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If  the  transformation  just  described  does  linearize  the  Normal  distri¬ 
bution  function,  it  will  plot  as  a  straight  line  on  a  Cartesian  paper.  A 
graphical  example  demonstrating  the  linearization  of  Normal  distribution  is 
shown: 


The  Normal  distribution  function  with  a  mean  of  y  -  28  and  a  «  4.25  was 
plotted  on  a  Normal  probability  paper  (See  Figure  6.7)  from  which  the  follow¬ 
ing  values  of  x  and  its  corresponding  values  of  P  were  selected: 


x 

15 

20 

25 

35 

For  the  corresponding  values  of  P, 
Tables  of  Normal  distribution: (9) 


X 

P 

15 

.0015 

20 

.035 

25 

.025 

35 

.95 

.0015 

.035 

.25 

.95 

the  values  of  z  were  found  from  the 
T{P)  -  z  -  2=11 

g 

2.97 

1.82 

0.67 

-1.65 


The  values  of  x  on  the  abscissa  and  z  on  the  ordinate  were  plotted  on  a  Carte¬ 
sian  coordinate  paper,  and  the  resultant  plot  was  found  to  be  a  straight  line 
(See  Figure  6.8).  Hence,  this  proved  that  this  method  did  linearize  the  Nor¬ 
mal  distribution  and  the  linearizing  transformation  equations  are  Equations 
(6.29)  and  (6.30); 


Y  -  z 

and  X  -  x  . 


6.4.3  Extreme  Value  Distributions 


The  general  expression  for  the  cumulative  distribution  function  of  the 
Largest  Extreme  Value  distribution  is,  from  Equation  6.20): 


P  -  e 


-e 


-8  (x-M) 


By  taking  the  In  twice,  this  equation  reduces  to  a  form  of  a  straight  line: 


or 


-lnln  (±)  -  fix  -  &M 
Y  -  gx  +  a 


(6.31) 
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Percent  Cumulative  Probabilit 
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Plgurs  6.7  Plot  of  x  vs.  P  on  Normal  Probability  Papar 
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where  Y  ■  -lnln  (^) 

(6.32) 

X  -  x 

(6.33) 

a  -  -6M 

(6.34) 

Equations  (6.32)  and  (6.33)  are  linearizing  transformation  equations 
for  Largest  Extreme  Value  distributions.  For  Smallest  Extreme  Value  dl stribu' 
tions,  the  linearizing  transformation  Equations  (6.35)  and  (6.36)  were  found 
in  a  similar  manner.  They  are: 

Y  -  lnln  1  l  p 

(6.35) 

X  -  x 

(6.36) 

a  -  -6M 

(6.37) 

6.5  FITTING  THE  LINEARIZED  DISTRIBUTION  TO  THE  STRENGTH  DATA 


In  the  previous  sections,  the  parameters,  the  linearizing  transforma¬ 
tion  equations,  and  the  linear  mode  variables  of  the  distribution  functions 
used  in  this  study  were  presented.  The  best  fitting  distribution  is  the  one 
which  has  the  highest  correlation  coefficient,  R.  In  order  to  determine  the 
values  of  R  for  each  distribution,  they  should  be  fitted  to  Probability- 
Strength  (P,x)  data  by  fitting  a  least  squares  line  (Y  »  a  +  BX)  to  the  linear 
mode  variables  of  these  data.  Ibis  was  done  for  all  the  distributions  in  the 
following  manner: 


6.5.1  Weibull  Distribution 


x). 


The  original  variables  of  the  data  are  Probability  and 
From  the  linearizing  transformation  Equations  (6.24)  and 


Strength  ()P  and 
(6.25); 


Y  - 


and  X  •  In  (x-Xc)  ,  the  linear  mode  variables 

and  In  (x-X0),  where  x  is  the  original  variable  (strength),  P 
probability  from  the  median  rank  table,  Xq  is  the  lower  bound 
X  and  Y  are  the  linear  mode  variables. 


are  lnln  yT'p~ 

is  the  cumulative 
of  strength,  and 


Zy  substitution  of  X  and  Y  in  Equations  (6.2)  and  (6.3),  the  values  of 
the  parameters  o  and  6  for  the  least  squares  line  (Y  -  a  +  BX)  were  computed: 


6 


(lnln 


1-P 


)  in(x 


i-V  - 

.1  ..  v . 


n  r  ri 

n  l  |ln(xi-X  )}2  -  l  InH  -X  ) 

1-1  U-i  *  J 


(6.38) 
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and  a 


n 


(6.39) 


(lnln 


8  Jl  ln  <xl-X0) 


In  order  to  compute  a  and  8  it  was  first  necessary  to  determine  X0  from 
the  Probability-Strength  (P,  x)  data.  Once  the  X  was  determined,  the  values 
of  a  and  8  were  computed,  and  the  correlation  coe?ficient  R  was  found  by  sub¬ 
stituting  the  values  of  a  and  8  in  Equation  6.6.  In  the  past  investigation^), 
X0  was  determined  by  a  graphical  method;  in  the  present  investigation  this  was 
done  by  a  computer  approach,  as  follows: 

1.  A  value  for  XQ  was  assumed  which  should  be  somewhere  between  zero 
and  the  lowest  strength  value  of  the  test  data.  Let  this  X_  be 

VI). 

2.  Once  the  value  of  Xq  is  known  or  assumed,  the  original  variables 
of  the  data  points  can  be  transformed  to  the  linear  mode  variables 
by  using  the  transformation  equations 

Y  -  lnln  ~ 

X  -  ln  (x-X0). 

3.  By  substitution  of  each  value  of  X  and  Y  corresponding  to  each 
value  of  the  original  data  points  (P,x)  in  Equations  (6.2)  and  (6.3), 
the  parameters  a  and  8  of  the  least  squares  line  Y  -  a  +  3X  were 
then  computed. 

4.  Using  the  values  of  a  and  8,  the  value  of  correlation  coefficient  R 
for  this  fit  was  computed  from  Equation  (6.6).  Let  this  correla¬ 
tion  coefficient  be  R(l). 

5.  Another  Xo  was  then  assumed  close  to  (greater  or  less  than)  Xq(1). 
Let  this  be  Xo<2).  The  steps  l  through  4  were  repeated.  Let  the 
correlation  coefficient  thus  computed  be  R(2). 

6.  The  next  value  of  Xq-X^(3)  was  then  computed  in  the  following  man¬ 
ner:  (See  Figure  (6.9) 

Whenever  the  ratio  *8  positive,  the  value  of 

XqO)  was  assumed  to  be  larger  than  Xq(2);  and  if  the  ratio  is  ne¬ 
gative,  the  value  of  Xq(3)  was  assumed  to  be  smaller  than  Xo(2). 

The  larger  the  standard  deviation  of  the  strength  data,  the  larger 
will  be  the  difference  [Xq(3)  -  Xq(2)]  that  can  be  assumed. 

7.  After  Xo(3)  was  found,  the  steps  2,3,  and  4  were  repeated;  the  cor¬ 
responding  value  of  correlation  coefficient  R(3)  was  computed. 

8.  This  was  repeated  twelve  times,  and  twelve  Xq' s  and  their  corres¬ 
ponding  values  of  R  were  computed.  (Note:  The  number  of  repeti¬ 
tions,  twelve,  was  selected  erbltrsrlly) . 
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Correlation  Coefficient 


Out  of  these  12  lower  bound  X0’s,  the  one  which  gave  the  maximum  correlation 
coefficient  (R)  (the  best  fit)  was  taken  as  the  "true"  lower  bound  of  strength 
X0.  Once  the  true  value  of  Xq  was  determined,  this  and  the  values  of  x  and  P 
were  substituted  in  Equations  (6.38)  and  (6.39)  to  compute  a  and  6;  ahd,  sub¬ 
sequently,  the  correlation  coefficient  R  was  determined  from  Equation  (6.6). 

6.5.2  Normal  Distribution 


To  fit  a  Normal  distribution  to  a  set  of  data  means  fitting  a  straight 
line  (least  squares  line  Y  =>  a  +  8X)  such  that: 

Y  *  z,  the  standardized  Normal  variate 
and  X  =  x,  the  fatigue  strength. 

The  values  of  z  correspond  to  the  values  of  median  ranks,  and  they  were  found 
from  the  Normal  Distribution  Table. (9) 

For  example,  if  the  total  number  of  test  points  is  three,  their  core¬ 
sponding  values  of  median  ranks  are  (2):  .20.63,  .5000  and  .7937  as  read  from 

the  median  ranks  table.  These  values  represent  the  areas  under  the  Normal 
distribution  curve  (See  Figure  6JL0).  Hence,  the  values  of  z  corresponding  to 
the  values  of  the  median  ranks,  .2063,  .5000,  and  .7937  are  -.8194,  0,  and 
+.8194  as  read  from  the  Normal  Distribution  Table.  The  values  of  z  or  median 
rank  z  values  were  found  in  the  similar  manner  for  various  values  of  median 
ranks.  The  table  of  median  rank  z  values  similar  to  the  table  of  median  ranks 
was  constructed,  see  Table  6.1.  These  tables  of  median  ranks  and  median  rank 
z  values,  which  are  in  matrix  form,  were  stored  into  the  computer  for  the  pur¬ 
pose  of  the  data  analysis.  The  values  of  X  •  x  (fatigue  strength)  and  Y  -  z 
(the  corresponding  value  of  the  standardized  Normal  variate)  were  substituted 
in  Equations  6.2  and  6.3,  and  a  and  8  were  computed  to  give  the  equation  of 
the  least  square  line,  Y  ■  o  +  BX,  fitted  through  these  points. 

The  next  step  was  then  to  determine  the  correlation  coefficient  R. 

This  was  done  using  the  Equation  6.6. 

Since  the  Normal  distribution  is  a  little  more  involved  than  the  other 
distributions  (Weibull  and  Extreme  Valued  an  example  is  given  to  illustrate 
the  procedure  of  fitting  the  Normal  distribution  to  data. 

Suppose  the  fatigue  strength  distribution  data  at  10s  life  cycles  are: 
27.0,  24.0,  32.5,  and  29.0  ksi.  (Although  the  number  of  test  points  were  at 
least  10  or  more  in  our  study,  four  test  points  are  given  to  make  the  correla¬ 
tion  simple) 

The  data  wara  then  arranged  in  an  Increasing  order  of  value  and  the 
appropriate  median  ranks  z  values  for  sample  size  n  ■  4  were  read  from  Table 
6.1  as  follows. 


Area  =  .  2063  ““5 

(Median  Rank)  \  S 


Area  ■  .7937 
(Median  Rank) 


-.8194  0  +.8194 


Figure  6.10  Illustration  for  the  Median  Rank  z  Values 


MEDIAN  RANKS  z  VALUES 


iber 


Sample  Size,  n 
3 


-0.5450 

-0.8194 

-0.9982 

-1.1290 

40.5450 

0.0000 

-0.2887 

-0.4826 

+0.8194 

+0.2887 

0.0000 

+0.9982 

+0.48(26 

+1.1290 

Table  62 

Median 

Ranks  z  Values 

x.  ksi 


Median  Ranks  z  Values 


24.0 

27.0 

29.0 

32.5 


-0.$982 

-0.2887 

+0.2887 

+0.9982 


By  substitution  of  X  ■  x  and  Y  -  z  values  in  Equations  (6.2)  and  (6.3), 
they  reduced  to:  n  '  n  n 

n  l  ziXi  -  l  xi  l  z± 

6  .  i-1  i-1  i-1 


n 


'n  1  xi  ~  I  I  xi 


i-1 


i-1 


and 


n  n 

I  H  "  .8  [ 

i-1  i-1 


The  values  of  c  and  8  were  then  computed  to  determine  the  equation  of  the 
least  squares  line,  Y  -  o  +  BX.  Knowing  the  value  of  8,  the  correlation  coef¬ 
ficient  R  was  determined  from  Equation  6.6: 


R  -  (8) 


2Sx 


(8) 


2  Sx 


where 


n  l  \  ’  (  T  XJ 

i-1  1  li-1  » 


n(n-l) 


and 


l  \2  -  ( r  ■  J 

i-i  1  \i-i  i 


n(n-l) 


Sample  Variance 
of  x 

(abscissa) 


Sample  Variance 
of  z 

(ordinate) 


6.5.3  Extreme  Value  Distributions 


The  linearising  transformation  equations  for  Largest  Extrema  Value  dis¬ 
tribution  (Equations  6.35  and  6.36)  are: 


and 


T  -  -lnln  (±) 
X  -  x  . 


By  substitution  of  X  and  T  in  Equations  (6.2)  and  (6.3),  the  values  of 
the  parameters  a  end  6  for  the  least  square  line  (T  -  o  +  BX)  were  computed: 
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6.6  DETERMINATION  OP  WE  BEST  FITTING  DISTRIBUTION  AND  ITS  PARAMETERS 

In  the  previous  sections,  the  ststisticsl  distributions,  their  para~ 
aeters,  end  their  linearising  transformation  equations  were  discussed.  The 
distributions  acre  then  fitted  to  the  data  by  fitting  a  least  squares  line 
(Y  »  a  ♦  8X)  to  the  linear  node  variablea  of  these  data  by  coaputing  the 
values  of  a  and  0  from  Equations  (6.2)  and  (6.3).  Knowing  o  and  0,  the 
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the  correlation  coefficients  R  were  determined  from  Equation  (6.6).  The  best 
fit  distribution  was  the  one  which  had  the  largest  value  of  correlation  coef¬ 
ficient  R.  This  distribution  could  be  taken  as  the  "true"  representation  of 
the  scatter  of  the  strength  data.  The  linear  mode  of  the  distributions  [for 
example,  ln(x-Xc)  and  (lnln  £?£)  for  Weibull]  are  not  used  in  practice,  and 
therefore,  the  original  parameters _(b,  Xq,  and  9, for  Weibull)  of  the  distri¬ 
bution  should  be  computed.  These  parameters  were  determined  from  the  slope  6 
and  the  interception  a  of  the  fitted  least  squares  line.  This  could  be  done 
by  using  the  inverse-transformation  equation  as  mentioned  in  Section  6.4.  The 
specific  inverse-transformation  equations  to  compute  the  parameters  of  the 
four  distributions  used  in  this  study  are  given  in  the  following  sections. 

6.6.1  Weibull  Distribution 


As  mentioned  before,  Weibull  distribution  has  three  parameters.  They 

are: 

XQ  *  lower  bound  of  strength 
b  ■  the  slope 

9  -  characteristic  strength 

The  method  of  determining  the  value  of  Xq  was  explained  in  Section  6.5.1. 
b  is  also  the  slope  of  the  least  squares  line  whose  parameters  are  a  and  3. 
Therefore,  b  ■  3.  9  was  determined  from  Equation  6.27: 

u  -  -bln  (9  -  Xq) 

i»  <•-*„>-  -  b 

a 

Therefore,  9  ■  e  11  +  Xq 

and  b  ■  3 

6.6.2  Normal  Distribution 

The  Normal  distribution  has  only  two  parameters.  They  are: 


(6.44) 

(6.45) 


o  ■  standard  deviation 

The  linear  form  of  the  distribution  is  expressed  as 

«  •  J  , 

00 

and  the  equation  of  the  least  squares  line  is  i 

Y  -  8X  +  o  . 

Hence,  the  parameters  of  this  line  are: 


a 


3 


1 

a 


.  _  a  . 

a  ’ 

a  and  y  were  computed  from  rearranging  these  equations : 

1 

°  “  3 

a 

u  -  -aa  -  -  j 


(6.46) 

(6.47) 


6.6.3  Extreme  Value  Distributions 

Both  of  these  distributions  have  two  parameters,  and  they  are  referred 
to  as  extremal  parameters.  However,  in  this  study,  these  parameters  are  re¬ 
ferred  to  as: 


b  ■  slope 

M  ■  Mode 

The  slope  b  is  the  slope  of  the  least  squares  line  (as  it  was  the  case 
for  Weibull  distribution).  Hence,  b  ■  0. 

The  mode  M  was  determined  from  Equation  6.34: 

M  -  -  f  (6.48) 

and  b  ■  8  (6.49) 

These  equations  would  apply  to  determine  the  parameters  (M,3)  for  both 
Largest  and  Smallest  Extreme  Value  distributions  even  though  the  numerical 
values  of  a  and  B  would  be  different  in  each  case. 

Using  the  statistical  and  analytical  tools  discussed  here,  s  computer 
program  was  developed  whereby  the  computer  prints  out  the  fatigue  strength  dis¬ 
tribution  functions,  their  parameters,  and  the  corresponding  degree  of  fit 
(correlation  coefficient  R)  for  several  preassigned  lives. 

Data  on  various  non-ferrous  materials  undar  various  conditions  were 
collected  and  analysed  by  computer-aided  procedure  as  described  above.  The 
results  are  tabulated  in  Appendix  1  in  terms  of  the  best  fitting  distribu¬ 
tion  function  and  its  parameters  for  several  lives.  The  most  representative 
parameters  were  then  plotted,  as  shown  in  Figures  6.14  to  6.64.  The  discussion 
as  related  to  these  plots  is  given  in  Section  6.7. 

6.7  GRAPHS  OP  THE  FATIGUE  STRENGTH  DISTRIBUTION  PARAMETERS 

6.7.1  Expression  of  Data  in  Terms  of  Normal  Parameters 

By  aid  of  the  computer  approach  discussed  in  Section  6,  the  fatigue 
strength  data  were  analysed  where  different  distributions  (Weibull,  Normal, 
Largest  Extreme  Value,  and  Smallest  Extreme  Value)  were  fitted  to  a  set  of 
data.  The  computer  printed  out  the  parameters  and  the  degree  of  fit  (correla¬ 
tion  coefficient  R)  for  each  of  these  distributions  for  a  given  set  of  data. 
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This  was  done  for  various  materials  under  various  conditions.  A  representa¬ 
tive  computer  output  sheet  for  a  material  and  conditions  shown  below  is  given 
in  Table  6.2. 


7075  Aluminum  (Bar),  8^=86  ksi 

10"*  cycles 

Rotary  Beam  Bending 

Kt  -  2.0 

See  Page  233 
and  Code  Number  (240) 

The  distribution  with  maximum  value  of  correlation  coefficient  R  was 
considered  to  be  the  best  fitting  distribution.  For  example,  for  the  data 
given  in  Table  6.2,  the  maximum  value  of  R  is  0.948724,  and  this  corresponds 
to  Largest  Extreme  Value  distribution.  Hence,  Largest  Extreme  Value  distribu¬ 
tion  is  the  best  fitting  distribution  of  the  fatigue  strength  of  the  above 
material  under  the  given  conditions.  Its  parameters  are: 

BETA  (8)  »  .301634 
and  MODE  (M)  *>  35.155 

The  fatigue  strength  data  of  various  non-ferrous  materials  under  various  con¬ 
ditions  (effects)  were  analyzed  in  a  similar  manner.  The  results  were  tabu¬ 
lated  In  terms  of  only  the  best  fitting  distribution  (this  could  be  Largest 
Extreme  Value,  Smallest  Extreme  Value,  Weibull,  or  Normal)  and  its  parameters 
according  to  the  materials  in  Appendix  1 .  The  tabulated  values  were  rounded 
to  four  digits.  It  was  felt  that  the  representative  parameters  should  be  plot¬ 
ted  to  study  the  effects  of  such  factors  as  stress  concentration,  surface  fin¬ 
ish,  type  of  loading,  etc.  on  the  distribution  of  the  fatigue  strength.  This 
may  lead  to  confusion  as  illustrated  by  the  following  example.  Say,  for  a 
given  material  and  conditions,  the  Weibull  parameters  are  XQ,  0,  and  b  with 
Kt  =  2.4.  However,  for  the  same  material  and  conditions,  if  the  design  should 
be  changed  so  that  Kt  becomes  1.5,  the  distribution  function  may  change  from 
Weibull  to  Smallest  Extreme  Value  (with  parameters  6  and  M) ;  and  if  Kt  should 
become  1.0,  it  may  change  to  Normal  (with  parameters  p  and  a). 

When  this  information  is  plotted  to  determine  the  effect  of  stress  con¬ 
centration  on  the  distribution  parameters,  a  plot  such  as  in  Figure  6.11  re¬ 
sults.  This  obviously  cannot  make  much  sense  to  the  reader,  as  it  is  very 
difficult  to  compare  cne  distribution's  parameters  (say  XQ„  b,  and  9)  for  one 
value  of  effect  (say  Kt  =  2.4)  wich  another  distribution’s  parameters  (say  8 
and  M)  for  another  value  of  the  effect  (say  Kt  =  1.5).  This  problem  can  be  re¬ 
solved  by  plotting  all  the  data  in  terms  of  the  same  distribution.  This  means 
if,  for  example.  Normal  distribution  is  chosen,  then  the  parameters  of  the 
other  distributions  should  be  expressed  in  terms  of  Normal  parameters.  This 
can  be  done  by  using,  for  all  the  test  data,  the  values  of  the  Normal  para¬ 
meters  directly  from  the  computer  sheet  (a  sample  is  shown  in  Table  6.2).  Al- 


Material : 

Life: 

Type  of  Loading: 

Theoretical  Stress  Concentration  Factor: 
Other  Conditions : 
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****************************************************************************** 


WEIBULL 


SMALLEST  EXTREME  VALUE 


THETA  (6) 

Xd 

b 

R 


-  37.8210 

-  30.5332 

-  1.9047 

-  0.889976 


BETA  (8) 
MODE  (M) 

R 


.301694 

38.701 

0.863512 


****************************************************************************** 


LARGEST  EXTREME  VALUE 


NORMAL 


BETA  (6)  -  .301634  MEAN  (y)  -  36.9284 

MODE  (M)  -  35.155  STD.  DEVIATION  (o)  -  4.2207 


R  -  0.948724  R  -  0.882308 

****************************************************************************** 


Table  6.2  A  Computer  Output  Data  Illustrating  the  Distribution  Parameters  and 
the  Corresponding  Correlation  Coefficients  for  a  Given  Set  of  Data 


5S 


Effect  of  Stress  Concentration 


Actual  Distribution  Parameters  for  Various  Stress  Concentration  Factors 


though  in  some  cases  the  best  fitting  distribution  might  be  other  than  the 
Normal  distribution,  it  would  not  loose  the  significance  by  more  than  5  to  .10%. 

For  example,  from  Table  6.2,  if  Normal  distribution  with  R  ■  0.882308  (that  is, 

88.2%)  is  used  instead  of  Largest  Extreme  Value  (the  best  fitting  distribution 
for  this  particular  case)  with  R  ■  94.9%,  the  level  of  significance  of  the  fit 
is  decreased  by  only  6.7%  (94.9  -  88.2  ■  6.7%).  Although  this  approach  results 
in  somewhat  lower  accuracy,  it  has  the  advantage  of  greater  clarity. 

Out  of  the  four  distributions  considered  in  the  present  investigation. 

Normal  distribution,  with  parameters  p  and  a,  was  selected  for  making  these 
plots.  The  reasons  for  the  selection  of  the  Normal  distribution  are: 

1.  Normal  distribution  with  its  parameters  y  (mean  or  average  value) 
and  o  (standard  deviation  or  measure  of  scatter)  is  easier  to  understand  for 
engineers  with  little  statistical  background. 

2.  The  plot  of  these  two  parameters  can  be  represented  by  a  single 
line  for  each  effect  (See  Figur  6.12)  as  compared  to,  say,  Weibull,  where 
three  lines  (one  for  each  o'  tui  three  parameters  XQ,  b,  and  9)  are  needed  to 
represent  each  effect. 

3.  Largest  or  Smallest  Extreme  Value  distribution  parameters  (0  and  M) 

are  not  as  easy  to  understand  as  in  the  case  of  Normal  parameters.  ' 

It  should  be  noted  that:  1)  the  plots  are  meant  only  to  represent  the  effect 
of  different  factors  on  the  fatigue  strength.  For  the  purpose  of  solving  the 
actual  design  problem,  it  is  recommended  that  the  actual  distribution  and  its 
parameters  be  used  from  the  tables  in  Appendix  1;  2)  only  the  representative 
values  from  the  entire  analysis  are  plotted,  whereas  the  tables  in  Appendix  1 
give  all  the  values. 

6.7.2  Explanation  of  the  Graphs 

Mean  strength  p  was  plotted  against  life  on  a  log-log  paper,  one  line 
representing  only  one  effect  (see  Figure  6,12).  To  the  right  of  the  graph  a 
specific  value  of  o/y  accompanies  each  line,  this  value  being  independent  of 
life.  Therefore,  once  y  is  determined  from  a  graph  for  a  given  life  and  given 
effect,  a  can  be  calculated  from  the  value. 

The  meaning  of  o/u  is  discussed  below: 

For  any  given  set  of  conditions  at  any  given  life  the  values  of  y  and  o  were 
generated  by  the  computer.  Consider  the  illustrative  example  in  Section  6.7.1 
graphically  shown  by  the  line  corresponding  to  Kt  ■  2.0  in  Figure  6.12, 

This  is  repeated  in  Figure  6.13.  At  10^  cycles,  mean  strength  is  y,. 

Locate  at  10*  cycles  (y,  +  a.), and  through  this  point  draw  a  line  parallel  to 
the  original  line.  (The  assumption  for  drawing  parallel  lines  is  justified  in 
Section  6.2,  specifically  in  Figure  6.4.)  At  10*  cycles  £n2Sti9  quoted 

above,  d-  4.2207  and  y*  36.9284,  and,  therefore,  o/y  »  ■  .1145.  This 

in  Figure  6.12  corresponds  to  °i/tfi  at  10*  cycles.  At  any°*  *  other  life  such 
as  10°,  locate  P2  and  (<*2  +  U2).  Tits  following  demonstrates  thst 
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02/1*2  "  al/wi  ■  .1145  and  is  thus  independent  of  life.  Since  the  lines  are 
parallel  (Figure  6  ..13), 


yr°i  »2-°2 

or  by  simple  algebra, 

°1  °2 

—  -  —  ■  Constant 


6.7.3  Use  of  Graphs 


In  this  manner  representative  data  were  plotted  for  various  materials 
and  under  various  conditions.  These  are  shown  in  Figure  6.14  to  Figure  6.64 
The  use  of  these  graphs  can  be  illustrated  by  a  representative  example  shown 
in  Figure  6.12.  For  a  life,  say,  105  cycles  and  Kt  -  2.0,  the  mean  strength 
y  -  28.0  ksi  and  o/y  ■  0.1145.  Therefore,  o  ■  0.1145  x  28  «  3.2  ksi. 
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INDEX  TO  GRAPHS  OF  NORMAL  PARAMETERS 
(Effect  of  Various  Factors  on  Mean  and  Standard  Deviation  of  Strength) 

Material  Normal  Parameters  can  be  found: 

Page  No,  Figure  No. 

ALUMINUM  ALLOYS 

7.5  Zn  -  2.5  Mg  Aluminum 


Effect  of  Environment 

67 

6.14 

Effect  of  Heat  Treatment 

67 

6.14 

M-257 

Effect  of  Stress  Concentration 

68 

6.15 

Effect  of  Test  Temperature 

68 

6.15 

M-276 

Effect  of  Stress  Concentration 

69 

6.16 

Effect  of  Test  Temperature 

69 

6.16 

2014  Aluminum 

Effect  of  Frequency  (T-4) 

70 

6.17 

Effect  of  Frequency  (T-6) 

70 

6.17 

Effect  of  Stress  Concentration 

71 

6.18 

Effect  of  Grain  Direction 

71 

6.18 

Effect  of  Manufacturing  Processes 

72 

6.19 

2014  Alclad  Aluminum 

Effect  of  Frequency 

73 

6.20 

2024  Aluminum  (Sheets) 

Effect  of  Specimen  Configuration 

74 

6.21 

Effect  of  Frequency 

74 

6.21 

Effect  of  Stress  Concentration 

75 

6.22 

Effect  of  Notch  Radius 

75 

6.22 

2024  Aluminum  (Bars) 

Effect  of  Stress  Concentration  (T-4) 

76 

6.23 

Effect  of  Stress  Concentration  (T-3) 

76 

6.23 

2024  Aluminum  (Extrusions) 

Effect  of  Manufacturing  Processes 

77 

6.24 
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Page  No.  Figure 

2025  Aluminum 


Effect  of  Heat  Treatment 

78 

6.25 

Effect  of  Surface  Finish 

78 

6.25 

Effect  of  Surface  Treatment 

79 

6.26 

Effect  of  Stress  Concentration 

79 

6.26 

2219  Aluminum 

Effect  of  Heat  Treatment 

80 

6.27 

Effect  of  Stress  Concentration 

80 

6.27 

2618  Aluminum 

Effect  of  Stress  Concentration 

81 
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SECTION  7 


THE  STATISTICAL  DISTRIBUTION  OF  STRESS 


7.1.  STRESS  SPECTRUM  VS  STRESS  DISTRIBUTION 

The  problem  of  stress  distribution  in  the  Interference  Theory  appears 
to  be  much  more  involved  than  the  problem  of  strength  distribution.  Consider, 
for  example,  the  problem  of  an  aluminum  cylinder  head  in  a  reciprocating 
engine.  Because  of  the  variation  in  hardness,  surface  finish,  etc.,  the 
fatigue  strength  will  vary  from  one  cylinder  head  to  another.  This  will  re¬ 
sult  in  a  distribution  curve  in  which  the  strength  will  be  plotted  on  the 
abscissa  and  the  number  of  cylinder  heads  having  a  given  strength  (i.e.,  fre¬ 
quency  of  occurrence)  on  the  ordinate. 

Consider  now  the  stress  distribution  in  the  cylinder  heads.  The 
stresses  in  the  head  result  from  the  gas  pressure  loading.  If  the  attention 
is  now  focussed  on  a  single  cylinder  head,  then  the  variation  in  the  gas  pres¬ 
sure  will  produce  a  distribution  of  stresses  in  this  particular  head.  The  re¬ 
sultant  curve  will  be  a  plot  of  the  stresses  in  the  head  on  the  abscissa  and 
the  number  of  times  that  this  stress  occurs  in  this  particular  head  on  the 
ordinate  [Figure  7-1  (a)]. 

This,  however,  is  not  what  is  wanted  in  the  application  of  the  Inter¬ 
ference  Theory,  because  this  distribution  of  stresses  cannot  be  matched  with 
the  distribution  of  strength.  In  the  strength  distribution  the  ordinate  gives 
the  number  of  cylinder  heads  having  a  given  strength.  Therefore,  in  the 
stress  distribution  the  ordinate  must  be  the  number  of  cylinder  heads  having 
a  given  stress  (and  not  the  number  of  times  a  given  stress  occurs  in  a  single 
head) .  This  can  be  obtained  by  considering  the  fact  that  different  engines 
will  be  subjected  in  service  to  different  operating  conditions,  and  therefore, 
the  distribution  of  gas  pressure  loading  will  vary  from  engine  to  engine.  As 
pointed  out  in  Section  7-2,  a  spectrum  of  stresses  must  be  converted  to  an 
equivalent  stress  for  the  purpose  of  Interference  Theory.  Therefore,  if  a 
spectrum  of  loading  due  to  different  service  conditions  varies  from  engine  to 
engine,  in  a  population  of  cylinder  heads,  the  equivalent  stress  will  vary 
from  head  to  head.  Thus  the  statistical  stress  distribution  desired  for  the 
Interference  Theory  may  be  obtained  [Figure  7.1  (b)].  In  this  distribution 
the  equivalent  stress  will  be  plotted  on  the  abscissa  and  the  number  of  cylin¬ 
der  heads  (frequency  of  occurrence)  having  that  stress  on  the  ordinate.  This 
distribution  then  can  be  compared  with  the  strength  distribution  to  obtain 
the  probability  of  Interference. 

7.2  CONVERSION  OF  A  STRESS  SPECTRUM  TO  AN  EQUIVALENT  STRESS  (8  ) 

By  definition,  the  equivalent  stress  is  a  completely  reversed  stress  of 
constant  amplitude  which,  when  imposed  on  a  part,  should  cause  failure  at  the 
same  life  as  if  the  stress  spectrum  was  imposed  instead.  Thus,  the  damage 
accumulated  at  any  given  life,  due  to  this  equivalent  stress,  will  be  the  same 
as  if  due  to  the  spectrum  of  stresses. 

The  first  step  towards  converting  the  spectrum  to  a  single  stress  (SeqU) 
is  to  convert  the  operating  stresses,  which  may  have  seme  mean  stress  associ¬ 
ated  with  them,  to  zero  mean  stress,  that  is,  to  a  completely  reversed  stress. 

v 
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Number  of  Times  that  the 

'umber  of  Cylinder  Heads  Cylinder  Head  is  Sub- 

Having  a  Given  Stress _ _ jected  to  a  Given  Stress 


(Figure  7.2).  This  can  be  done  by  means  of  the  modified  Goodman  diagram.  Draw 
the  Goodman  diagram  as  shown  in  Figure  7.3.  From  the  spectrum  of  operating 
stresses  plot  each  stress  cycle  on  this  diagram  as  shovm,  for  example,  line  AB. 
Connect  CA  and  CB  and  extend  sto  the  vertical  line  where  the  mean  stress  is 
equal  to  zero.  Hence,  XY  is  the  zero  mean  stress  equivalent  to  AB.  After  re¬ 
ducing  all  such  stress  cycles  to  zero  mean  stress,  the  stress  spectrum  will 
have  all  the  stress  cycles  completely  reversed.  The  magnitude  of  XY  will  be 
different  for  different  stress  cycles.  Therefore,  the  original  operating 
stress  spectrum  [Figure  7.2  (a)],  with  various  mean  stress  levels,  is  thus  re¬ 
duced  to  a  stress  spectrum  with  zero  mean  stress,  that  is,  to  a  completely 
reversed  stress  (Figure  7.4). 

This  spectrum  can  then  be  reduced  to  a  single  equivalent  stress  of  con¬ 
stant  amplitude  by  means  of  Miner's  or  Corten-Dolan's  Rule. 

7.2.1  Miner's  Rule 

Miner's  Rule^A^  assumes  that  the  total  life  of  a  component  can  be  esti¬ 
mated  by  simply  adding  the  fraction  of  life  consumed  by  each  over-stress  cycle. 
Overstress  can  be  defined  as  the  stress  above  the  endurance  limit  of  the  mate¬ 
rial  which,  if  applied,  will  damage  the  part.  This  rule  is  expressed  as: 


i-k  n 

or  I  r-1  <7-» 

i-1  i 

where  ni,  n2,  03... nk  represent  the  number  of  cycles  at  specific  overstress 
levels;  Nj.,  N2,  N3...Nk  thenlife  cycles  to  failure  at  these  levels  as  read 
from  the  S-N  curve,  and  E  — -  is  the  fraction  of  the  total  life  consumed  by 
l  cycles.  Ki  ' 

The  equivalent  life  (the  total  life)  of  a  part  (N  )  under  a  spectrum 
of  stresses  may  be  founu  by  proportion:  equ 


1 


i-k 

I 

i-1 


i-k 

I 

i-1 


Suppose,  for  example,  there  are  three  stress  levels,  90,  70,  and  SO  ksi, 
in  a  given  spectrum.  With  reference  to  the  curve  in  Figure  7.3,  1/(6  r.  10^) 
la  the  fraction  of  life  consumed  by  each  90  ksi  stress  cycle,  1/(5  x  10^)  by 
each  70  ksi  cycle,  1/(8  x  10^)  by  each  55  ksi  cycle,  etc.  Using  Equation  (7.2): 

N«qu  -  — j — - - —  -  1.5  x  105  cycles 

6x10^  +  5x105“  +  8x10^ 
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Stress 


Figure  7.2  Conversion  of  a  Stress  Spectrum  to  An  Equivalent  Stress 


Figure  7.3  A  Modified  Goodman  Diagram 


Thus,  the  life  of  the  part  under  the  above  spectrum  of  stresses  will  be 
equivalent  to  a  life  of  1.5  x  105  cycles.  The  completely  reversed  stress  equi¬ 
valent  to  this  life  is  (from  Figure  7.5)  75  ksi.  Hence,  the  damage  that  the 
part  accumulates  due  to  the  above  spectrum  of  varying  stress  amplitude  will  be 
the  same  as  if  stress  cycles  of  constant  amplitude  equal  to  Sequ  (in  this  case 
75  ksi)  were  imposed  for  Nequ  (1.5  x  10^  cycles)..  Thus,  the  spectrum  of  stres¬ 
ses  can  be  replaced  by  a  single  stress. 

Miner's  Rule,  as  stated  in  Equation  (7.1) ,  gives  one  (1.0)  as  the  cri¬ 
terion  for  failure.  Miner's  original  tests  showed  that  the  value  for  the  sum¬ 
mation  in  Equation  7.1  actually  varied  between  0.61  and  1.45.  His  more  recent 
dataCU)  gives  a  range  of  0.7  to  2.2.  Other  sources^'  quote  a  range  as  high 
as  0.18  to  23.0.  In  view  of  all  this  scatter,  it  is  generally  agreed  that  the 
value  of  one  (1.0),  originally  proposed  by  Miner,  is  probably  the  best  overall 
estimate  that  can  be  made  at  this  time. 

7.2.2  Corten-Dolan's  Rule 


The  application  of  this  rule  in  converting  the  stress  spectrum  to  a 
single  equivalent  stress  (Sequ)  is  identical  to  that  of  Miner's  rule,  except 
that  the  S-N  curve  used  to  obtain  the  life  values  Ni,  N2,...Njt  is  modified. 
This  modification  is  done,  as  shown  in  Figure  7.6,  by  changing  the  slope  of 
the  S-N  curve.  A  line  is  drawn  with  an  inverse  slope,  d,  passing  through  the 
point  Ni  on  the  S-N  curve,  having  the  maximum  stress  amplitude  (in  this  case, 
S^)  occurring  in  the  stress  spectrum.  This  new  line  is  known  as  the  Corten- 
Dolan  line. 

From  available  data^  ^  it  appears  that  for  specimens  having  no 
stress  concentration  (Kf  ■  1),  the  value  of  d/d’  =  0.8  is  a  reasonable  esti¬ 
mate.  A  recent  study  by  Harris  and  Lipson'^ '  indicates  that  when  stress 
concentrations  are  present  the  following  relationship  can  be  used 

d/d'  -  (0.73  +  0.07  Kf)  (7.3) 

This  can  be  graphically  expressed  as  in  Figure  7.7.  It  will  be  noted  that  if 
Kf  ■  3.5,  d/d'  *  1,  and  this  becomes  equivalent  to  the  criterion  obtained  from 
Miner's  Rule. 
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Stress,  ksl  (log  scale) 


Stress,  ksi  (log  scale) 


Life,  cycles 


Figure  7.7  Corten-Dolan's  Lines  for  Various  Stress*,  .w 
Concentration  Factors,  According  to  Harris  and  Lipson^1  ' 
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SECTION  8  THE  INTERFERENCE  OF  A  STRESS  DISTRIBUTION 
WITH  A  STRENGTH  DISTRIBUTION 

After  the  strength  distribution  and  the  stress  distribution  are  deter¬ 
mined  (Sections  6  and  7  respectively) ,  the  two  ate  compared  and  the  percent 
interference  is  determined  as  discussed  in  Section  2,  Section  5,  and  in  detail 
in  Section  9.  For  a  given  strength  distribution  the  percent  interference  will 
depend  on  the  distribution  of  the  equivalent  stress  Sequ.  A  search  through 
literature  and  other  sources  produced  a  considerable  amount  of  data  leading 
to  the  distribution  of  strength  but  very  little  information  on  the  distribu¬ 
tion  of  stress. 

In  some  engineering  applications  there  is  very  little  scatter  in 
stresses.  This  leads  to  a  stress  distribution  with  a  standard  deviation  equal 
to  zero.  This  distribution  can  be  represented  by  a  straight  line,  as  in 
Figure  8.1,  and  the  interference  can  be  determined  as  shown. 

For  a  given  S6qU,  interference  may  increase  or  decrease  if  the  life  to 
which  the  components  are  designed  is  changed.  .This  is  shown  in  Figure  8.2, 
and  in  terms  of  S-N  diagram  in  Figure  8.3.  The  shape  of  the  distribution 
curve  in  Figure  8.2  is  different  from  those  in  Figure  8.3  because  the  former 
are  plotted  on  a  linear  scale  and  the  latter  on  a  log-log  scale. 

In  those  engineering  applications  where  the  scatter  in  stresses  is 
appreciable,  the  above  approach  will  obviously  not  apply.  On  the  basis  of 
past  experience,  in  the  present  investigation  the  stress  distribution  (SeqU) 
was  assumed  to  be  normal  and  the  range  of  standard  deviations  to  be  not  less 
than  .Oiu  and  not  more  than  ,10p  where  u  is  equal  to  SeqU.  The  resulting 
interference  is  represented  qualitatively  in  Figure  3.4. 

Examples  of  design  problems  employing  this  method  are  giver  in 
Section  9. 


Stress 


S 

equ 

Stresi*  or  Strength,  ksi 


Figure  8.1  Interference  with  Standard  Deviation 
of  Stress  equal  to  Zero 


Stress ,ksi  (log  scale) 


I 


S-N  curves  with 


lift*  cycles 


Figure  8.3  S-N  Diagram  Representing  the  Dependence 
of  Interference  on  Life 
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Frequency  of  Occurrence 


Stress  (Strength) 

Figure  8.4  Interference  of  Se(JU  with  Strength  Distribution 
for  Verioue  Values  of  Standard  Deviation 
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SECTION  9  APPLICATION  OF  INTERFERENCE  THEORY 
TO  DESIGN  PROBLEMS 

Once  the  parameters  of  the  strength  distribution  (X  ,Q,b;  or  p,a;  or  $, 
M)  and  the  stress  distribution  (p=SeqU  and  a=Kp,  where  K  represents  a  frac¬ 
tion  of  the  average  stress)  are  determined,  as  shown  in  Sections  6  and  7  res¬ 
pectively,  the  percent  interference  and  consequently  the  reliability  (Relia¬ 
bility*  1-  Interference)  can  be  computed  with  the  aid  of  the  tables  in  Appen¬ 
dix  2.  The  specific  steps  to  be  taken  in  solving  particular  design  problems 
are  illustrated  by  the  following  examples. 

9.1  EXAMPLE  PROBLEM  NO.  1  (EFFECT  OF  STRESS  CONCENTRATION  ON  RELIABILITY) 

The  present  design  for  an  aircraft  part  specifies  the  following: 


Material : 


Design  Life: 


2024  Aluminum  S  =70  ksi,  S  =50  ksi 
u  y 

10^  cycles 


Type  of  Loading: 


Axial 


Size: 


.125"  sheet 


Surface  Finish: 


Electropolished 


Stress  Concentration  Factor:  Kt  *  2.4  (Hole) 


Operating  Temperature: 


Room  Temperature 


An  alternative  design  calls  for  reducing  the  stress  concentration  at 
the  critical  section,  at  additional  cost,  by  milling  an  edge  notch  instead  of 
drilling  a  hole.  This  alternative  design  reduces  the  theoretical  stress  con¬ 
centration  factor,  Kt,  to  1.5.  Will  the  resultant  reduction  in  stress  concen¬ 
tration  affect  the  reliability  of  the  part? 

9.1.3  Strength  Parameters 

The  first  step  is  to  determine  the  strength  distribution  in  terms  of  its 
parameters  corresponding  to  the  design  conditions.  From  the  table  on  Page  194 
(Code  No.  100),  at  10^  cycles  the  strength  distribution  of  the  present  design 
is  found  to  be: 


Strength  Distribution: 
Strength  Parameters: 


Weibull 

X  -  13.76  ksi 
o 

b  =  1.808 
e  -  16.31  ksi 


For  the  alternative  design,  the  strength  distribution  at  10  cycles  is  found 
from  the  table  on  Page  194  (Code  No.  102  )  to  be: 


133 


Strength  Distribution:  Smallest  Extreme  Value 
Strength  Parameters:  8  ■  .8348 

M  -  19.98  ksi 


9.1.2  Stress  Parameters 


In  order  to  determine  the  parameters  of  the  stress  distributions  (y  ■ 
SeqU  and  o  ■  Ky),  a  prototype  of  each  design  was  instrumented  and  the  stress 
spectrums  were  recorded  as  shown  in  Columns  1  and  2  of  Tables  9.1  and  9.2.  To 
determine  the  equivalent  stress  (SeqU) ,  Miner's  rule  is  used.  From  the  S-N 
curves  of  the  material  (Figure  9.1  and  9.2)  the  number  of  cycles  to  failure, N, 
corresponding  to  the  stresses  in  Column  1  of  Tables  9.1  and  9.2  are  determined. 
These  are  shown  in  Column  3  of  Tables  9.1  and  9.2. 


Using  Miner's  Rule, 
data  in  Table  9.1,  NeqU  is 

N 

equ 


as  expressed  in  Equation  (7.2),  and  the  tabulated 
determined  for  the  present  design: 


1  x 


In^ 

ynl 


V  '  1  *  255?7a5=S  ‘  4'05  *  1()6 

From  the  S-N  curve  (Figure  9.1)  the  stress  corresponding  to  Nequ  ■  4.05  x  106 
cycles  is  found  to  be  Sequ  -  13.7  ksi.  Hence,  the  application  of  a  completely 
reversed  stress  of  13.7  ksi  can  be  substituted  for  the  recorded  stress  spec¬ 
trum  in  Columns  1  and  2  of  Table  9.1.  Similarly,  using  the  data  from  Table 
9.2  for  the  alternative  design: 

"equ  ■  1  x  ersoJlO-S  ’  1>H  *  1<)7  “l*1" 

From  the  S-N  Curve  (Figure  9.2)  the  stress  corresponding  to  NeqU  ■  1.51  x  107 
cycles  is  found  to  be  Sequ  *  13.5  ksi. 

/ In  some  engineering  applications,  the  scatter  in  operating  stresses  is 
very  small,  therefore,  the  standard  deviation  of  the  equivalent  stress  can  be 
assumed  to  be  zero.  In  those  engineering  applications  where  the  scatter  in 
stress  is  appreciable,  the  standard  deviation  lies  in  the  range: 

O.Oly  *  o  *  O.lOy  . 

In  the  absence  of  any  specific  information,  an  average  value  of  o  ■  .05u  can 
probably  be  assumed.  For  the  present  problem,  interference  will  be  calculated 
for  the  two  cases:  o  -  0  and  a  ■  .05y. 


Thus,  the  stress  parameters  for  the  present  design  are: 

Case  1.  u  •  S#qu  -  13.7  ksi  and  o  -  0  and, 

Case  2.  y  -  13.7  ksi  and  a  -  .05y  -  (.05) (13.7)  -  .685  ksi. 
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PRESENT  DESIGN 


1  Stress  Spectrum 

Miner ' 8 

Rule  Data 

Completely*  Reversec 

Number  of 

Cycles  to 

ni 

Axial  Stresses,  ksl 

Occurrences,  Ni 

Failure,  Ni 

Ni 

1 

2 

3 

4 

11.7 

240 

1.5  x  107 

16.0  x  10"6 

12.5 

217 

9  x  106 

24.1  x  10“6 

13.0 

176 

6.2  x  106 

28.4  x  10"6 

13.8 

150 

3.9  x  106 

38.4  x  10~6 

14.1 

110 

2.9  x  106 

37.9  x  10’6 

14.9 

75 

1.8  x  106 

41.6  x  10~6 

15.7 

52 

1.15  x  106 

45.2  x  10“6 

15.9 

20 

1.05  x  106 

19.0  x  10“6 

16.0 

5 

1.0  x  106 

5.0  x  10"6 

"  1035 

£%  •255-7  *  10'6 

Table  9.1  Stress  and  Life  Data  for  Miner's  Rule 


*Actually,  the  stress  was  not  completely  reversed.  It  was  reduced  with  the 
aid  of  Goodman  Diagram  to  a  completely  reversed  stress  using  the  procedure 
given  In  Section  7.2. 


ALTERNATIVE  DESIGN 


|  Stress  Spectrum 

Miner’s 

Rule  Data 

Completely*  Reversed 

Number  of 

Cycles  to 

ni 

Axial  Stresses,  ksl 

Occurrences,  ni 

Failure,  Ni 

Ni 

1 

2 

3 

4 

li.7 

240 

5.1  x  107 

4.70  x  10-6 

12.5 

217 

3.0  x  107 

7.24  x  10"6 

13.0 

176 

2.1  x  107 

8.39  x  10~6 

13.8 

150 

1.3  x  107 

11.52  x  10~6 

14.1 

110 

1.1  x  107 

10.00  x  10"6 

14.9 

75 

7.2  x  !06 

10.40  x  10~6 

15.7 

52 

4.8  x  106 

10.81  x  10~6 

15.9 

20 

4.4  x  106 

4.55  x  10'6 

16.0 

5 

4.2  x  106 

1.19  x  10"6 

[nt  -  1035 

n.  f, 

l  ~  -  68.80  x  10_o 

Table  9.2  Stress  and  Life  Data  for  Miner's  Rule 


*Actually,  the  stress  was  not  completely  reversed.  It  was  reduced  with  the 
eld  of  a  Goodman  diagram  to  a  completely  reversed  stress  using  the  procedure 
given  In  Section  7.2. 
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Ufa, cycles 

Figure  9.1  S-N  Relationship  for  the  Present  Design 


JNOT  REPRODUCIBLE 


Life,  cycles 

Relationship  for  the  Alternative  Design 


Similarly,  the  stress  parameters  for  the  alternative  design  are: 

1.  y  ■  S6qU  -  13.5  ksi  and  a  *  0 
and  2.  y  -  13.5  ksi  and  o  =  .05y  -  (.05)  (13.5)  -  .675  ksi. 

9.1.3  Percent  Interference 

Once  the  strength  and  stress  distribution  parameters  are  determined, 
the  Percent  Interference  and  the  Reliability  (%Reliability  *  100-%  Interfer¬ 
ence)  can  be  determined. 


Case  (1)  when  stress  standard  deviation  (o)  *  0: 

For  this  case  where  stress  is  Normally  distributed  with  y  *  Sequ  and 
0,  the  percent  interference  can  be  determined  as  follows: 

Present  Design 

For  Weibull  distributed  strength, 

.x-Xo.b 

"  vg-Y  ) 

Interference  F(x)  ■  1  -  e  ;  0  «  shaded  area  under  Figure  8.1 

where  x  -  Sequ  -  13.70  ksi 
Xq  -  13.76  ksi 
b  -  1.808 
9  -  16.31  ksi 


F(x)  -  1  -  e 


.13. 70-13. 76v 1.308 
”^16.31-13.76* 


Percent  Interference  -  0% 

It  can  also  be  seen,  from  Figure  A-2,1  in  Section  A-2.2.1,  that  when  Seq  is 
less  than  X0  there  la  no  Interference.  In  the  case  when  SeqU  la  greater  qthan 
Xq ,  the  Percent  Interference  can  be  determined  as: 


Find, 


S  -X  b 
I  equ  o 


9  -  X, 


Corresponding  to  the  value  of  X  read  Interference  F(x)  from  the  table  In  Section 
A-2.2.1,  Pages  371  or  372. 

Alternative  Design 

For  Smallest  Extreme  Value  distributed  strength, 

+8 (x-M) 

Interference  F(x)  ■  1  -  e  »  shaded  area  under  Figure  8.1 
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where  x  -  S 


13.50  ksi 


F(x) 


equ 

6  =  .8348 
M  =  19.93 

.8348(13.50-19.98) 

=  1  -  e 


=  1  - 


-5.404 

-e 

e 


=  1  - 


-.0045 

e 


=  .0045 


Percent  Interference  =  .45% 

This  can  also  be  read  directly  from  the  table  in  Section  A-2.5.1,  Page  419. 

Find  X  =  -3(Sequ-M)  =  -.8348  (13.50-19.98) 

-  5.404 

Corresponding  to  X  =  5.404  read  interference  F(x)  =  .0045  from  the  above 
table.  Therefore,  Percent  Interference  =  .45%. 


Case  (2)  when  stress  standard  deviation  (a)  ^  0.*  (in  this  case  a-.05vi) 

For  this  case  where  stress  is  Normally  distributed  with  a  i  0,  addi¬ 
tional  interference  parameters  must  be  calculated  to  determine  the  percent 
interference. 


Present  Design 


Strength  (Weibull) 


Stress  (Normal) 


X  -  13.76  ksi 
o 

b  -  1.308 
e  «  16.31  ksi 


U  ■  SeqU  ■  13.70  ksi 
a  ■  0.05u 
■  0.685  ksi 


For  the  above  data,  interference  parameters  C,  A,  and  B(x)  to  be  used  in  the 
interference  tables  were  computed. 


o-x 

0 

16.31-13.76 

3.72 

a 

0.685 

\>z 

13.76-13.70 

.0875 

0 

0.685 

B(x)  -  b  -  1.808 

From  the  tables  in  Section  A-2.2.2  on  Pages  588  and  590,  the  interference  val¬ 
ue  corresponding  to  these  parameters  is  read  using  linear  interpolation  (for 
a  higher  order  anti  more  accurate  interpolating  technique,  see  Section  A-2.1). 


Interference  =  .0398 
or  Percent  Interference  *  3.98% 


Alternative  Design 

Strength  (S.E.V.)  Stress  (Normal) 

8  =  .8348  p  =  Sequ  =  13.50  ksi 

M  =  19.98  ksi  a  =  .05u  ~  .675  ksi 

For  the  above  data,  interference  parameters  a  and  y  to  be  used  in  the  inter¬ 
ference  tables  are  computed. 

a  *  6a  *»  (.8348)  (.675)  =  .564 
y  =  8(SeqU-M)  *(.8348) (-6.48)  *  -5.4 

From  the  table  in  Section  A-2 .5 -2  on  Page422,  the  interference  value  correspond¬ 
ing  to  these  parameters  is  read  using  linear  interpolation  (for  a  higher  order 
and  more  accurate  interpolating  technique,  see  SectionA-2.l). 

Interference  *  .0059 
or  Percent  Interference  *  .59% 

A  summary  of  the  results  is  given  in  Table  9.3. 

Thus,  the  reduction  in  the  stress  concentration  factor  should,  in 
effect,  result  not  only  in  lower  stress  but  also  in  somewhat  higher  reliability. 
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9.2  EXAMPLE  PROBLEM  NO.  2  (EFFECT  OF  TEMPERATURE  ON  RELIABILITY) 

A  part  In  service  has  been  performing  satisfactorily  under  the  follow¬ 
ing  conditions : 


Material : 


T1-4A1-3MO-1V;  .125'  •heat;  Su  •  16.0  ksi 
Axial  8y  “  13,2  kil 


Loading:  Axial  ? 

Surface  Finish:  Ground 

Stress  Concentration  Factor 

at  Critical  Section:  Kt  *  1.0 

Operating  Temperature:  400*F 

C 

Reliability  Requirement:  99.99%  at  5  x  10  cycles 

Stress  Parameters  for  the  y  -  28.2  ksi 

Critical  Section:  o  -  1.41  ksi  determined) 

If  the  efficiency  of  the  system  can  be  increased  by  increasing  the  operating 
temperature,  how  will  the  reliability  of  the  component  part  be  affected  if  the 
temperature  is  increased  to  600*F,800*F,  and  the  stress  parameters  remain  the  same? 


9.2.1  Strength  Parameters 


From  the  table  on  Page  303  (Code  Nos  •  586  and  587)  the  strength 
parameters  for  600#F  end  800*7  at  5  x  106  cycles  ere: 

600*F 

Strength  Distribution:  Largest  Extreme  Value 
Strength  Parameters:  8  ■  .2066 


40.79  ksi 


Mil 


Strength  Distribution:  Largest  Extreme  Value 
Strength  Parameters:  6  ■  .2836 

M  -  28.74  ksi 
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9.2.2  Stress  Parameters 


The  stress  parameters  have  been  previously  determined  and  are: 

Stress  Distribution:  Normal 

Stress  Parameters:  y  =  SeqU  =  28.2  ksi 

a  =  1.41  ksi 


9.2.3  Percent  Interference 

With  the  strength  and  stress  distribution  parameters  established,  the 
percent  interference  can  be  determined  at  the  various  temperatures: 

600aF 

Strength  Distribution:  Largest  Extreme  Value  (B,M) 

Stress  Distribution:  Normal  (y,cr) 

Interference  Parameters: 

o  ■  So  *  .2066  (1.41)  *  .292 

y  =  S(y-M)  -  .2066  (28.2-40.79)  -  -2.6 

Using  parameters  (a,y),  the  interference  at  5  x  10^  cycles  is  obtained 
by  linear  interpolation  from  the  table  in  Section  A-2.4.2  on  Page  4ll(For  a 
more  accurate  interpolation  method,  see  Section  A-2.1). 

Interference  ■  .00092 
Percent  Interference  -  .092% 

Percent  Reliability  ■  (100-%  Interference) 

-  99.908% 

800*F 

Strength  Distribution:  Largest  Extreme  Value (B,M) 

Stress  Distribution:  Normal  (u,o) 

Interference  Parameters : 

o  -  8a  -  (.2836) (1.41)  -  .400 
y  -  8(u-M)  -  .2836  (28.20-28.74)  -  -.153 
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Using  parameters  (a,y) ,  the  interference  at  5  x  10°  cycles  is  obtained 
by  linear  interpolation  from  the  table  in  Section  A-2.4.2  on  Page  411  (For  a 
more  accurate  interpolation  method,  see  Section  A-2.1). 

Interference  =  .33426 
Percent  Interference  =  33.426% 

Percent  Reliability  =  (100-%  Interference) 

=  66.574% 

A  summary  of  the  results  is  given  in  Table  9.4. 

Thus,  an  increase  in  the  operating  temperature  from  600°F  to  800°F 
results  in  an  appreciable  reduction  in  reliability. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


A  method  was  developed  for  employing  Stress-Strength  Interference  Theory 
as  a  practical  engineering  tool  to  be  used  for  designing  and  quantita¬ 
tively  predicting  the  reliability  of  mechanical  parts  and  components,  made 
of  non-ferrous  materials,  subjected  to  mechanical  loading. 

This  method  is  based  on  considerable  empirical  data  gathered  (Appen¬ 
dix  1)  and  it  also  has  sound  theoretical  basis  (Appendix  3).  This  method 
eliminates  the  concept  of  a  Factor  of  Safety  and  substitutes  Percent  In¬ 
terference  (Probability  of  Failure).  Tables  of  interference  values  are 
given  in  Appendix  2  for  a  variety  of  stress  and  strength  conditions. 

Although  a  great  deal  of  data  was  gathered  and  analyzed  in  the  course  of 
the  present  study,  no  data  were  found  to  permit  the  establishment  of  con¬ 
fidence  intervals  on  the  probability  of  interference. 

In  our  past  investigation^) ,  because  of  the  limitations  of  the  graphical 
approach  employed,  this  method  was  limited  to  the  following  three  cases: 

Stress  Distribution  Strength  Distribution 

Normal  Normal 

Normal  tfeibull 

Weibull  Weibull 

In  the  present  investigation,  based  on  the  computer  approach,  the  method 
was  extended  to  the  following  three  additional  cases: 

Stress  Distribution  Strength  Distribution 

Normal  Weibull  (Extended-see  Sec¬ 

tion  4.2) 

Normal  Largest  Extreme  Value 

Normal  Smallest  Extreme  Value 

5.  For  these  combinations  of  stress  and  strength  distributions,  interference 
values  were  calculated  by  an  Integral  Method  and  these  are  tabulated  In 
Appendix  2.  These  values  were  found  to  be  highly  non-linear.  Therefore, 
for  the  purpose  of  interpolating  these  values  in  a  given  table  or  between 
the  tables,  a  higher  order  interpolation  (Section  A-2.1)  should  be  used. 

6.  A  computer  approach  (referred  to  above)  was  developed  for  determining  the 
statistical  distribution  function  of  fatigue  strength  at  a  given  life. 

When  the  raw  data  (conventional  S-N  type  data)  were  fed  into  this  computer 
program,  the  computer  printed  the  correlation  coefficient  (a  degree  of 
fit)  and  the  parameters  for  Normal,  Weibull,  Largest  Extreme  Value  and 
Smallest  Extreme  Value  distributions,  for  any  given  sec  of  strength  date. 
This  method  has  the  advantage  of  high  accuracy  and  time  saving  over  the 
graphical  method  used  in  the  previous  investigation. (1) 

7.  The  effect  of  type  of  loading,  surface  finish,  stress  concentration,  haat 
treatment,  Temperature,  environment  and  other  factors  (Section  4.4)  on  the 
statistical  distribution  function  of  fatigue  strength  of  various  non-fer- 


rous  materials  was  studied.  These  effects  were  expressed  in  terms  of  the 
distribution  parameters  which  were  tabulated  (Appendix  1,  pages  155-350) 
and  the  representative  ones  plotted  (Section  6.7,  pages  53-117). 

8.  In  most  cases,  the  effect  of  the  factors  listed  under  item  7  above  was  to 
change  the  distribution  function.  That  is,  for  a  given  material  and  sur¬ 
face  finish  the  distribution  function  might  be  Weibull  but  for  the  same 
material  and  different  surface  finish,  the  distribution  could  change  to 
Largest  Extreme  Value.  (This  refers  to  the  best  fitting  distribution  for 
each  particular  material  and  set  of  conditions,  as  described  in  Section 
6.6).  The  only  exception  is  the  effect  of  life.  Life  was  found  to  retain 
the  same  distribution  and  change  only  the  parameters  of  that  distribution 
(see  the  basic  assumption  given  in  Section  6.2  on  page  31). 

9.  For  the  non-ferrous  materials  studied,  the  effect  of  life  on  the  fatigue 
strength  distribution  parameters  is  shown  in  the  table  below: 


Fatigue  Strength  Distribution  Effect  of  Life  on  Distribution 

Distribution _ Parameters _ Parameters _ 


Weibull 

Xo 

6 

b 

Linearly  decreases  with  life  on 
log-log  scale 

Linearly  decreases  with  life  on 
log-log  scale 

Increases  or  decreases  with  life 
depending  on  conditions 

li 

Linearly  decreases  with  life  on 

Normal 

log-log  scale 

c 

Decreases  with  life 

M 

Linearly  decreases  with  life  on 

Largest  Extreme 

log-log  scale 

Valua 

B 

Increases  with  life 

M 

Linearly  decreases  with  life  on 

Smallest  Extreme 

log-log  Beale 

Value 

P 

Increases  with  life 

10.  As  to  the  problem  of  stress,  the  stress  distribution  was  investigated  in 
more  detail  than  in  the  previous  study(^)  but  no  additional  information 
was  located. 

11.  In  order  to  verify  the  validity  of  the  Interference  Technique  developed 
in  the  present  study,  it  should  be  checked  against  an  actual  life  situ¬ 
ation.  That  is,  percent  Interference  should  be  computed  for  an  actual 
engineering  problem.  These  results  then  should  be  compared  with  actual 
service  failures. 
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APPENDIX  1  TABLES  OF  THE  DISTRIBUTION  PARAMETERS  OF  FATIGUE  STRENGTH 


A-l.l  INDEX  TO  DISTRIBUTION  PARAMETERS 
FOR 

VARIOUS  NON-FERROUS  MATERIALS 

Distribution  Parameters 
can  be  found: 


Material  Page  No. 

ALUMINUM  ALLOYS 

1.  1100  Aluminum  175 

2.  5  Mg  Aluminum  176 

3.  355  Aluminum  176 

4.  7.5  Zn-2.5  Mg  Aluminum  177 

5.  M-257  Aluminum  179 

6.  M-276  Aluminum  180 

7.  2014  Aluminum  182 

8.  2014  Alclad  Aluminum  187 

9.  R-303  Aluminum  189 

10.  2020  Aluminum  189 

11.  2024  Aluminum  (Sheets,  Bars,  Extrusions,  Forgings)  190 

12.  2024  Alclad  Aluminum  202 

13.  2025  Aluminum  203 

14.  2026  Aluminum  205 

15.  2219  Aluminum  206 

16.  2618  Aluminum  207 

17.  5052  Aluminum  210 

18.  5056  Aluminum  2 1C 

19.  5083  Aluminum  211 

20.  5086  Aluminum  211 


21.  5154  Aluminum 


212 


22.  6061  Aluminum 


212 


23. 

5456  Aluminum 

213 

24. 

7001  Aluminum 

214 

25. 

7039  Aluminum 

217 

26. 

7075  Aluminum  (Plates  and  Sheets, 

Bars,  Extrusions  and  Forgings) 

218 

27. 

7075  Alclad  Aluminum 

228 

28. 

7076  Aluminum 

229 

29. 

7106  Aluminum 

230 

30. 

7079  Aluminum 

231 

31. 

7178  Aluminum 

234 

COBALT  ALLOYS 

32. 

Stellite-31 

237 

33. 

S-816 

241 

COPPER  ALLOYS 

34. 

Pure  Copper 

244 

35. 

70-30  Brass 

245 

36. 

Cu-7.3  A1  Bronze 

247 

37. 

5.6  A1  Bronze 

247 

38. 

Al-Ni  Bronze 

248 

39. 

Phosphor  Bronze 

248 

40. 

Muntz  Metal  (Brass) 

249 

41. 

Cu-6.5  Al-2.4  Fe 

249 

MAGNESIUM  ALLOYS 

42. 

2.5  Al  Magnesium 

251 

43. 

AZ31A  Magnesium 

252 

44. 

AZ31B  Magnesium 

252 

253 


45. 

AZ61A  Magnesium 

253 

46. 

AZ80A-F  Magnesium 

254 

47. 

AZ81  Cast  Magnesium 

255 

48. 

HM21  Magnesium 

255 

49. 

ZK60A  Magnesium 

256 

NICKEL  ALLOYS 

50. 

"A"  Nickel 

258 

51. 

FS-27  Nickel  Base  Cermet 

259 

52. 

Hastelloy  C 

259 

53. 

Hastelloy-R235 

260 

54. 

Incoloy-901  (AMS-5560A) 

262 

55. 

Inconel 

263 

56. 

Inconel  713C 

264 

57. 

Inconel-X  (AMS-5667) 

265 

58. 

Inconel-718 

266 

59. 

Waspalloy 

269 

60. 

Inconel-751  (INCO-X550) 

270 

61. 

INOR-8 

271 

62. 

Monel 

271 

63. 

Nicro-Tung 

272 

64. 

Nlmonlc-95 

274 

65. 

Rene-41 

276 

66. 

Nlmonic-80 

277 

67. 

Udimet-500 

278 

68. 

Udimet-650 

280 

69. 

6  Mo-Waspalloy 

281 
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TITANIUM  ALLOYS 


70.  T1-A55  284 

71.  T1-75A  288 

72.  T1-150A  290 

73.  Ti-0.2  02  291 

74.  Ti-0.2  C  292 

75.  Ti-N2  293 

76.  Ti-0.07  N2-0.2  02-0.2  C  293 

77.  K-151A  Ti  Base  Cermet  294 

78.  K-162B  Ti  Base  Cermet  ,  294 

79.  K-183A  Ti  Base  Cermet  295 

80.  Ti-4  Al-3  Mo-1  V  296 

81.  Ti-4  Al-4  Mn  307 

82.  Ti-5  Al-2.5  Sn  308 

83.  Ti-5  Al-2.5  Sn-0.07  N2  308 

84.  Ti-5  Al-2.5  Sn-0.2  02  309 

85.  Ti-5  Al-2.5  Sn-.2  C  309 

86.  Ti-5  Al-2.5  Sn-.07  N2-.2  02-.2  C  3i0 

87.  Ti-6  A1  311 

88.  Ti-6  Al-4  V  313 

89.  Ti-6  Al-4  V-0.07  N2  315 

90.  Ti-6  Al-4  V-0.2  02  316 

91.  Ti-6  Al-4  V-0.2  C  316 

92.  Ti-6  Al-4  V-.07  N2-.2  02-.2  C  317 

93.  Ti-7  Al-4  Mo  317 

94.  Ti-7  Al-3  Mo  318 
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95 


319 


.  Ti-8  Al-1  Mo-1  V 

96.  Ti-3  Mn-0.2  02  320 

97.  Ti-3  Mn-0.07  N2  321 

98.  Ti-3  Mn-0.2  C  32i 

99.  Ti-3  Mn  Complex  322 

100.  Ti-3  Mn-.07  N2-.2  02-.2  C  325 

101.  Ti-4  Mn  325 

102.  Ti-8  Mn  326 

103.  Ti-30  Mo  328 

104.  Ti-13  V-ll  Cr-3  A1  329 

105.  Ti-16  V-2.5  Sn  334 

OTHER  NON-FERROUS  ALLOYS 

106.  Berilco  #25  347 

107.  Commercially  Pure  Lead  347 

108.  0.5  Ti  Molybdenum  348 

109.  Commercially  Pure  Tantalum  349 

110.  10  W-Tantalum  350 

111.  Pure  Tungsten  350 


ICO 


A-1.2  distribution  functions  and  their  parameters  used  in  tables 

A-I.2.1  Weibull  Distribution  |\ 
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b  i«  the  Weibull  Slop.  (elope  of  the  lineal  fora  of  the  distribution). 

«  the  ch.reet.ri. tic  .treagth  (.here  63.2!  of  th.  population  ha*e  fatigue  strengths  less  then  or 

equal  co  chis  value). 

X  is  the  lower  bound  of  fatigue  strength. 


A-l.2.2  Normal  Distribution 
Density  Function 


population  mean  (an  average  value), 
standard  deviation  (a  measure  of  scatter) . 

standardized  normal  variate  (number  of  standard  deviations  from  the  mean). 


s  intensity  function  (slope  of  the  linear  form  cf  the  distribution). 
!  mode  (the  most  probable  value) . 


Smallest  Extreme  Value  Distribution 
Density  Function 


intensity  function  (slope  of  the  linear  form  of  the  distribution) . 
mode  (the  most  orobable  value). 


A-1.3  list  of  abbreviations  and  symbols  used  in  tables 
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A-1.4  COMPOSITION  OF  THE  NON-FERROUS  MATERIALS 


I.  ALUMINUM  ALLOYS 


Alloy 

Percent  Composition 

1.  1100  ALUMINUM 

99.7  A1 

Nominal 

2.  5  Mg  ALUMINUM 

.21  Fe,  .58  Si,  5.02  Mg,  .26  Mn, 
Balance  A1 

Actual 

3.  355  ALUMINUM 

5  Si,  1.3  Cu,  .5  Mg,  Balance  A1 

Nominal 

4.  7.5  Zn  -  2.5  Mg  ALUMINUM 

7.5  Zn,  2.5  Mg,  Balance  A1 

Nominal 

5.  M-257  ALUMINUM 

Not  Reported 

6.  M-276  ALUMINUM 

Not  Reported 

7.  2014  ALUMINUM 

4.46  Cu,  .43  Fe,  .95  Si,  .8  Mn, 

.44  Mg,  .072  Zn,  Balance  A1 

Nominal 

8.  2014  ALfLAD  ALUMINUM 

4.46  Cu,  .43  Fe,  .95  Si,  .8  Mn, 

.44  Mg,  .072  Zn,  Balance  A1 

Nominal 

9.  R-303  ALUMINUM 

6.4  Zn,  2.5  Mg,L2  Cu,  Balance  A1 

Nominal 

10,  2020  ALUMINUM 

Not  Reported 

11.  2024  ALUMINUM 

4.5  Cu,  1.5  Mg,  .6  Mn,  Balance  A1 

Nominal 

12.  2024  ALCLAD  ALUMINUM 

4.5  Cu,  1.5  Mg,  .6  Mn,  Balance  A1 

Nominal 

13.  2025  ALUMINUM 

4.28  Cu,  .36  Fe,  .77  Mn,  .76  Si, 
Balance  A1 

Actual 

14.  2026  ALUMINUM 

.5  Mg,  .75  Si,  .75  Mn,  4.5  Cu, 
Balance  A1 

Actual 

15.  2219  ALUMINUM 

.2  Si,  .3  Fe,  6  Cu,  .3  Mn,  .02  Mg, 
.1  Zn,  Balance  A1 

Actual 

16.  2618  ALUMINUM 

2.4  Cu,  .25  Si,  1  Fe,  1.5  Mg, 

.08  Ti,  .1  Ni,  Balance  A1 

Actual 

17.  5052  ALUMINUM 

2.5  Mg,  .25  Cr,  Balance  A1 

Nominal 

18.  5056  ALUMINUM 

5.2  Mg,  .1  Ms,  .1  Cr,  Balance  A1 

Nominal 

19.  5083  ALUMINUM 

.27  Fe,  .15  Si,  .06  Cu,  .76  Ma, 

4.48  Mg,  .01  Zn,  .09  Cr,  .01  Ti, 

Balance  A1  Actual 
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Nominal 


Alloy 

Percent  Composition 

20. 

5086  ALUMINUM 

4.0  Mg,  .5  Mn,  Balance  A1 

Nominal 

21. 

5154  ALUMINUM 

3.5  Mg,  .25  Cr,  Balance  A1 

Nominal 

22. 

6061  ALUMINUM 

1.  Mg,  .6  Si,  .25  Cu,  .25  Cr, 

Balance  A1 

Nominal 

23. 

5456  ALUMINUM 

5.  Mg,  .7  Mn,  .15  Cu,  .15  Cr, 

Balance  A1 

Nominal 

24. 

7001  ALUMINUM 

.11  Si,  .12  Fe,  2.1  Cu,  3.03  Mg, 

7.59  Zn,  .2  Cr,  .01  Ti,  Balance  A1 

Actual 

25. 

7039  ALUMINUM 

Not  Reported 

26. 

7075  ALUMINUM 

5.5  Zn,  2.5  Mg,  1.5  Cu,  .3  Cr, 
Balance  A1 

Nominal 

27. 

7075  ALCLAD  ALUMINUM 

5.5  Zn,  2.5  Mg,  1.5  Cu,  .3  Cr, 
Balance  A1 

Nominal 

28. 

7076  ALUMINUM 

.6  Cu,  7.6  Zn,  1.6  Mg,  .6  Mn, 

.1  Ti,  .5  Fe,  .25  Si,  Balance  A1 

Actual 

29. 

7106  ALUMINUM 

.35  Si,  .35  Fe,  .1  Cu,  .2  Mn,  2  Mg, 

4  Zn,  .1  Cr,  Balance  A1 

Actual 

30. 

7079  ALUMINUM 

4.3  Zn,  3.3  Mg,  .6  Cu,  .2  Mn, 

.2  Cr,  Balance  A1 

Nominal 

31. 

7178  ALUMINUM 

6.8  Zn,  2.7  Mg,  2.0  Cu,  .3  Cr, 
Balance  A1 

Nominal 

II. 

OOBALT  ALLOYS 

32. 

STELLITE-31 

4.5  C,  24.8  Cr,  .42  Mn,  10.4  Nl, 

.93  Si,  1.39  Fe,  .01  P,  .010  S, 

.03  Al,  7.26  W,  Balance  Co 

Actual 

33. 

S-816 

.397  C,  1.12  Mn,  .50  Si,  19.42  Cr, 
20.62  Nl,  4.1  Mo,  .018  S,  .012  P, 
2.86  Cb,  4.03  V,  2.99  Fe,  1.09  Ta, 
42.9  Co 

Actual 

III. 

COPPER  ALLOTS 

34. 

PURE  COPPER 

99.98  Cu 

Nominal 

35. 

70-30  BRASS 

70  Cu,  30  Zn 

Nominal 

Alloy 

Percent  Composition 

36. 

Cu-73  A1  BRONZE 

92,7  Cu,  7.3  A1 

Nominal 

37. 

5.6  A1  BRONZE 

94.4  Cu,  5.6  A1 

Nominal 

38. 

Al-Ni  BRONZE 

81.75  Cu,  .71  Mn,  2.86  Fe, 

4.69  Ni,  9.90  A1 

Actual 

39. 

PHOSPHOR  BRONZE 

94.05  Cu,  .18  P,  4.2  Sn,  .05  Pb 

Actual 

40. 

MUNTZ  METAL  (BRASS) 

60.3  Cu,  39.7  Zn 

Actual 

41. 

Cu  -  6.5  A1  -  2.4  Fe 

91.1  Cu,  6.5  Al,  2.4  Fe 

Nominal 

IV. 

MAGNESIUM  AlLOYS 

42. 

2.5  A1  MAGNESIUM 

2.5  Al,  97.5  Mg 

Nominal 

43. 

AZ31A  MAGNESIUM 

3  Al,  1  Zn,  .2  Mn,  Balance  Mg 

Nominal 

44. 

AZ31B  MAGNESIUM 

3  Al,  1  Zn,  .2  Mn,  Balance  Mg 

Nominal 

45. 

AZ61A  MAGNESIUM 

5.5  Al,  1  Zn,  .2  Mn,  Balance  Mg 

Nominal 

46. 

AZ80A-F  MAGNESIUM 

8.5  Al,  .5  Zn,  Balance  Mg 

Nominal 

47. 

AZ81  CAST  MAGNESIUM 

7.5  Al,  .7  Zn,  .15  Mn,  Balance  Mg 

Nominal 

48. 

HM21  MAGNESIUM 

2  Th,  .5  Mn 

Nominal 

49. 

ZK60A  MAGNESIUM 

5.6  Zn,  .66  Zr,  Balance  Mg 

Nominal 

V. 

NICKEL  ALLOYS 

50. 

"A"  NICKEL 

99.35  Ni  +  Co 

Nominal 

51. 

PS-27  CERMET 

34.3  Ti,  50  Ni,  6.2  Cr,  9.5  C 

Actual 

52. 

HASTELLOY  C 

.05  C,  .48  Mn,  .64  Si,  56.6  Ni, 

15.5  Cr,  16.2  Mo,  5.7  Fe,  3.8  W 

Actual 

53. 

HASTELLOY-R235 

.17  C,  .017  S,  2.06  Al,  2.55  Ti, 
5.31  Mo,  15.51  Cr,  21  Co,  9.95  Fe, 
.22  Si,  .03  Mn,  Balance  Ni 

Actual 

54. 

INCOI.OY-901 

12.41  Cr,  6.14  Mo,  J8  Co,  2.67  Ti, 
.13  Al,  43.20  Ni,  Balance  Fe 

Actual 

55.  INCONEL 


02  C,  .13  Mn,  6.81  Fe,  .01  S, 

17  S1,j06  Cu,  75.85  Ni,  15.54  Cr  Actual 


Alloy 

56.  INCONEL-713C 

57.  INCONEL-X 

58.  INCONEL-718 

59.  WASPALLOY 

60.  XNCONEL-751 

61.  INOR-8 

62.  MONEL 

63.  NICRO-TUNG 

64.  NIMONIC-95 

65.  RENE-41 

66.  NIMONIC-80A 

67.  UDIMET-500 

68.  UDIKET-650 


Percent  Composition 

11.9  Cr,  .86  Fe,  .49  Si,  .13  Mn, 

.11  C,  5.0  Mo,  5.6  Al,  .52  Ti, 

2.1  Cb,  Balance  Ni 

.05  C,  .57  Mn,  6.77  Fe,  .29  Si 
.08  Cu,  72.85  Ni,  14.98  Cr,  .72  Al, 
2.54  Ti,  1.12  (Cb  +  Ta) 

.04  C,  .20  Mn,  .20  Si,  52.5  Ni, 

19  Cr,  3  Mo,  18  Fe,  .8  Ti,  5.2  Cb, 
.6  Al 

.1  Si,  .1  Mn,  19.5  Cr,  .55  Fe, 

4.39  Mo,  13.5  Co,  2.87  Ti, 

1.29  Al,  Balance  Ni 

.05  C,  .73  Mn,  .28  Si,  14.92  Cr, 

2.5  Ti,  1.16  Al,  6.59  Fe,  .03  Cu, 
1.03  (Cb  +  Ta),  Balance  Ni 


Not  Reported 
31.2  Cu,  66.9  Ni 

.1  C,  .05  B,  .05  Zr,  12  Cr,  10  Co, 

8  W,  4  Al,  4  Ti,  Balance  Ni 

.08  C,  .36  Si,  .06  Cu,  .31  Fe, 

,08  Mn,  IS. 8  Cr,  2.98  Ti,  1.74  Al, 
17.8  Co,  Balance  Ni 

19  Cr,  11  Co,  10  Mo,  3  Ti,  1.5  Al, 
Balance  Ni 

.08  C,  .66  Si,  .05  Cu,  .3  Fe,,.06Mn, 
19.7  Cr,  2.4  Ti,  1.18  Al,  l.OlCo, 
Balance  Ni 

18.75  Cr,  .57  Fe,  4.25  Mo,  17.61.  Co, 
3.09  Ti,  3.15  Al,  Balance  Ni,  .1  Cu 

.61  Al,  1.05  Ti,  2.90  Mo,  U5  Cr, 

.1  Co,  17.2  Fe,  6.20  Ta,  2.75  W, 

•1  Si,  .21  Mn,  Balance  Ni 


Actual 

Actual 

Actual 

Actual 

Actual 

Actual 

Actual 

Actual 

Nominal 

Actual 

Nominal 

Actual 


69.  6  Mo  WASPALLOY 


.136  C,  .56  Mn,  .31  Si,  19.68  Cr, 
6.32  Mo,  13.5  Co, 2.59  Ti,  1.04  Al, 
2.5  Fe,  Balance  NI 


Actual 


VI 


Alloy 

.  TITANIUM  ALLOYS 

70.  T1-A55 

71.  T1-75A 

72.  T1-150A 

73.  Ti-0.2  02 

74.  Ti-0.2  C 

75.  Ti-N2 

76.  T1-.07  N2-.2  02-.2  C 


Percent  Composition 


.051  C,  .058  02,  .01  Fe,  Balance  Ti 
Not  Reported 

2.6  Cr,  1.3  Fe,  0.2  0^  Balance  Ti 

.25  02,  1.0  W,  .061  C,  Balance  Ti 

.2  C,  Balance  Ti 

.2  N2,  1.0  W,  Balance  Ti 

.26  C,  .08  N2,  .1  Fe,  .19  02, 
Balance  Ti 


77.  K-151A  CERMET 

78.  K-162B  CERMET 

79.  K-183A  CERMET 

80.  Ti-4  Al-3  Mo-1  V 

81.  Ti-4  Al-4  Mn 

82.  Ti-5  Al-2.5  Sn 

83.  Ti-5  Al-2.5  Sn-.07  N2 

84.  Ti-5  Al-2.5  Sn-.2  02 

85.  Ti-5  Al-2.5  Sn-.2  C 

86.  Ti-5  Al-2.5  Sn  COMPLEX 

87.  Ti~6  A1 

88.  Ti-6  Al-4  V 

89.  Ti-6  Al-4  V-.07  N2 


58  Ti,  19  Ni,  7.5  Nb,  .5  Ta,  15  C 

52  Ti,  25  Ni,  5  Mo,  4.5  Nb,  3  Ta, 
13.2  C 

41  Ti,  32  Ni,  3  Mo,  2.5  Al, 

7.5  Nb,  2.5  Cr,  .5  Ta,  12  C 

4  Al,  3  Mo,  Balance  Ti 

4  Al,  4  tta,  Balance  Ti 

5  Al,  2.5  Sn,  Balance  Ti 

.04  C,  .07  N. ,  .08  Fe,  .12  02 
5.04  Al,  2.5  sn,  Balance  TI 

•05  C,  ,007  N2,  .1  Fe ,  .2  02, 

.5  Al,  2.5  Sn,  Balance  Ti 

.2  C,  .05  02,  .007  N2,  ,1  Fe, 

.5  Al,  2.5  Sn,  Balance  Ti 

Complex  of  the  previous  three 

.06  C,  .035  N,,  .578  Al, 

Balance  Ti 

.02  C,  .18  Fe,  .011  N,,  6  Al 
4  V,  Balance  Ti 


Actual 

Nominal 

Actual 

Nominal 

Actual 

Actual 

Actual 

Actual 

Actual 

Nominal 

Nominal 

Nominal 

Nominal 

Nominal 

Nominal 

Nominal 

Actual 

Nominal 


Not  Reported 


Alloy 


Percent  Composition 


90.  Ti-6  Al-  4V-.2  02 


91.  Tl-6  Al-4  V-.2  C 


.07  C,  .005  N2,  .12-. 13  Fe,  67- 
82  H2  ppm,  201  02,  6.03-6.37  Al, 
4.18  V,  Balance  Ti 

.23-. 24  C,  .014-. 016  N2,  .10- 
.13  Fe,  75-92  H2  ppm,  .113-. 120  02, 
6.10-6.35  Al,  4.12-4.24  V,  Bal.  T1 


92. 

Ti-6  Al-4  V-.07  N2- 
.2  02-.2  C 

.19-. 22  C,  .052-. 053  N2,  .09- 
.10  Fe,  62-63  H2  ppm,  .194-. 201  02 
5. 74-6J8  Al,  4.12-4.30  V,  Bal.  Ti 

93. 

Ti-7  Al-4  Mo 

Not  Reported 

94. 

Tt-7  Al-3  Mo 

.020  C,  .03  Fe, 
.004  H2,  Balance 

7.0  Al,  2.7  Mb, 

Ti 

95. 

Ti-8  Al-1  Mo-1  V 

.023  C,  .09  Fe, 

.013  N2,  7.6  Al, 

1.0  V,  1.1  Mb,  .010-. 014  H2, 
Balance  Ti 


f 


96.  Ti-3  Mn-.2  02  .02-. 03  C,  .016-. 034  N2,  .80- 

.83  Fe,  36-39  H.  ppm,  .150-. 175  02, 
1.00-1.06  V,  3.0-3. 2  Mn,  .90- 
1,04  Cr,  1,15-1.20  Mo,  Balance  Tl 


97.  Ti-3  Mn-.07  N2  .03  C,  .052-. 053  N2,  .88-1.06  Pe, 

43  H2  ppm,  .122-1.44  02,  1.00-1.06 
2,50-2.90  Mn,  .96-1.08  Cr,  1.03  Mo, 
Balance  Tl 


98.  Ti-3  Mn-. 2  C 

99.  Ti-3  Mn  COMPLEX 

100.  Ti-3  Mn-.07  N2- 
.2  02-.2  C 

101.  Ti-4  Mn 

102.  Ti-8  Mn 


.21-. 22  C,  .014-. 017  N2,  1.08- 
1.12  Fe,  96-108  H2  ppm,  .103-. 148  02 
1.29-1.35  V,  3.00-3.05  Mn,  .94- 
.98  Cr,  1.35-1.36  Mb,  Balance  Tl 

.02-. 07  C,  .011-. 012  N2,  .86-.S8  ft, 
48-64  H2  ppm,  .091-. 082  02,  1.06- 
1.24  V,  3.08  m,  .80-. 96  Cr,  1.12- 
1.20  Mo,  Balance  Tl 

.11-. 15  C,  .056-. 057  N2,  .98-1.0  ft, 
54-58  H,  ppm,  .189-. 216  02,  1.00- 
1.18  V,  2.5-2. 6  Mn,  .94-. 96  Cr, 
l.r-*  1.38  Mo,  Balance  Ti 

.14  C,  .02  N2,  4.20  Mo,  Bal.  Tl 

8  Mn ,  92  Ti 


Actual 

Actual 

Actual 

Actual 

Actual 

Actual 

Actual 

i 

Actual 

Actual 

Actual 

Actual 

Nominal 


171 


Percent  Composition 


Alloy 

103. 

Ti-30  Mb 

104. 

Ti-13  V-ll  Cr-3  Al 

105. 

Ti-16  V-2.5  Sn 

VII. 

OTHER  NON-FERROUS  ALLOYS 

106. 

BERILCO  #25 

107. 

PURE  LEAD 

108. 

0.5  Ti-MOLYBDENUM 

109. 

PURE  TANTALUM 

no. 

10  W-TANTALUM 

in. 

PURE  TUNGSTEN 

.011-.C16  C,  .065-. 069  N2,  29.88- 
30.17  Mo,  Balance  Ti 

.029  C,  .14  Fe,  .028  N2,  3.0  Al, 
13.5  V,  11.2  Cr,  .010  H2,  Bal.  Ti 

16  V,  2.5  Sn,  Balance  Ti 

Not  Reported 
93.995  Pb 

0.46  Ti,  0.02  C,  Balance  Mo 

.08  Cb,.01  W,  .01  Fe,  Balance  Ta 

90  Ta,  10  W 
Not  Reported 


Actual 

Nominal 

Nomina.. 

Actual 

Actual 

Nominal 
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Sheet  -  .125  In.  thick  or  less.  Symmetrical  Edge  Notch,  Flank  Angle  0  Degree  Grain 
direction  parallel  to  lengthwise  axis 


41  'O 

<u  <u 
in  u  , 
<0 


w 


I  ■ 


H  J 


W 


•J*  ng 


301X1 

ana 


•mtaiM 

‘JLSia 


i»H30l 


gfc.d-il 

•b»ii 


.6885* 

6.652 

1.863 

8.423 

6.625 

3.374* 

7 . 863 
5.230 

1.719* 

9.579 

8.123 

3.512* 

8.680 

6.007 

1.329* 

10.44 

.9326 

9.012 

I. 863 

II. 20 
8.813 

3.433 

10.33 

6.828 

3.553 

11.25 

7.761 

1.618 

12.71 

BEfli 

n  O  sT 
kO  CM  O 

CO  •  vO 
•  <N  • 
r— (  *H  CTs 

39 

nno 

I. 893 
14.90 

II. 69 

H  r>»  r-« 

rH  m  CT> 

m  •  oo 

•  CO  • 

C0  H  00 

1.808 

16.31 

13.76 

3.628 

14.59 

9.994 

1.384 

13.37 

rIiIH 

3.539 

14.73 

9.630 

<r  p<*  co 

^  VO  00 

00  •  • 

• 

H  H  H 

rH  CO  ON 
vt  cr> 

•  • 

•  in  o 

n  H  H 

2.092 

16.43 

AXIAL  (Completely  Reversed) 
Effect  of  Frequency 

1.711 

16.53 

. :  n 

RSjjpSjjl 

1.936 

21.29 

17.84 

3.512 

18.92 

13.10 

2.401 

18.86 

1.875 

18.12 

3.463 

19.35 

12.75 

1.980 

23.07 

19.28 

2.547 

20.01 

2.318 

22.40 

1.862 

26.37 

20.74 

3.560 

23.42 

15.27 

r>-  vO  oo 

O  N  ^ 
p*  •  • 

•  ps.  cn 

H  (N  N 

3.640* 

24.55 

16.78 

2.924* 

22.97 

3.141* 

30.35 

1.781* 

35.07 

27.76 

3.549* 

30.77 

20.08 

1.984* 

36.25 

30.29 

3.591* 

31.83 

21.87 

3.561* 

27.98 

o  a 

BSi 

CEEjfl 

i9399 

029 

o  a 

xeauoN 

TTn<lT9ft 

Tinqi®rt 

nnqi*tt 

n»qi9ft 

-[BttUON 

E-dS 

t-dS 

E-as 

E“dS 

E-dS 

E-dS 

rz 

9*3 

rz 

rz 

Z'Z 

VZ 

08 

08 

08 

08 

08 

08 

9 

(8 

9 

OE-VZ 

9 

0E-*Z 

p»qtn°d 

04139T3 

Q 

m 

E-l 

E-l 

E-J. 

E-l 

SP-5  Sheet-. 125  in.  thick  or  less.  Hole  Notch,  grain  direction  parallel  to  lengthwise  axis 
SP-6  Sheet-. 125  in.  thick  or  less.  Unnotched,  Constant  stress  test  length 


h|  *uTi 

g|  ’Jins 


■K 

t"  h  cm 
so  in  o 
i —  r~-  -<j- 


i— I  so  i'* 

St  (O  rl 
00  •  iA 

•  O  • 
H  .-I  00 


H  H  vO 

r-»  <o  m 
oo  •  cm 


*  * 

r-i  o  m  os 
tn  os  os  oo  m 

1/1  H  O  00  sO 

•  •  •  so  » 

rl  MO  •  sO 


O  sj  Os 

so  r*  co 
sO  so  H 


* 

cn  -o  «n 
>-t  so  ro 

N  H  N 


oo 

m 

m 

os 

CM 

m 

ao 

• 

• 

• 

O 

sO 

cm 

H 

CM 

ro 

lA 

cn 

cn 

CA 

• 

• 

• 

CM 

SO 

CM 

« 

vj-  r-~  oo  os  O'  O 

CO  O'  vT  >A  -O  © 

rs  <  «  'O  •  • 

■  sj  o  •  r"i  o 

f-4  fs|  fsj  H  N  N 


*  * 

N  IN  I/I  sO  H  S' 

sO  -4  Cl  H  sO  ifl 

A*  *  •  iO  *  1 

i  A  IN  '  H  vO 

H  SI  (N  M  n  (N 


•O  » 


TinqTart  TTnqi»n  T»««H 


o 

CM 

r* 

sC 

CM 

(A 

• 

• 

• 

cn 

in 

CM 

* 

Os  r*  00 
stNO 
•A  •  * 

•  Q  O 

MAN 


* 

o^oo 

Sf  tA  I" 

vO  •  • 

•  <r  so 
A  N  H 


« 

i-«  <*S  A* 

os  ao  ao 

SA  •  • 


* 

N  N  'I 
OS  Os  IA 
sO  •  * 

•  ^  o* 

S3  A  H 


lnmaiMj] 


C-dS  C-dS  €-dS 


9  *Z 

rfc 

Li 

’Z 

09 

os 

09 

p*i4JT!Iodoi50»T3 


Sheet  -  .125  In.  thick  or  less.  Fillet  notch,  grain  direction  parallel  to  lengthwise  axis 
Plate  -  over  .125  in.  thick.  Unnotched 


<  w 

N 

g  d 

a  e 

£  as 

esfi  m 


M  W 

s  J 


* 

* 

r-  cn  | 

cn  cn  m 

on  oo 

sO  CM  CM 

cm  m 

co  -a*  vO 

•  • 

•  •  • 

i—l  00 

H  00  vO 

in 

cn  o  cn 

OO  H 

\D  CN  H 

vC  • 

00  •  00 

•  H 

•  rl  • 

pH  H 

H  H  00 

m  r*» 

m  o  sr 

CM  O 

vON  O 

00  * 

00  •  NO 

•  CN 

•  CN  • 

rl  rl 

rH  rl  ON 

cn  o 

cn  o  ON 

on  m 

O'  O'  VO 

H  • 

00  •  • 

•  ^ 

•  ^  rl 

CM  H 

rl  r|  rl 

cn  o 

V0  <f  CM 

h*  r^. 

HNS 

cn  • 

ON  •  • 

•  m 

•  >0  CN 

CM  <-4 

rH  rl  rl 

.-4  no 

r>*  cn  m 

in  oo 

S  «  4 

00  • 

On  •  • 

•  oo 

•  on  in 

CM  t-4 

*— 1  r|  rl 

NO  -4 

in  on  o 

oo  -4 

in  ri 

O  • 

r^.  •  # 

•  o 

•HS 

cn  cm 

-4  CM  r4 

* 

0\  ON  cn 
H  N 
MAH 
•  #  • 
H  ON  00 


On  ON  On 
iH  <f.  m 
r*.  *'  • 

•  CM  o 


^ 

cn  m  <r 

r*>.  •  • 

•  co  *H 

H  rl  H 


00  H  NO 
O  fA 
00  •  • 
•  V0  CO 


Sf  N  00 

vO  oo 
00  •  • 
•  <r 


vO  On  ^ 
cn  cm  oo 
Cn  •  • 
•  »-i  r> 
HMH 


O  N00 
00  O  N 
On  •  • 

•  CO  ON 
H  CM  H 


HNnO 
00ON 
N  •  • 

•  w*  f* 
HO  N 


«*■  cn  Hrt 
©  a 


« 

Sin  O' 
CM  CM 
On  • 
•  'O  O 


* 

CM  Hf  00 
CM  CM  00 

rs  cm  cn 


* 

CM  0^  CO 

m  o  vo 

00  O  VO 


* 

oo  o  m 
o  o  m 
ov  h  in 


* 

m  in  m 

00  CM  CM 
N  •  • 

•  CM  O 


1*5 

cm 

NN4 

rs.  NO  00 

« 

«»  r- 

O  cn 

£  ■£  5 

Q  W 

O  *n 

vo  n  n 

NO  rv  sf 

CM  O' 

5  >4 

O  r-»  so 

o 

GO  •  • 

rv  ♦  • 

ON  • 

in  • 

On  ♦  • 

S/  %4 

•  ^ 

•  NO  © 

•  r*^  cn 

•  CM 

St  NO 

•  oo  cn 

o 

cn  cm 

-4  CM  CM 

H  CM  CM 

CM  CN 

•  cn 

rH  CM  CM 

* 

cn  pn.  m 
CM  ON  #«4 
•  • 
•  rt  n 
h  cn 


TWU0N  Tinq^tfi  Iinq|»ft|lwuoN  ‘A'd'S  ITnqi*M 


C-dS  C-dS  C-dS 


9’3  9*3 


S3  IR 


/r  |  0093  I  OC-93 


p»qt|Xodoai3®ia 


2024  ALUMINUM  (Cont'dN 

(SHEET  OR  PLATES)  Composition:  See  Page  166 

FATIGUE  STRENGTH  DISTRIBUTIONS  AND  THEIR  PARAMETERS 
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HT-53:  1450°F  8  hrs.,  F.C.,  L850°F  5  min.  A.C.,  1375°F,  15  rain.  A.C 

HT-54 :  1250°F  1  hr.,  F.C.  to  900°F,  A.C. 
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HT-60:  Annealed  450#F,  0.5  hrs.,  Shot  peened  with  .010  Steel  Shot  to  Almen  #10 
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HT-63:  Solution  H.T 
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A-2.1  A  METHOD  FOR  INTERPOLATING  INTERFERENCE  VALUES 


The  values  of  Interference  given  in  Appendix  2  are  highly  non-linear. 
As  a  result,  when  interpolation  is  required  to  obtain  an  interference  value,  a 
higher  order  interpolating  procedure  Will  be  more  accurate  than  a  linear  one. 

One  method  for  interpolating  at  a  higher  order  is  to  construct  what 
is  called  an  "interpolating  polynomial."  The  interpolating  polynomial  gives  an 
nth  order  interpolation  using  n  +  1  interference  values  from  the  table. 

A-2.1.1  Higher  Order  Interpolation  Procedure 


There  are  several  types  of  interpolating  polynomials,  one  of  the  easi¬ 
est  to  use  is  constructed  using  Neville's  Algorithm.  This  method  and  the 
symbols  used  are: 


X 


is  the  value  of  an  interference  parameter  which  does 
not  appear  in  the  interference  tables. 


X0 ,Xi , . . . ,Xn  are  any  n+1  consecutive  values  of  the  interference 

parameters  which  do  appear  in  the  interference  tables. 

d=0,l,...,n  are  the  consecutive  indices. 


Pjt  0(X)  K“0,l,...,n  are  the  n+1  interference  values  correspond¬ 

ing  to  the  values  of  the  interference  parameter  X0,Xj_, 

' • • *^n* 

pn,nOO  is  the  desired  n**1  order  interpolated  interference  value 

for  the  interference  parameter,  X. 


(X) 


pn,nOO  is  che  last  entry  of  the  triangular  array  of  polynomials 
(as  shown  in  Table  A-2.1  where  Pk  d(X)  is  defined  as: 


pk,d+lOO 


(xk-x)  Pfe-i^CX)  -  (Xk_d-1-X ) Pjj t d (X) 

*k  -  xk-d-l 


where  d  -  0,  1,  ...,  n 

k  ■  d+1,  d+2,  . . . ,  n  , 

It  is  important  to  note  that  the  value  of  the  Interference  parameter, 
X,  for  which  the  interference  must  be  interpolated  can  lie  anywhere  in  the 
interval  X0  «  »  <  Xn.  However,  the  value  of  X  may  also  lie  outside  this  inter¬ 
val,  in  which  case  the  problem  becomes  one  of  extrapolation  rather  than  inter¬ 
polation. 


When  the  interpolation  is  required  for  only  one  interference  para¬ 
meter,  say  X,  it  is  called  one-way  interpolation;  and  when  the  interpolation 
is  done  for,  say  two  or  three  parameters,  X  and  Y  or  X,Y,  and  Z,  it  is  called 
two-way  and  three-way  interpolation  respectively. 


Interference 
Values  From 


Table  A-2.1  Construction  of  the  Triangular  Array  of  Polynomials  for  nfch  Order  1 


Interference 

Parameter 

Interference 
Values  From 
Table 

Interference 
When  U-.025 
and  y»-2.5 

-y-x* 

m 

Hi 

1 

u 

1 

m&m 

X  -1.0 

0 

X.-2.0 
X“2,5)  X2-3.0 
X3-4.0 

po,o--3079 

P,  „-.0486 

pi,i-°559 

P,  ,-.0876 

?3  j».0638 

P2,2“-0747 

?3’2-.0817 

P3  3-.0782 

Table  A-2.2  Polynomials  for  Third  Order  Interpolation 


I 


The  following  examples  will  illustrate  the  application  of  this  higher 
order  interpolating  method  for  one,  two,  and  three-way  interpolation  problems. 
In  the  example  chosen,  linear,  third  order  and  fifth  order  interpolations  were 
used  but  this  method  applies  equally  well  to  even  and  odd  order  interpolations. 

A-2.1.2  Example  Problem  for  One-Way  Interpolation 


A  particular  part  has  Smallest  Extreme  Value  distributed  strength 
with  parameters  ($,M)  and  is  subjected  to  Normally  distributed  stresses  with 
parameters  (p,cr).  The  calculated  interference  parameters  are: 

a  ®  3a  -  .025 
y  =  0(p-M)  -  -2.50 

What  is  the  %  Interference? 


From  the  tables  in  Appendix  2  on  Page  422,  it  is  noted  that  a  ■  .025 
appears  in  the  table  but  y  ®  -2.5  does  not;  therefore,  the  value  of  interfer¬ 
ence  must  be  interpolated.  The  values  of  interference  as  read  from  the  above 
tables  are: 

Interference 
-y  When  a® •  025 


2.0  .1266 

3.0  .0486 

Linear  Interpolation: 

If  the  quick  and  easy  linear  interpolation  is  used,  the  interference 
corresponding  to  a®. 025  and  y  ®  -2.5  is  found  as: 

Interference  ■  .1266  +  (3. 0-2. 5) 

-  .1266  -  .0390 

-  .0876 

or  Percent  Interference  ■  8.76% 


Third  Order  Interpolation: 


Using  a  third  order  interpolation,  by  constructing  the  appropriate 
polynomials,  Pk  d(X),  the  value  of  P3  3^)  is  the  interpolated  value  of  inter¬ 
ference.  Where* P3  3(X)  is  calculated’as  follows  and  the  results  are  tabulated 
in  Table  A-2.2: 


P  ,(X) 
MH 


cya^.dW  ■  (xk-d-rx)pk.d(x) 

*k  ‘  \-d-i 


pm<2,5)  " 


2-1 


.0359 


555 


p  (2,5)  .  (3-2.5K.1266)  -^(2-2.5)  t.0486)  .  ,0876 

P  (2.5,  ,  (4-2,5) (.0486)  -  (3-2. 5) (.0182),  _  ^ 

3 )  X  A  “  3 

P  /,  c.  (3-2.5) (.0359)  -  (1-2.5)  (.0876)  = 
r2,2'‘ '  '  3-1 

p  (2.5)  „  (4-2.5)  (.0876)  -  (2-2.5)  (.06381  =  >081? 

3,2  *  4-2 

_  (4-2.5) (.0747)  -  (1-2.5) (.0817) 

Interference  =  P3  3(2. 5)  =  - ^ - 4  - 


=  .0782 


or  Percent  Interference  =  7.82% 

Thus,  it  c.an  be  seen  that  the  third  order  interpolation  gives  signifi¬ 
cantly  more  accurate  results  than  the  linear  interpolation. 


Fifth  Order  Interpolation; 


Using  a  fifth  order  interpolation,  by  constructing  the  appropriate 
polynomials  (P^  ,(X)),  the  value  of  P^  _(X)  is  the  interpolated  value  of  inter¬ 
ference  where  *aP_  -(X)  is  calculated3 as  follows  and  the  results  are  tabulated 
in  Table  A-2.3.  3,3 


pi,i<2-5> 

P2,i(2.5) 

P3,l(2.5) 

P,,l(2.5) 

P5,l(2,5) 


(1-2.5) (.6321)  -  (0-2.5) (.3079)  17R.n 

1-0  ~ 

(2-2.5) (.3079)  -  (1-2.5) (.1266)  _  0i03595 


(3-2.5) (.1266)  -  (2-2,5) (.0486) 
3  -  2 

(4-2,5) (.0486)  -  (3-2.5)  (.0182) 

4- 3 

(5-2.5) (.0182)  -  (4-2.5) (.0067) 

5- 4 


0.08760 

0.06380 

0.03545 


P2,2(2.5) 

P3j2(2*5) 

P4>2(2.5) 

P5>2(2.5) 


(2-2.5)  (-.17840)  -  (0-2.5)  (.03595)  _  i0q9^ 

(3-2.5) (.03595)  -  (1-2.5) (.08760)  .  <0746875 

(4-2.5) (.08760)  -  (2-2.5) (.06380)  „  .0816500 
4-2 

(5-2. 5)(.06380)i, -13,^,5).C.03545),  .  j070g875 
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p  (2.5,  ,  Si-2, jj, (.0746875)  -  (l-:2. 5) (.0816500)  , 

4,3  4-1 

P5>3(2.5)  -  ^-liX-WSOO^-  (2-2-5)  (-070887^  .  ,07985625 

P4  4(2.5)  .  l^^C-OT.^JL.-....(P-2.,3M.0781687a  .  .0777914 
P5j4(2.5)  -  i-^a^7816875).-_11.215)i.0798,5625)  .  ,07880156 

Interference  -  P,  ,(2.5)  ■  i5--?.-..5.)-!: ^iji  .-.(^2 ■  5 >( ■  0788M5 62.. 

or  Percent  Interference  »  7-83# 

Thus,  the  fifth  order  interpolation  does  not  significantly  improve  on 
the  third  order  result  for  this  example. 


A- 2 .1.3  Example  Problem  for  Two-Way  and  Three-Way  Interpolation 
Two-Way  Interpolation : 

The  preceeding  example  demonstrates  one-way  higher  order  interpolation 
The  method  can  be  easily  extended  to  the  problem  of  two-way  and  three-way  inter¬ 
polation  where  values  of  2  or  3  different  interference  parameters  are  not  in 
the  tables.  [For  example,  a  and  y  or  C,A,  and  B(x)J 

For  two-way  interpolation  of  order  n,  the  value  of  interference  cor¬ 
responding  to  interference  parameters  X  ■  X*  and  Y  -  Y*  (as  shown  in  Table  A- 
2.4)  is  determined  as  follows: 

(1)  Determine  the  interference  values  ...,  In*  for  X  -  X0, 

...  Xn,  at  Y  ■  Y*  by  performing  the  n+1  one-way  interpolations  as  discussed  in 
Section  A-2.1.2. 

(2)  Determine  the  value  of  interference,  I*,  at  X  ■  X*,  Y  ■  Y*  by 
performing  the  one-way  interpolation  on  the  interference  values  Ig*,I]*,  .... 
In*  obtained  in  (1). 

The  next  section  on  three-way  Interpolation  contains  several  examples 
of  two-way  interpolation. 

Three-Wav  Interpolation: 

The  need  for  a  three-way  interpolation  may  arise  when  the  fatigue 
strength  is  Weibull  distributed  and  the  stress  is  Normal  distributed.  If  the 


\.>l 


Interference  Parameter 


interference  parameters  have  values,  say  B(x)  ■  b*,  A  ■  A*  and  C  ■  C*,  and 
none  of  these  values  are  available  in  the  tables,  then  a  three-way  interpola¬ 
tion  is  required  to  obtain  the  interference  value. 

(1)  Determine  the  interference  value  I  **  for  A  -  A*,  C  =  C*,  and  for 

a  given  value  of  B(x)  *  bfi  appearing  in  the  tahles  by  performing  a  two-way 
interpolation  as  discussed  above.  By  repeating  this  for  various  values  of 
B(x)  ■  bf.,b1,b_,  ...  b  ,  appearing  in  the  tables,  n+1  interference  values  I_** 
1^** .  1^,  will  be  generated  as  shown  in  Table  A-2.5  (a)  and  (b). 

(2)  Determine  the  value  of  interference  I**  at  B(x)  *  b*,  C  «  C*, 

A  ■  A*  by  performing  the  one-way  interpolation  on  the  interference  values 
IQ**,  •••>  obtained  in  (1)  as  shown  in  Table  A-2.6. 

As  an  example  of  three-way  interpolation,  consider  the  Zxample 
Problem  No.  1,  given  in  Section  9,  page  140,  where  the  interference 
parameters  are  b  «•  1.808,  C  **  3.72,  A  ■  .0875.  None  of  these  values  ap¬ 
pear  exactly  in  the  tables,  so  a  three-way  interpolation  is  required. 

As  it  involves  a  great  number  of  calculations  for  high  order 
three-way  interpolations,  a  second  order  interpolation  will  be  performed 
to  demonstrate  the  technique. 

For  the  three-way  second  order  interpolation  the  following  in¬ 
terference  values  were  selected  from  the  Tables  in  Appendix  2,  pages  388 
through  392 . 

The  first  step  is  to  perform  a  total  of  three  "two-way"  inter¬ 
polations  where  each  two-way  interpolation  requires  four  "one-way"  in¬ 
terpolations.  The  necessary  calculations  are  as  follows. 

(1)  Two-way  interpolation  within  Table  A-2.7(a)  for  b  ■  1.0, 

A-.0875,  03. 72 

(i)  One-way  interpolation  for  A  •  -.2,  b  ■  1.0,  C  -  3.72 


xk 

-  ck 

(interference 
values  from 
tables) 

pk,0 

Pktl  Pk,2 

Xo 

-  3.00 

.1372 

X1 

-  3.50 

.1210 

.113872 

(X  -  3.72)  v 
*2 

-  4.00 

.1081 

.115324  .11491744 

B(x)  -  be . -V*  *  I  *  B(x) 


Table  A-2.6  Determination  of  the  Interference  Value 
Corresponding  to  A  *  A*,  C  =  C*,  and 
B(x)  -  b*  by  One-Way  Interpolation. 
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-0.2 

.0359 

.0166 

.0088 

b  -  3.0 

0.0 

.0254 

.0116 

.0062 

0.2 

.0176 

.0080 

.0042 

Table  A-2.7  Values  of  Intarferenca 
a a  Read  from  the  Tables  In  Appendix  2 
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.113872 


Pk.d+100 

Pn(3.72) 

P21(3.72) 


(Xk  -  X)  Pk-l,d(X)  “  <xk-d-l-X)Pk>d(X) 

Xk  -  xk_d_Y  " 

(3.5  -  3.72)  (.1372)  -  (3.0  -  3.72)(.1210) 
'3.5  -3.0 

(4.0  -  3.72) (.1210)  -  (5*5-3.72) (.1081) 

4.0  -  3.5 


=  .115324 


p22(3.72) 


(4. 0-3. 72)(. 113872)  -  (3.0  -  3. 72) (. 115324) 

_____ 


.11491744 


In  a  similar  manner  the  following  results  are  obtained! 


(ii)  One-way 

interpolation  for  A 

*  0,  b  «  1.0 

C  =  3.72 

Ck 

pk,0 

Pk,l 

PM 

3.00 

.1095 

3.50 

.0963 

.090492 

(3.72) 

4.00 

.0860 

.091688 

.09135312 

(iii)  One-way 

interpolation  for  A 

-  .2,  b  -  1.0,  C  -  3.72. 

Ck 

Pk,0 

Pk,l  Pk.2 

3.00 

.0850 

3.50 

.0749 

.070456 

(3.72) 

4.00 

,0668 

.071336  .07108960 

364 


(iv)  One-way  interpolation  on  the  values  obtained  in  (i) , 


(ii) ,  and 
b  =  1.0, 

(iii),  to  obtain  interference  for  A  =  .0875 
C  =  3.72 

values  obtained 
in  (i),  (ii), 
and  (iii) 

Ak 

Pk,0  Pk,l 

Pk,2 

-0.2 

.11491744 

0.0 

(.0875) 

0.2 

.09135312  .08104375 

.07108960  .08248783  .082081676 

Two-way  interpolation  within  Table  A—  2.7; 
b  =  2.0,  C  =  3.72 

A  =  .0875, 

(i)  One-way  interpolation  for  A  *  -.2,  b  * 

■  2.0,  C  -  3.72 

Ck 

Pk,0  Pk.l 

pk,2 

3.0 

.0639 

3.5 

(3.72) 

4.0 

.0490  .042444 

.0386  .042444 

.042444 

(ii)  One-way  interpolation  for  A  ■  0,  b  - 

2.0,  C  -  3.72 

Ck 

Pk,0  pk,l 

Pk,2 

3.0 

.0477 

3.5 

.0364  .031428 

(3.72) 

4.0 

.0286  .032968 

.03253680 

(iii)  One-way 

interpolation  for  A  ■  ,2,  b 

-  2.0,  C  -  3.72 

Ck 

pk,o  \.l 

pk,2 

3.0 

.0348 

3.5 

(3.72) 

4.0 

.0265  .022848 

.0207  .023948 

.023640 
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(iv)  One-way  interpolation  on  the  values  obtained  in  (i) , 
(ii) ,  and  (iii) ,  to  obtain  the  interference  for 
A  -  .0875,  b  =  2.0,  C  =  3.72. 


Ak 

Pk.O 

Pk,l 

V  . 

-0.2 

.042444 

0.0 

.0325368 

.02820240 

(.0875) 

0.2 

.023640 

.02864445 

.028520123 

Two-way  interpolation  within 

Table  A-2.7(c)  for 

A  *  .0875,  b 

=  3.0,  C  =  3.72. 

(i)  One-way 

interpolation  for  A  =  -.2,  b 

=  3.0,  C  =  3.72 

Ck 

Pk,0 

Pk,l 

Pk,2 

3.0 

.0359 

(3.72) 

4.0 

.0166 

.022004 

5.0 

.0088 

.018784 

.02084480 

(ii)  One-way  interpolation  for  A  ■  0,  b 

-  3.0,  C  =  3.72 

Ck 

Pk.O 

Pk,l 

Pk,2 

3.0 

.0254 

(3.72) 

4.0 

.0116 

.015464 

5.0 

.0062 

.013112 

.01461728 

(iii)  One-way  interpolation  for  A  ■  .2, 

b-  3.0,  C  -  3.72 

Ck 

pk,0 

pk,l 

pk,2 

3.0 

.0176 

(3.72 

4.0 

.0080 

.010688 

5.0 

.0042 

.009064 

.01010336 

(iv)  One-way  interpolation  on  the  values  obtained  (i) , 
(ii) ,  and  (iii)  to  obtain  the  interference  for 
A  -  .0875,  b  -  3.0,  C  -  3.72. 


-0.2  .0208448 


(.0875) 


0.0 

0.2 


.01461728 

.01010336 


.01189274 

.01264244  .012431587 


In  the  second  step,  the  final  result  is  obtained  by  performing 
the  one-way  interpolation  for  b  =  1.808  on  the  interference  values  cal¬ 
culated  in  part  (iv)  of  (1),  (2),  and  (3)  above. 


bk 

(1.808)  1,0 

2.0 


.082081676 

.028520123 


.038803941 


3.0  .0- 2431587  .031609122  .035897  -  .0359 


Thus  the  interference  for  b  =  1.808,  A  »  .0875,  C  -  3.72  as  ob¬ 
tained  by  a  three-way  second  order  interpolation  is: 

Interference  *  .0359  or, 

Percent  Interference  ■  3.59% 

Interference  obtained  by  the  three-way  linear  interpolation 
(as  shown  in  Example  Problem  No.  1  in  Section  9,  page  141)  is: 

Percent  Interference  ■  3.98% 

Hence,  a  substantial  improvement  is  provided  by  using  the  higher  order 
interpolation. 
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A-2.2  STRESS  DISTRIBUTION  -  NORMAL 

STRENGTH  DISTRIBUTION  -  WEIBULL 
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TABLES  OF  INTERFERENCE  [F(X)] 


Stress  Distribution:  Normal  Strength  Distribution:  Weibull 

y  ■  SeqU  (Equivalent  X  Depend  on  the 

Stress)  material  and  the 

o  *•  0  b  operating  conditions 


6 


Figure  A-2.1  Interference  with  Standard  Deviation  of  Stress 
Equal  to  Zero 


Weibull  density  function  is:  j 

f(x)  -  -  (jL.7  *g$ml  e^0_X°) 

0  -  x0  x0  ; 

Cumulative  distribution  function,  F(x) 


f(x)  dx 


ut  7  ‘(^t)  .  dy  ’  (-H£-) 

,(H^1 

F(x)  ■  Jf(x)dx  -  )  e“y  dy  ■  -e”y  J 


-  1  -  • 


where  X  •(jLzJii 

ve  -  x, 
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INTERFERENCE,  F(X) 
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STRESS  DISTRIBUTION  -  NORMAL  (<J=0) 


TABLES  OF  INTERFERENCE 


Stress  Distribution:  Normal 


Strength  Distribution:  Weibull 


V  m  Sequ 
a  0 


(Equivalent 

Stress) 


X 

o 

b  or  B(x) 


Depend  on  the 
material  and  the 
operating  conditions 


6 


f(*) 


Figure  A-22  Interference  with  Standard  Deviation  of  Stress 
Not  Equal  to  Zero 
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Figure  A-2, 3  Interference  With  Standard  Deviation  of  Stress  Equal  to  Zero 
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TABLES  OF  INTERFERENCE  [F(X) ] 


Stress  Distribution:  Normal 


U  -  S 


equ 


(Equivalent 

Stress) 


a  ■  0 


Strength  Distribution:  Smallest  Extreme 

Value 

M  Depend  on  the 

material  and  the 
6  operating  conditions 


f(x) 


Figure  A-2.7  Interference  with  Standard  Deviation  of  Stress  Equal 
to  Zero 


Density  Function  of  Smallest  Extreme  Value  is: 
f (x)  -  e  e 

Cumulative  Distribution  Function,  F(x)  ■  ^  f(x)  dx 
Lst  y  -  -e*8^  dy.-S.*6^  dx 


«(S.,U-M)  «S(S.qu-tf) 

** *  ’  ^  J  '  1  *  * 

-e’X 

F(X)  -  1  -  • 
where  X  ■  -8(S*qU-M) 
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..A-'2,5.2  Stress  Standard  Deviation  j*  0 

fi 


A20 


TABLES  OF  INTERFERENCE 


Stress  Distribution:  Normal  Strength  Distribution:  Smallest  Extreme 

Value 

Vi  =  S  (Equivalent  M  Depend  on  the 

ec*u  Stress)  material  and  the 

of  0  6  operating  conditions 


Figure  A-2.8  Interference  with  Standard  Deviation  of  Stress  Not  Equal 
to  Zero 
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:th  distribution  -  SMALLEST  extreme  value 


APPENDIX  3  EVALUATION  OF  THE  INTEGRALS  BY  NUMERICAL  ANALYSIS 


A-3.1  THE  PROBLEM 


In  this  report  stress/strength  interference  theory  was  used  to  pre¬ 
dict  failure  probabilities  for  several  important  cases  not  included  in  [1]. 
The  two  cases  of  interest  here  are  (1)  the  smallest  extreme  value  distributed 
strength  and  normally  distributed  stress  and  (2)  the  largest  extreme  value 
distributed  strength  and  normally  distributed  stress.  The  integral 


Pr [Failure] 


Pr(X-Y)  - 


F(y)g(y)dy 


(1) 


was  evaluated  where  F(*)  is  the  distribution  function  of  X  (the  random  varia¬ 
ble  associated  with  the  strength)  and  g(*)  is  the  density  function  of  Y  (the 
random  variable  associated  with  the  stress).  That  is,  we  suppose  that  Y  has 
a  probability  density  function  : 


.  i^2 


*Y(y)  mJIn  e 


dy ,  -oo  5  y  1  +00 


and  X  has  a  probability  distribution  function 

B(x-M) 


Fx(x)  -  1-e 


—os  —  x  —  +°° 


if  it  is  distributed  as  the  smallest  extreme  value.  If  X  is  distributed  as 
the  largest  extreme  value,  its  probability  distribution  function  is: 


Fx(x)  -  e' 


-B(x-M) 


i  x  * 


One  sees  that  if  z.  ■  1-Fx(x)  is  plotted  on  a  lnln  scale  it  will  be 
linearly  related  to  x  on  an  arithmetic  scale  as 

z  ■  P(x-M)  or  z  ■  -B(x"M) 

Hence  one  can  call  M  the  mode  of  the  distribution  and  S  (or  -0)  the  "slope"  or 
"shape"  parameter. 

Inserting  (2)  and  (3)  or  (A)  into  (1)  will  give  (5)  or  (6). 


Pr  (X*Y) 


(5) 


J 


Pr(X*Y) 


•![■ 


_e-6(y-M)-|  j 


A-3.2  SIMPLIFICATION  OF  THE  EQUATIONS 

Equations  (5)  and  (6)  are  quite  difficult  to  work  with  numerically 
since  each  one  involves  four  parameters.  It  is  desirable  to  define  a  new 
variable  of  integration  and  new  parameters.  To  this  end  we  let: 

u  “  ,  a  m  Bo  ,  and  y  ■  S(p-M) 


Then  (5)  of  Section  A-3.1  becomes 


Pr(X*Y) 


‘&T 


*  >  djr 


^_eS(ua+p-M)^- 

a  » 


and  this  further  reduces  to  ^ 

.  ou4y  -  ru* 

Pr(X^Y)  -  1-  -jgv  J  e  6  e  du  (7) 

Also,  (6)  of  Section  A-3.1  becomes 
1  1*  -(cu+y)  -  |u2 

Pr(X*Y)  -  I  e  e  du  (8) 

Inspection  of  (7)  and  (8)  will  show  that  the  integral  to  be  calcu¬ 
lated  is  of  the  same  form  in  each  case  except  that  the  uppermost  exponent  is 
negative  in  the  Largest  Extreme  Value  case.  Therefore,  to  further  Isolate  the 
numerical  analysis  problem  involved,  define  new  constants  A  and  G  such  that 
for  (7),  the  Normal-Smallest  Extreme  Value  Integral, 

A  •  a  and  G  ■  y 

but  for  (8),  the  Normal-Largest  Extreme  Value  integral, 

A  »  -a  and  G  •  -y 

Then,  the  calculation  of  the  probabilities  given  by  Equations  (5)  and 
(6)  of  Section  A-3.1  reduces  essentially  to  the  integration  of 

i  "  T*2  Au+G 

X  *  •'  du  <9) 
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Since  this  integral  does  not  exist  in  closed  form,  numerical  techniques  must 
be  employed.  The  ranges  of  a  and  y  of  practical  interest  and  the  associated 
values  of  A  and  6  are 

A  ■  where  oranges:  .001,  .005,  .010,  .025,  .05,  .075,  .1,  .2, .. .  3., 
3.5,...  9.5,10 

G  ■  Ay  where  y  ranges:  0. ,-l .  -40 

Note  that  a  itself  is  always  positive  and  y  is  always  negative  (or 
zero).  A  negative  a  would  imply  a  negative  standard  deviation.  This  situa¬ 
tion  is  without  mathematical  meaning.  If  y  were  positive,  the  material  in 
question  would  have  a  mean  strength  less  than  the  mode  of  the  stress.  Such  a 
situation  is  of  dubious  engineering  importance.  (6  is  always  positive.) 

Note,  also,  that  for  the  Normal-Smallest  Extreme  Value  case,  the  pro¬ 
bability  of  failure  is  one  minus  the  number  given  by  (9). 

A-3.3  NUMERICAL  METHODS  CONSIDERED  FOR  THE  PROBLEM 

The  integral  (9)  of  the  previous  section  could  be  evaluated  in  sever¬ 
al  ways.  Its  form  Immediately  suggests  Gauss ian-Hermlte  integration.  This 
method  can  be  summarized  as  follows: 

-  Au-HS 

Let  F(u)  ■  e  and  define  v  ■  u//2  , 

-4s  2 

F(u)du  -  j"  e"^  F(v)dv 
and  by  the  Gauss ian-Hermite  formulation* 

7?  J  *"y2  ™iv '  7*j,  V'V 

where  the  arguments  are  the  zeros  of  the  nth  Hermlte  polynomial 

Vv>  -  (-1)“  e*2  <e“v  ) 

and  the  coefficients  A^  are 


The  numbers  v.  and  A.  are  available  in  tables.  Ths  approximation  error  using 
this  method  depends  on  the  higher  derivatives  of  F.  For  the  n-point  case. 
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,  ~eo  <  £  <  +oo 


error 


niy?  F^2n^  (£) 

2n(2n) J 


(The  formula  for  computing  the  higher  derivatives  of  F  is  given  in  Table  A-3.2) 

Errors  in  Integrating  (9)  of  Section  A-3.2  were  computed  for  the  n  « 
2,  4,  and  8  point  formulas  for  the  entire  range  of  the  parameters.  It  was 
found  that  the  n  ■  8  formulas  would  give  acceptable  errors  for  small jA|  and 
large  Id.  Theoretically  at  least,  this  scheme  could  be  used  to  integrate  (9) 
for  any  A  and  G  if  the  number  of  points  taken  were  sufficiently  large. 

Despite  its  advantages,  the  Gausslan-Hermite  Integration  method  was 
not  used  to  prepare  the  tables  of  this  report  for  the  following  reasons: 

1)  The  higher  derivatives  of  F  become  quite  difficult  to  compute 
and  hence  the  errors  Involved  are  difficult  to  evaluate. 

il)  It  proved  useful  for  most  values  of  the  parameters  for  which 
the  approach  presented  below  proved  even  simpler. 


ill)  It  was  decided  that  using  only  one  method  (that  below)  was  more 
efficient  than  programming  different  methods  for  different 
values  of  A  and  G. 

If  users  of  the  tables  wish  to  compute  probabilities  of  failure  for  a's  and 
y's  not  included  ir.  this  report,  the  Gausslan-Hermite  method  may  prove  accept¬ 
able,  but  the  method  employed  here  is  recommended. 

The  tables  of  probabilities  Included  in  this  report  were  generated  by 
a  FORTRAN  program  run  on  the  University  of  Michigan's  IBM  360/67.  For  every  a 
and  y»  the  program: 

1)  Determined  limits  of  integration  following  the  logic  of  \ 

Section  A-3.,5,  ' 

11)  Used  Simpson's  rule  to  Integrate  (9)  within  these  limits,  and 


ill)  Added  to  the  above  value  the  integral  of  the  normal  curve 
neglected  by  using  the  calculated  limits. 

This  approach  is  Illustrated  in  Figure  A-3.1.  Its  major  advantages,  which 
will  be  made  clearer  in  later  sections,  can  be  listed  as 

1)  Automatically  preventing  computing  difficulties  by  limiting 
the  exponentials  to  be  computed  to  acceptable  values, 

11)  Having  an  easily  computed  error, 

ill)  Reducing  computing  time  by  using  the  error  function  subroutine 
wherever  possible  and  by  automatically  recognising  where  the 
integral  is  approximately  aero,  and 
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Au+G 

DERIVATIVES  OF  e~* 


F 


Au+G 

e 


Au+G 

■e 


Fd) 

F<2) 

.(3) 


AXF  [P1] 

A2F  [P1  +•  P2] 

3  1  2  3 

A  F  [P~  +-  3  P  +-  PJ] 

A4F  [P-1  +  7  P2  +-  6  P3  +■  P4] 


AnF  (Cj^P1  +•  C2nP2  +■  ....  CknPk+- 


If  we  let  (^n+])n!a  0,  then  the  kth  coefficient  is  given  by 


CKi  -  k'c>‘(„-l)  +  C<k-U(„-1) 


FOURTH  DERIVATIVE  OF  e 


1  2  Au+G 
-'pj  -e 


U  a  e 


TTIV 

U  =  e 


1  2  Au+G 
-ju  -e 


F  -  e 


-e 


Au+G 


1  2 


0-6u2+-u4)F  +  (12u-Au3)Fa)  +  C-6+-6u2)Ff2)  +•  (-Au)F(3)+-  F(4) 


UIV  -  U  / (3-6u2+u4)  +■ 

l 


(_eAu+^) [A(12u-Au3)  +■  A2(-6+-6u2)  +•  A3(-Au)  +-  A4] 


(+eZ'Au+G))[ 
(_e3(Au+G)j j 

(+eA(Au+G))[ 


A2(-6+-6u2)  +  3A3(-Au)  +-  7A4] 
A3(-Au)  +>  6A4] 
A4]' 


Table  A-3.1  Expressions  for  UIV  arid 


A  32 


Normal-Smallest  Extreme  Value 


Curve  Approximately  Normal 


Region  Integrated 
Using  Simpson's  Rule 


Curve  gives  approximate¬ 
ly  zero  contribution  to 
integral  above  this 
point 


»u  Axis 


PR (Failure) 


S1  v  us2 

1. -(total  area  under  curve) 


Normal-Largest  Extreme  Value 


PR(Feilure)  ■  (Total  area  under  curve) 


Plgure  A-3.1  Limits  of  Integration 


iv)  Being  applicable  to  any  range  of  A  and  G  simply  by  changing  the 
Simpson's  integration  step-size. 


A-3.4  SIMPLIFYING  APPROXIMATIONS 


The  functions  comprising  the  integrand  of  (9)  of  Section  A-3.2  are 
quite  well  behaved  and  bounded  between  0  and  1.  Their  exponential  form  leads 
to  a  number  of  important  approximations.  Let  U  equal  the  integrand  of  (9)  and 
note  that  U  is  the  product  of  the  normal  curve  and  a  "double  exponential." 


The  integral  of  the  normal  curve  is  equal  to  1  within  .5734  x  10 
over  the  interval  -5  to  +5.  Therefore,  since 

Au+G  -  -  Au+G  .  -  -^u^ 

0  *  e  6  M  implies  e  )  e  e  4  •  e  ) , 


Au+G 

Let  F  -  e"*e 


to  at  least  6  places. 


(10) 


If  the  function  F  differs  from  one  over  some  interval  [uj_,  u£]  by 
less  than  5  x  10“5f  then  the  integral  (10)  can  be  evaluated  to  four  places 
simply  by  integrating  the  normal  curve  (actually  the  FORTRAN  error  function 
subroutine  was  used).  That  is, 


)(1-. 00005) du 


1  2 
V 


du 


worst  error  »  .0005  x 


7?re 


du  ■  .0005 


where  the 


The  combination  of  u,  A,  and  G  which  causes  F  to  be  within  5  x  10  5 
of  one  is  easily  determined: 

Au+G 


-e 

e 

> 

.99995 

for 

Au+G 

0 

< 

5  x  10’5 

Take 

au-k; 

o 

< 

1  x  10'5 

Then 

Au+G 

< 

-11.5  j 

(ID 

If  the  function  F  differs  from  sero  over  the  Interval  [-5f+5j  by  less 
than  5  x  10"3,  then  the  lntegrel  can  be  approximated  as  sero  to  four  places. 
The  above  argument  concerning  the  error  applies  here  also;  that  is,  the 
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t*'.*s>*»*i  no<»tj»'>r««>1- 


yi’t&K 


integral  differs  from  zero  by  less  than  5  x  10“5  when  F  differs  from  zero  by 
less  than  5  x  10“5.  it  will  prove  useful  to  know  the  value  of  the  upper  expo¬ 
nent  for  which  this  occurs. 


Au+G 

e"e  <  .00005 

Au+G  . 

Take  e-e  <  1  x  10”5 


Au+G  , ,  c 
-e  <  -11.5 


Au+G 


>  In  11.5 


Au+G 


>  2.44 


(12) 


A-3.5  LIMITS  OF  INTEGRATION 


Equation  (10)  defined  the  limits  of  integration  for  our  problem  as  -5 
to  +-5.  The  discussion  below  shows  that  it  is  possible  to  further  Isolate  the 
interval  on  the  u-axis  in  which  numerical  integration  of  (9)  is  necessary. 

A-3,5.1  Limits  for  Normal-Smallest  Extreme  Value  Case 

For  a  given  A  and  G,  expression  (11)  defines  a  u-value,  call  it  us. , 
below  which  the  integrand  can  be  approximated  as  a  normal  curve. 


U8l 


-G-11.5 


Likewise  (12)  defines  a  u,  say  u^,  above  which  the  contribution  to  the  inte¬ 
gral  is  zero. 


U®2  " 


-G+-2.44 

A 


Therefore,  numerical  Integration  need  only  be  used  on  [u8.,  u82].  This  ie  il¬ 
lustrated  in  Figure  A-3.1.  In  some  Instances  u8.  and  u82  may  be  outside  the 
-5  to  +-5  limits  established  by  (10).  The  limits  should  then  be  taken  as  -5  to 
+■5.  Note  that: 

When  u8^  >  +-5  the  Integrand  can  be  approximated  to  sufficient 
accuracy  by  considering  it  to  be  the  normal  density  function. 

When  us2  <  -5  the  integrand  can  be  approximated  by  considering 
it  equal  to  zero. 

A-3.5. 2  Limits  for  Normal-Largest  Extreme  Value  Case 

The  results  of  the  previous  article  apply  here  directly  except  that 
for  this  case  A  ■  -o  (recall  that  a  Is  always  positive)  and  the  senses  of  both 


'•  • 


m 


.Vgj&p> 
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f 


(il)  and  (12)  are  reversed.  This  makes 

-G+2.44 

uLi - - 


and  for  u  <  the  integral  is  "zero". 
-G-11.5 


also. 


UL2 


and  for  u  >  uj^  the  curve  is  "normal". 

When  uj^  >  +5  the  entire  curve  is  "zero". 
When  UL2  <  -5  the  entire  curve  is  "normal". 


Of  course,  in  those  situations  where  uj.  <  -5,  the  lower  limit  should  be  taken 
as  -5  and  where  uj^  >  +5,  the  upper  limit  becomes  +5. 

A-3.6  DETERMINATION  OF  SIMPSON'S  RULE  STEP-SIZE 


The  truncation  error  for  Simpson's  integration  scheme  is  given  by 


| Error |  -  UIV(£)(Au)4,  Uj  s  5  ^  «2  .  (13) 


The  fourth  derivative  in  this  expression  should  be  evaluated  at  the  £  which 
makes  UIV  a  maximum.  The  complexity  of  UIV  (shown  in  Table  A-3.l)  made  a  com¬ 
plete  analysis  of  its  behavior  impossible.  Study  of  the  form  of  UIV  together 
with  careful  scrutiny  of  a  large  number  of  computed  values  gave  what  was  be¬ 
lieved  to  be  good  estimates  for  the  maximum  of  U*v  for  all  combinations  of  A 
and  G.  These  values  were  substituted  into  (13)  and  Au's  were  determined 
assuming  a  desired  error  of  less  than  5  x  10-5, 

Actually,  only  three  Au's  were  used:  .1,  .05,  and  .01.  Although 
the  decision  of  which  to  use  for  a  given  set  of  the  parameters  was  based  upon 
inspection  of  calculated  values  of  UlV,  an  understanding  of  the  behavior  of 
UIV*  with  changes  in  A  and  G  is  Important.  The  double  exponential  coeffi¬ 
cient  of  UIV  Win  damp  out  any  influence  by  the  other  exponentials  to  make  UIV 
large.  For  A-+0,  UIV  will  be  small.  In  fact,  for  A  small,  U*v  should  behave 
similarly  to  the  fourth  derivative  of  the  normal  curve  since  G  only  acts  as  a 
"scaling  factor"  and  does  not  effect  the  "curvature"  of  U.  Tables  show  the 
fourth  derivative  of  the  normal  curve  to  be  approximately  1. 

As  A+0,  U1^  *  1,  and 

1 5  x  10— 5 1  -  (l)(Au)4 

(Au)4  -  9  x  10-4 

Au  a  .17 
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Second.  Quadrant : 
Normal-Largest  Extreme  Value 


Note:  is  only  given  for  the  interval  where 

Simpson' 8  integration  was  used. 


Figure  A-3.2  for  Various  A  and  G  and  Simpson's  Integration  Step-Size 
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For  large  |A)'s  (neglect  G)  the  U-^mnv  will  behave  as  A4, 
gest  A  of  interest,  A  -  10.0; 

,-4 


Then  for  the  lar- 


(4u>4  -  9  *  10 


104 


-  9  x  10 


-8 


Au  «  .0017 


,IV 


It  can  also  be  seen  from  the  form  of  l)'-'  that  for  large  | G { * s  the  value  of  UIV 
will  be  driven  to  zero  by  one  of  the  exponentials.  (Note  also  that  the  beha¬ 
vior  of  UlV  is  only  of  interest  on  that  interval  where  Simpson's  integration 
is  to  be  used.) 

Recall  that  for  each  of  the  two  cases  being  studied  and  for  each 
pair  of  the  parameters  «  and  y,  the  integration  of  (5)  or  (6)  of  A-3.1  requires 
use  of  Simpson's  integration  on  an  interval  identified  by  Section  A-3.5.  The 
error  produced  by  the  use  of  Simpson's  rule  depends  on  the  derivative  UIV  and 
the  fourth  power  of  the  step-size.  Figure  A-3,2  gives  the  value  of  U*Vmax  at 
various  points  in  the  A-G  plane  and  shows  the  regions  in  the  plane  where  the 
three  step-sizes  (.01,  .05,  and  .1)  could  be  used.  In  the  A-G  plane,  the 
second  quadrant  is  the  region  of  interest  for  the  Normal-Largest  Extreme  Value 
case  (A  ■  -a,  .001  s  a  *  10.;  G  *»  -y,  -40.  *  y  *  0).  The  fourth  quadrant  is 
the  A-G  region  of  interest  in  the  Normal-Smallest  Extreme  Value  case.  In 
Figure  A-3.2,  choose  an  abscissa  and  an  ordinate;  the  number  indicated  is  the 
maximum  of  U*v  found  by  calculating  U  for  many  u's  on  the  interval  in  which 
Simpson's  rule  was  used.  The  cross-hatching  defines  the  regions  where  the 
Indicated  step-size  will  maintain  the  error  less  than  5  x  10_5.  (These  re¬ 
gions  were  delineated  more  carefully  for  the  actual  computer  runs.) 

Example : 

Assume  that  we  wish  to  find  the  probability  of  failure  for  the  Normal 
Largest  Extreme  Value  case  for  a  ■  .7  and  y  -  -5.  For  this  case,  the  para¬ 
meters  A  and  G  in  the  integral  (  9  )  of  A-3.2  become  A  -  -a  —  -w 7  and  G  ■  -y.  ■ 
+5.  Figure  A-3.2  shows  that  for  A  ■  -.7  and  G  ■  +5,  UIV,nax  s  2.0  x  10^,  The 
cross-hatching  shows  that  a  step-size  of  Au  -  .01  is  acceptable  for  the  inte¬ 
gration. 


A-3.7  ERROR  ANALYSIS 

The  possible  errors  produced  in  generating  the  tabulated  values  are 
listed  below: 


Error  Source 

Truncation  of  Integral  on  [-5, +51 
Approximations  of  Section  A-3.5 
Simpson's  integration  truncation 
Computer  round-off 
Round-off  to  4  places 

-A 

The  tabular  values  are  correct  to  2  x  10 


Maximum  Error 

i-7 


-6  x  10 
-1  x  10 
*5  x  10 
*1  x  10 
*1  x  10' 


”5 

-5 

“7 

-4 


A-3.8  SUMMARY 


It  was  shown  that  calculation  of  the  probabilities  of  failure  re¬ 
duces  to  the  numerical  evaluation  of  a  single  integral.  The  Integral  was  com¬ 
puted  by  using  Simpson's  rule  and  by  employing  several  approximations. 
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APPENDIX  4  FLOW  CHART  OF  THE  COMPUTER  PROGRAM  USED  FOR  THE  DETERMINATION  OF 
THE  FATIGUE  STRENGTH  DISTRIBUTION 

A-4.1  EXPLANATION  OF  THE  SYMBOLS  USED  IN  THE  FLOW  CHART 

Ci,  C2,  C3,  Ci,  ...  Cm  are  the  lives  at  which  the  distributions  of 

fatigue  strength  are  to  be  determined.  In  the  present  study,  these 
are  103,  104,  5xl04,  105,  5xl05,  106,  5xl06,  107,  and  108  cycles. 

m  is  the  total  number  of  lives  which  in  this  study  is  equal  to  9. 

U  is  /the  code  number  for  a  mode  of  stress,  for  example; 

U  *  0  is  for  the  alternating  (completely  reversed)  stress, 

U  =  1  is  for  the  maximum  stress, 

U  -  2  is  for  the  range  of  stress  (maximum-minimum). 

US  is  the  identification  number  for  the  units  of  stress,  for  example; 

US  -  1  is  for  the  units  in  psi, 

US  ■  2  is  for  the  units  in  ksi, 

US  *  3  is  for  the  units  in  short  tons/in^, 

US  ■  4  is  for  the  units  in  kilograms/mm^ , 

US  -  5  is  for  the  units  in  long  tons/in^. 

UL  is  for  the  identification  number  for  the  units  of  life,  for  example; 
UL  -  1  is  for  the  units  in  cycles, 

UL  *  2  is  for  the  units  in  reversals  (two  reversals  ■  one  cycle), 

UL  ■  3  is  for  the  units  in  hours  where  frequency  (cycles/hr)  was 
known. 

Sn,  is  the  mean  stress  which  in  the  present  study  was  always  equal  to 
zero. 

t 

Si,  S2,  S3,  ...  Si,  ...  Sm  are  the  values  of  the  applied  stresses,  and, 

Ll>  l2 *  L3*  •••  Li»  •••  Im  are  the  values  of  the  life  to  failure  corres¬ 

ponding  to  the  above  stresses. 

1  is  the  subscript  to  designate  the  ith  data  point  (for  example,  Li 

and  Si  are  the  coordinates  of  the  1th  data  point  of  the  S-N  type 

data) . 

n  is  the  number  of  data  points. 

MRn,i,  2*  •••»  ^®n  n  are  the  median  rank  values  for  the  sample 
size  of’n(for  example,  Ml^i,  MR5  2,  ^5,3,  MR5  4,  and  MR5  5  are  the 
median  rank  values  when  n  »  5).  *  "  “ 

Zn,i,  Zn  3*  •••»  Zn,n  are  the  standardized  normal  variates  corres¬ 

ponding  to  the  median  ranks  of  sample  size  n.  (These  are  also  called 
the  median  ranks  z  values  as  given  in  Table  6.1  on  Page  49). 
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YV±  is  the  linear  mode  variable  of  the  ordinate  scale  for  Weibull  dis¬ 
tribution  for  the  ith  data  point,  and  it  is  also  the  same  for 
Smallest  Extreme  Value  distribution  (See  Equations  6.24  and  6.25 
on  Pages  40  and  45). 

YEi  is  the  linear  mode  variable  of  the  ordinate  scale  fcr  Largest  Ex¬ 
treme  Value  distribution  for  the  ith  data  point. 

YLi  is  the  linear  mode  variable  of  the  ordinate  scale  for  Logistic  dis¬ 
tribution  for  the  i6"  data  point. 

YNi  is  the  linear  mode  variable  of  the  ordinate  scale  for  Normal  distri¬ 
bution  for  the  ith  data  point. 

a  is  the  intercept  of  the  least  squares  line  on  the  Y-axis.  (See 
Equation  6.3  in  Section  6.1.1  on  Page  28). 

cN,  aL,  cM,  oE,  and  oN  are  the  intercepts  of  the  least  squares  lines 
fitted  to  Normal,  Logistic,  Largest  Extreme  Value,  Smallest  Extreme 
Value  and  Weibull  distributions  respectively. 

8  is  the  slope  of  the  least  squares  line.  (See  Equation  6.2  in  Section 
6.1.1  on  Page  28.) 

BN,  BL,  BM,  BE,  and  BW  are  the  slopes  of  the  least  squares  line  fitted 
to  Normal,  Logistic,  Largest  Extreme  Value,  Smallest  Extreme  Value, 
and  Weibull  distribution  respectively. 

R  is  the  correlation  coefficient.  (See  Equation  6.6  in  Section  6.1.3 
on  Page  29) 

RN,  RL,  RM,  RE,  and  RW  are  the  correlation  coefficients  of  Normal,  Lo¬ 
gistic,  Largest  Extreme  Value,  Smallest  Extreme  Value,  and  Weibull 
distributions  respectively, 

1C  is  an  identification  number  for  each  data  set. 

LCj  Is  the  logarithm  of  xltn  j  at  which  the  fatigue  strength  distribu¬ 
tion  is  desired. 

LSTjti  is  the  logarithm  of  fatigue  strength  at  a  given  life  J  for  a 
given  date  point  1. 

SN  is  the  scatter  of  fatigue  strength  at  a  given  life. 

DEV  is  the  standard  deviation. 

DX  and  DR  are  the  differences  between  the  two  consecutive  assumed 
values  of  Xq  for  the  Weibull  case. 

ADX  is  the  ahaolute  value  of  DX. 


DR  Is  the  difference  between  the  two  consecutive  values.  of  correla¬ 
tion  coefficients  corresponding  to  the  two  assumed  values  of  Xq 
for  the  Weibull  case. 

A-4.2  EXPLANATION  OF  THE  FLOW  CHART. 

The  flow  chart  of  the  computer  program  used  to  determine  the  fa¬ 
tigue  strength  distribution  was  made  up  of  three  parts.  In  the 
first  part  of  this  program,  the  scatter  in  the  life  was  converted 
into  the  scatter  in  the  fatigue  strength,  and  the  y-ordinate  of 
non-linear  mode  variables  into  the  linear  mode  variables  for  each 
distribution  (For  example,  JP  was  converted  to  In  In  y—p  for  the 
Weibull  case).  In  the  second  part,  the  distribution  —  functions 
(Weibull,  Normal,  Logistic,  Largest  Extreme  Value  and  Smallest 
Extreme  Value)  were  fitted  to  these  strengths  data  simultaneously, 
and  the  one  with  the  highest  degree  of  fit  was  the  best  fitting 
distribution.  The  part  three  is  an  external  Function  Fit  which 
was  used  in  part  one  and  part  two  to  fit  the  least  squares  lines 
to  the  data.  As  shown  in  Section  A-4.3,  blocks  were  made  around 
different  sections  of  each  of  these  three  parts.  The  specific 
functions  of  each  of  these  blocks  in  the  total  program  are  discus¬ 
sed  in  the  following  sections. 

A-4.2.1  Part  One  -  Conversion  of  Life  Data  to  Strength  Data  and 
Non-Linear  Mode  Variables  to  Linear  Kode  Variables. 

Block  JL:  The  function  of  this  block  is  to  store  all  the 
necessary  information  and  the  various  symbols  used  in 
this  computer  program.  The  program  will  be  rejected  in 
either  case  when  number  of  data  points  (N)  is  larger  than 
40  or  number  of  lives  (m),  at  which  the  strength  distrib¬ 
ution  is  determined,  is  larger  than  11.  (See  the  flow 
chart  in  Section  A-4.3). 

Block  2t  for  soma  data,  where  the  units  of  stress  (US) 
values~are  In  tona/in2,  killograma/ln2  etc.,  and  the  modes 
of  stress  (U)  In  terms  of  range  of  stress,  etc.,  these 
are  converted  to  the  units  In  psl  and  the  mode  of  stress 
in  completely  reversed  stress.  This  is  done  for  each 
data  point  in  a  given  data  set. 

Block  2:  Non-linear  mode  of  y  co-ordinate  of  each  data 
point  (in  this  case  the  median  rank  value  (£)  assigned 
to  each  data  point)  is  converted  into  the  linear  mode 
of  y  co-ordinate  (e.g.  In  In  for  the  Weibull  case). 
This  is  done  for  all  the  four  **  distribution  functions. 

A*  In  some  data  if  the  units  of  life  are  in  terma 
of,  say,  hours  or  number  of  reversals,  these  are  convert¬ 
ed  into  cycles.  This  is  done  for  each  data  point.  The 
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values  of  stress  and  life  are  converted  into  the  loga¬ 
rithm  of  stress  (log  S)  and  logarithm  of  life  (log  F) . 

Block  5_:  Using  the  External  Function  Fit  (as  discussed 
in  Section  A-4.2.3  and  the  flow  chart  shown  in  A-4.3.3) 

1  least  squares  line  (Y  =  a  +  6][)  is  fitted  through  these 
data  points  log  (S)  and  log  (N)  found  in  Block  4^  and 
the  values  of  a  and  8  is  determined  in  terms  of  log  a 
and  log  3. 

Block  j5:  Using  the  values  of  e  and  B  as  found  by  Block 
5,  and  the  procedure  as  discussed  in  Section  6.2,  page 
31,  the  scatter  in  strength  in  terms  of  logarithm  of 
strength  (LSTj^)  at  a  given  life  j  is  determined. 

A-4 .2.2  Part  Two  -  Determination  of  the  Best  Fitting  Distribution. 


Block  7 :  The  values  of  logarithm  of  strengths  [log  (s)] 
at  a  given  life  (as  found  in  Block  6)  are  arranged  in  an 
increasing  order,  and  they  are  then  converted  from  log 
(S)  to  S. 

Block  jh  Using  the  strength  values  as  found  from  Block 
7  and  the  information  from  Block  3,  the  distributions 
such  as  Normal,  Logistic,  Largest  Extreme  Value  and  Small¬ 
est  Extreme  Value  are  fitted  (as  discussed  in  Section 
6.5)  by  aid  of  the  External  Function  Fit  discussed  in 
Section  A-4.2.3.  The  values  of  the  correlation  coef¬ 
ficients,  R,  for  each  of  these  distribution  are  then 
computed. 

Block  9:  As  discussed  in  Section  6.5.1  (on  Weibull 
distribution,)  it  is  necessary  to  determine  the  value 
of  Xjo,  the  lower  bound  of  strength,  before  fitting  the 
Weibull  distribution  to  the  data.  In  this  block,  the 
initial  value  of  ^  is  assumed  in  the  manner  shown. 

Block  10;  Other  values  of  XQ  subsequent  to  the  one  as¬ 
sumes  in  Block  9  are  assumed,  and  the  magnitude  of  the 
new  values  of  Xo  depends  on  the  magnitude  of  the  pre¬ 
vious  ones,  as  shown  in  this  block. 

Block  11;  Once _a  value  of  Xq  is  assumed,  the  values  of 
Xj  ■  In  (xi  -  Xq)  -  In  (SN^  -  Xq)  are  found.  Using  the 
value  of  Xj.  and  YW|  (as  found  from  Block  3).  The  Weibull 
distribution  is  fitted  to  these  data  by  External  Func¬ 
tion  Fit  and  the  corresponding  values  of  the  parameters 
a  and  8  of  the  least  squares  line  (Y  ■  a  +  8X)  and  the 
correlation  coefficient,  R,  are  computed.  This  is 
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repeated  twelve  times  and  twelve  values  of  3^  and  their 
corresponding  values  of  a,  3,  and  R  are  computed.  (Note: 
The  number  of  repetitions,  twelve,  was  chosen  arbitrarily). 

Block  12:  Out  of  these  twelve  values  of  the  correlation 
coefficients  which  has  the  maximum  value  (the  best  fit) 
is  selected,  and  its  corresponding  value  of  3^  is  taken 
as  the  "true"  value  cf  the  lower  bound  of  strength  3T0. 

Block  13:  The  fatigue  strength  distributions  (such  as, 
Weibull,  Normal,  Logistic,  Largest  Extreme  Value  and 
Smallest  Extreme  Value),  the  values  of  their  parameters 
(in  units  of  KSi  where  appropriate)  and  their  correspond¬ 
ing  values  of  correlation  coefficients  are  printed  out. 

From  this  information  computer  also  prints  out  the  name 
of  the  best  fitting  distribution. 

A_4*2*3  Part  Three  -  External  Function  Fit 

Block  14 :  External  Function  Fit  is  used  in  Part-one  and 
Part-two  of_the  flow  chart  to  fit  a  least  squares  line 
(Y  -  a  +  BX)  to  a  set  of  data  points  where  the  values 
of  o  and  6  are  computed  using  Equations  (6.2)  and  (6.3) 
given  on  page  2S.  The  corresponding  value  of  the  cor¬ 
relation  coefficient  is  then  calculated  using  equations 
(6.6),  (6.7),  and  (6.8)  as  given  on  page  29. 
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